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The  report  "Predictability  of  LNG  Vapor  Dispersion 
from  Catastrophic  Spills  onto  water:  An  Assessment" 

was  prepared  for  the  Coast  Guard  by  Dr.  J.  A.  Havens, 
Associate  Professor  of  Chemical  Engineering,  University 
of  Arkansas,  while  on  sabbatical  leave  acting  as  Technical 
Advisor (-  Cargo  and  Hazardous  Materials  Division. 

The  report  reviews  the  methods  which  have  been  used 
to  predict  tlie  dispersion  of  flammable  gases  from  a very 
large  spill  of  liquefied  natural  gas  (LNG)  onto  water. 

It  provides  &■  comprehensive  description  of  these  methods 
and  contains  an  assessment  of  the  present  "state  of  the 
art"  in  predictability  of  vapor  dispersion  from  such 
spills.  The  report  recommends  further  work  which  should 
be  done  to  increase  confidence  in  such  predictions. 

This  report  should  help  to  answer  a number  of  ques- 
tions which  «*re  frequently  asked  by  persons  interested  in 
water  transportation  of  this  important  cargo,  and  provides 
a recommended  course  of  action  which  should  provide  answers 
to  questions  which  remain  regarding  the  dispersion  of 
flammable  vapors  from  large  spills  of  hazardous  materials. 
Comments  and  recommendations  are  solicited  and  should  be 
sent  to  the  Chief,  Cargo  and  Hazardous  Materials  Division 
(G-MHM/83) , U.  S.  Coast  Guard,  Washington,  D.C.  20590. 
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FOREWORD 


This  review  of  mathematical  models  which  have  been  used 
to  predict  the  downwind  travel  of  flammable  gas  mixtures  in 
the  event  of  a catastrophic  spill  of  liquefied  natural  gas 
onto  water  was  undertaken  while  the  author  was  on  sabbatical 
leave  from  the  Department  of  Chemical  Engineering,  University 
of  Arkansas,  serving  as  Technical  Advisor,  Cargo  and 
Hazardous  Materials  Division,  Office  of  Merchant  Marine  Safety, 
U.S.  Coast  Guard  Headquarters,  Washington,  D.C. 

The  motivation  for  this  review  resulted  from  two  needs 
of  the  Coast  Guard: 

1.  The  Coast  Guard  is  actively  developing,  through 
contract  research  and  in-house  efforts,  techniques  for 
the  assessment  of  hazards  associated  with  the  marine 
transportation  of  chemicals.  A significant  part  of  the 
hazardous  nature  of  some  chemicals  shipped  by  water 
relates  to  fire  and  explosion  behavior.  The  increasingly 
routine  marine  carriage  of  volatile  flammable  liquids 
and  liquefied  flammable  gases  in  large  quantities  carries 
with  it  the  risk  of  fire  and  explosion  phenomena  result- 
ing from  formation  of  large  flammable  vapor  clouds  in 
the  event  of  an  accident.  The  assessment  of  such  risks 
and  the  development  of  emergency  response  procedures 
requires  a methodology  for  predicting  the  extent  and 
nature  of  flammable  cloud  formation  in  a variety  of 
possible  accident  scenarios.  Thus  the  Coast  Guard  has 
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a general  need  for  accurate  vapor  dispersion  models. 


2.  The  proposed  large  scale  importation  of  liquefied 
natural  gas  (LNG)  into  the  United  States  is  the  subject 
of  intense  argument,  particularly  in  relation  to  the 
assessment  of  risk,  to  the  public  from  accidental  release 
of  LNG  as  a result  of  ship  collision.  In  the  event  of 
a catastrophic  release  it  iSv^onsidered  highly  likely 
that  an  immediate  fire  would  ensufr.  However,  in  the 
event  that  ignition  did  not  occur  immediately,  an  LNG 
vapor  cloud  would  form  over  and  downwind  of  the  spill 
site.  Wide  disagreement  regarding  the  extent  of  travel 
(and  the  accompanying  possible  public  exposure)  of  the 
flammable  portion  of  such  a cloud  has  contributed  to  an 
apparently  growing  concern  regarding  the  risks  associated 
with  LNG  importation.  The  Coast  Guard  is  responsible  for 
regulating  the  movement  of  LNG  by  water  in  the  United 
States  and  thus  has  a specific  interest  in  the  development 
<f  accurate  LNG  vapor  dispersion  models.  There  appeared 
to  be  a need  for  a review  and  assessment  of  vapor  disper- 
sion predictability  by  someone  not  immediately  involved 
in  LNG  safety  related  research. 

Since  the  Coast  Guard's  primary  interest  is  in  LNG  spills 
on  water,  this  review  was  immediately  restricted.  Several  models 
for  LNG  vapor  dispersion  which  have  been  used  primarily  for 
analyzing  vapor  dispersion  from  land  spills  are  essentially 
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identical  to  those  reviewed  herein.  However,  if  they  had 
not  been  used  as  a basis  for  published  predictions  of  vapor 


dispersion  from  LNG  spills  on  water  they  were  not  included. 
Furthermore,  the  scope  of  this  review  was  limited  to  the 
predictability  of  dispersion  from  a very  large  ING  spill  on 
water.  No  consideration  was  given  to  site-specific  factors 
which  may  have  an  important  bearing  on  the  assessment  of 
downwind  flammable  cloud  travel,  such  as  topographical  features 
and  structures.  Likewise,  no  consideration  was  given  to  the 
specific  applicability  of  weather  conditions,  since  this  would 
depend  on  the  site  involved  as  well  as  the  traffic  control 
measures  which  are  imposed.  For  example,  if  LNG  ship  movement 
is  restricted  to  daylight  hours,  the  probability  of  a very 
large  spill  during  stable  or  inversion  conditions  may  be 
remote  for  some  ports.  My  intent  was  to  review  published  models 
used  to  predict  downwind  travel  for  a very  large  spill  (25,000  M^) 
and  to  identify  and  explain  the  differences  in  those  models. 

I have  also  offered  recommendations  for  future  work  based  on 
the  assessment  of  the  models  reviewed.  There  may  also  be  other 
models  proposed  for  the  prediction  of  vapor  dispersion  from 
LNG  spills  on  water  which  I have  overlooked.  If  this  is  the 
case,  such  omission  is  due  only  to  my  time  constraints  for 
reviewing  the  literature. 

This  work  could  not  have  been  completed  without  the 
excellent  cooperation  received  from  all  of  the  parties  whose 
work  was  reviewed.  At  my  request,  all  of  the  groups  clarified 
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questions  which  I had  based  on  my  review  of  the  published 
models,  and  one  group  iCabot  Corporation)  provided  a computer 
program  which  I required  to  make  predictions  utilizing  their 
model.  However,  if  errors  appear  in  the  review  they  are  my 
own  produce. 

The  reader  should  be  aware  that  the  presentation  of 
these  models  gives  no  insight  into  the  historical  perspective 
in  which  they  might  be  viewed.  A close  look  at  the  literature 
cited  in  this  report  indicates  that  there  was  indeed  a 
"development"  process  involved  in  the  formulation  of  these 
models  for  LNG  vapor  dispersion.  It  is  not  surprisi  g that 
the  models  which  are  recommended  for  further  use  and  evaluation 
are  in  a real  sense  the  product  of  efforts  to  modify  or  build 
on  the  efforts  of  the  earlier  investigators  in  the  field. 

In  order  to  insure  accuracy  of  description  and  inter- 
pretation of  the  models  reviewed  herein,  a draft  of  this  report 
was  sent  to  all  parties  whose  work  is  discussed  v* rh  a request 
that  they  examine  the  description  of  their  model  for  technical 
and  interpretive  accuracy.  Comments  were  received  from  all  cf 
the  groups  and  were  carefully  considered  in  the  preparation  of 
the  final  report.  Corrections  and  revisions  of  the  draft 
report  were  made  in  several  instances  as  a result  of  the 
comments  received.  For  the  sake  of  completeness,  the  comments 
on  the  draft  report  are  appended. 


J.  A.  HAVENS 
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III.  EXECUTIVE  SUMMARY 

A number  of  predictions  of  LNG  vapor  cloud  formation  and 
dispersion  which  might  result  from  a catastrophic  LNG  spill 
onto  water  have  been  published.  The  predictions  of  the 
following  groups  have  been  repeatedly  cited  in  the  literature 
related  to  safety  of  marine  LNG  transportation: 

1.  U.S.  Bureau  of  Mines  - Burgess  et  al.  (1,  2) 

2.  American  Petroleum  Institute  - Feldbauer  et  al.  (3) 

3.  Cabot  Corporation  - Germeles  and  Drake  (4) 

4.  U.S.  Coast  Guard  CHRIS  (Chemical  Hazard  Response 
Information  System)  - Arthur  D.  Little,  Inc.  (5) 

5.  Professor  James  Fay,  Massachusetts  Institute  of 
Technology  (6) 

6.  Federal  Power  Commission  (7) 

7.  Science  Applications,  Inc.  (8) 

Crder  of  magnitude  differences  in  the  predictions , based  on 
these  models,  of  the  extent  of  flammable  vapor/air  mixtures 
following  a catastrophic  spill  are  significant  in  the  overall 
assessment  of  the  potential  risk  of  marine  transportation  of 
LNG. 

i'*  The  purpose  of  this  study  is:  fivefold : 

> ■ 1)  To  provide  a detailed  description  of  the  mathematical 

models  upon  which  published  predictions  of  LNG  vapor  travel 

downwind  of  catastrophic  LNG  spills  onto  water  have  been 

based.  -• tr  ‘ Ck  ■ ) 
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i 2)  To  estimate,  using  these  models, r the  maximum  downwind 
travel  of  flammable  LNG  vapor/air  mixtures  for  a ^standard* 
spill  scenario,  so  that  ,a  valid  comparison  can  be  made  of 
the  results  obtained  when  thd  different  models  are  used  to 
describe  the  same  event, 

3,  To  identify  the  reason  for  differences  in  predictions 
which  occur  when  the  models  are  used  to  describe  the  same 
event,  and  to  assess  the  technical  credibility  of  the  " ’ 

methodology  which  results  in  such  differences. 

,4,  To  define  the  present  “"state  of  the  art"  in  predicta- 
bility of  LNG  vapor  dispersion  from  catastrophic  spills 
onto  water,  with  emphasis  on  the  extent  to  which  the 
present  state  of  the  art  justifies  reliance  on  existing 
published  predictions  in  formulating  LNG  safety  management 
programs.  ^v. 

5.  To  provide  recommendations  for  further  work  which 
would  increase  confidence  in  the  predictability  of  vapor 
dispersion  from  catastrophic  LNG  spills  onto  water. 

The  models  used  by  the  groups  cited  above  for  prediction 
of  vapor  cloud  formation  and  dispersion  can  be  categorized  as 
follows : 

1.  Predictions  which  utilize  classical  air  pollutant 
dispersion  models  which  were  developed  to  describe  relative- 
ly near-field  dispersion  of  neutrally  buoyant  materials. 
These  models  are  based  on  the  general  observation  that  the 
concentration  profiles  downwind  of  a pollutant  source  can 


be  represented  by  a Gaussian  or  Normal  distribution. 

This  model  type  is  subdivided  to  describe  two  different 
dispersion  phenomena. 

a.  Dispersion  of  an  essentially  instantaneous 
release  of  a pollutant  into  the  atmosphere,  the 
dispersion  being  associated  with  the  growth  of  this 
instantaneously  released  "puff",  or  cloud,  as  it 

is  being  translated  by  the  wind.  The  predictions 
due  to  Fay,  Germeles  and  Drake,  and  CHRIS  utilize 
this  type  of  model. 

b.  Dispersion  of  material  which  is  being  emitted 
at  a continuous  steady  rate  forming  a "plume"  down- 
wind of  the  emission  source.  The  predictions  of 
Burgess,  Feldbauer,  and  the  FPC  utilize  this  type 
of  model. 

2.  Predictions  based  on  solution  of  the  combined  mass, 
momentum  and  energy  balance  equations.  (The  classical  air 
pollutant  dispersion  equations  of  category  1 above  are  a 
special  case  where  energy  effects  and. momentum  effects  are 
not  considered) . The  SAI  predictions  utilize  this  type 
of  model. 

The  "standard  scenario"  LNG  spill  which  is  assumed  in  this 
report  for  purposes  of  comparison  of  the  above  models  is  an 
.instantaneous  release  of  25,000  cubic  meters  of  LNG  onto  water. 
It  is  considered  that  such  an  event  provides  a conservative 
upper  limit  on  the  severity  of  a spill  which  might  conceivably 
occur . 
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Table  IXI-1  shows  the  maximum  downwind  distance  to  the 
time-average  5%  vapor  concentration  level  following  an  instan- 
taneous 25,000  cubic  meter  spill  as  predicted  by  the  models 
suggested  by  the  seven  groups  above.  The  distances,  with 
the  exception  of  the  estimate  attributed  to  SAI,  were  computed 
by  the  author  using  the  procedure  suggested  by  the  investigating 
groups  cited.  The  corresponding  distance  for  SAI's  model  could 
not  be  computed  due  to  the  proprietary  nature  of  the  SAI  com- 
puter model.  Table  III-l  therefore  includes,  for  comparison, 
the  distance  predicted  by  SAI  for  a 37,000  cubic  meter  spill  as 
described  in  their  risk  assessment  study  prepared  for  the  Western 
LNG  Terminal  Company  (8).  In  reviewing  Table  III-l,  it  should 
be  noted  that  the  meteorological  conditions  suggested  as 
applicable  by  the  groups  are  not  necessarily  the  worst  that 
might  have  been  assumed.  Specifically, 

1.  The  0.75  mile  distance  obtained  with  the  FPC  model 
reflects  the  assumption  of  neutral  atmospheric  stability 
values  (Pasquill  D) , as  recommended  by  the  FPC  staff. 


2.  The  5.2  mile  distance  obtained  with  the  American 
Petroleum  Institute  Model  assumes  vertical  dispersion 
characterized  as  Singer  and  Smith  D and  horizontal 
dispersion  characterized  as  Pasquill  C,  following  the 
procedure  suggested  by  Feldbauer  (J' 

3.  The  SAI  prediction  assumes  neutral  atmospheric 
stability;  however  the  result  is  claimed  to  be  essentially 


TABLE  III-l 


MAXIMUM  DOWNWIND  DISTANCE  TO  TIME-AVERAGE  5% 
CONCENTRATION  LEVEL  FOLLOWING  25,000  M^ 
INSTANTANEOUS  SPILL  OF  LNG  ONTO  WATER 


(Assumes  5 MPH  Wind  except  as  noted  and  Meteoroloaical 
Conditions  Considered  Applicable  by  Investigating  Groups) 


Model 

U.  S.  Bureau  of  Mines  (1,2) 
American  Petroleum  Institute  (3) 
Cabot  Corporation  (4) 

U.  S.  Coast  Guard  CHRIS  (5) 
Professor  James  Fay  (6) 

Federal  Power  Commission  (7) 
Science  Applications,  Inc.  (8) 


Distance  (MILES) 
25.2  - 50.3* 
5.2 
11.5 
16.3** 

17.4** 

0.75 

1.2*** 


* A range  was  presented  to  indicate  uncertainty  x.i  vapor 
evolution  rate 

**  Wind  velocity  not  considered  explicitly  in  model 

For  37,500  cubic  meter  instantaneous  release,  wind  velocity  = 
6.7  MPH 


*** 


the  same  when  stable  atmospheric  conditions  are 
assumed.  The  SAI  model  also  gives  longer  downwind 
travel  distances  when  higher  wind  velocities  are 
assumed, in  contrast  to  the  other  models. 

4.  The  predictions  obtained  using  the  Bureau  of 
Mines  (Burgess)  model,  the  CHRIS  model.  Fay's  model, 
and  the  Cabot  (Germeles  and  Drake)  model  assume 
stable  atmospheric  conditions. 

The  variation  in  the  predictions  shown  in  Table  III-l 
is  significant  in  assessing  the  potential  hazard  associated 
with  a large  accidental  release  of  LNG  on  water  if  the 
release  should  occur  without  immediate  ignition  of  the  flam- 
mable vapor  mixture  at  the  spill  site.  The  probability  of 
the  cloud  reaching  the  maximum  distance  at  which  a 5% 
concentration  just  persists  is,  however,  considered  very 
low  due  to  the  anticipated  contact  with  ignition  sources 
which  would  develop  as  a result  of  frictional  heating 
accompanying  such  catastrophic  accidents.  Even  if  the  cloud 
were  not  ignited  at  the  spill  site  it  is  unlikely  that  the 
cloud  would  travel  over  populated  areas,  to  the  extent  pre- 
dicted by  the  models  in  Table  III-l,  without  being  ignited. 
Nevertheless,  the  predictability  of  vapor  dispersion  from  a 
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catastrophic  release  of  LNG  onto  water  does  have  a bearing 
on  the  routing  and  traffic  control  of  vessels  as  well  as 
for  emergency  response  considerations.  Furthermore,  a 
reasonably  accurate  prediction  of  the  dispersion  process  is 
required  for  a characterization  of  cloud  burning  and  to 
assess  potential  damage  which  might  result  from  explosions 
of  vapor/air  mixtures,  if  such  explosions  are  possible. 
Although  the  experience  to  date  indicates  that  detonation 
of  unconfined  LNG  vapor/air  mixtures  is  not  likely,  a good 
method  of  vapor  dispersion  prediction  would  be  v>' 1 uable  in 
attempts  to  understand  the  circumstances  under  which 
detonations  of  vapor/air  mixtures  might  be  expected,  such 
as  partial  confinement  and  high  energy  initiation.* 

Analysis  of  the  models  and  the  results  predicted  for 
dispersion  from  a 25,000  cubic  meter  instantaneous  spill 
indicate  that  all  of  the  variation  in  results  shown  in  Table 
III-l  for  the  classical  air  pollutant  dispersion  models  can 
be  attributed  to  four  factors. 

1.  The  methods  used  to  estimate  the  rate  at  which 

the  vapor  enters  the  atmosphere  from  the  liquid  LNG 

pool  results  in  estimates  thereof  ranging  from  1.43  x 

105  ft^/sec  (FPC)  to  2.0  x 10^  ft^/sec  (at  atmospheric 

* The  Coast  Guard  is  presently  sponsoring  a test  program  at 
the  Naval  Weapons  Center,  China  Lake,  CA.,  to  determine  the 
burning  characteristics  of  large  vapor  clouds  (39) . 
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pressure,  70°F)  (Burgess) . This  factor  is  primarily 
responsible  for  the  very  short  distance  predicted  by 
the  FPC. 

2.  Four  of  the  classical  models  incorporate  effects 
associated  with  the  gravity  induced  spreading  of  the  cold 
LNG  vapors  (FPC,  Feldbauer,  Fay,  Germeles,  and  Drake); 
two  do  not  (Burgess,  CHRIS).  Those  models  which  do 
incorporate  such  effects  assume  a sequential  process; 
spill-pure  vapor  cloud  formation  - gravity  spread  with  or 
without  air  entrainment  - dispersion  by  atmospheric 
turbulence.  The  specific  method  of  treatment  differs 
widely,  and  the  resulting  differences  are  reflected  in 
the  varied  predictions  of  downwind  distance  shown  in 
Table  III-l. 

3.  Some  of  the  models  (Feldbauer,  CHRIS,  FPC, 

Germeles  and  Drake)  incorporate  corrections  for  the 
area  nature  of  the  source  (the  classical  equations 
used  in  all  of  the  models  are  derived  for  a point 
source  emission),  while  others  do  not.  (Burgess,  Fay). 

The  method  of  treatment  of  the  area  nature  of  the  source 
appears  relatively  unimportant  to  the  final  differences 
in  results,  except  for  Feldbauer's  model,  whose  result 
is  strongly  influenced. 

9 

4.  The  predictions  do  not  all  assume  the  same  atmos- 
pheric stability  categories.  Atmospheric  stability 
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considered  applicable  for  such  a prediction  varies  from 
neutral  to  very  i cable,  with  a strong  effect  on  the 
results . 

The  following  conclusions  are  drawn: 

1.  This  review  and  comparison  of  published  predictions  of 
the  downwind  travel  of  flammable  gas-air  mixtures  following 
the  instantaneous  release  of  25,000  cubic  meters  of  LNG 
onto  water  identifies  the  sensitivity  of  such  predictions 
to  the  following  factors. 

a.  Characterization  of  atmospheric  turbulence 
(stability) 

L.  Allowances  for  area-source  effects 

c.  Specification  of  vapor  release  rate 

d.  Allowances  for  gravity  spread/air  entrainment 
effects 

2.  The  shortest  distance  to  the  time  average  5%  concentra- 
tion level  for  a 25,000  cubic  meter  instantaneous  spill 
predicted  by  the  models  reviewed  herein  is  0.75  miles. 

This  distance,  predicted  by  the  FPC  model,  results 
primarily  from  the  use  of  an  unrealistically  low  vapor 
release  rate  and  the  use  of  neutral  atmosphere  stability 
characteristics.  The  FPC  estimate,  in  the  author's  opinion, 
is  not  justified. 
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3.  At  the  other  extreme,  distances  of  the  order  of  tens 

of  miles  are  predicted  for  a 25,000  cubic  meter  instantaneous 
spill  under  stable  weather  conditions  using  continuous  plume 
models  (Burgess)  which  do  not  account  for  any  heat  transfer 
or  momentum  transfer  effects.  Such  estimates  are  not 
justified  in  this  author's  opinion. 

4.  Intermediate  distances  to  the  5%  concentration  level 
are  predicted  for  a 25,000  cubic  meter  spill  during  stable 
weather  conditions  by  Germeles  and  Drake  (11.5  miles). 

Fay  (17.4  miles)  and  the  CHRIS  model  (16.3  miles).  The 
difference  in  predicted  downwind  distances  obtained  with 
the  CHRIS  and  Germeles-Drake  models  can  be  attributed 
primarily  to  the  inclusion  of  gravity  spread/air  entrainment 
effects  in  the  Germeles  and  Drake  Model.  The  rough  agree- 
ment of  these  two  predictions  with  the  value  of  17.4  miles 
predicted  by  Fay  must  be  considered  fortuitous  since  the 
modeling  process  assumed  by  Fay  is  quite  different  from 

the  other  two.  Professor  Fay  now  believes  (42)  that  his 
model  should  be  used  with  different  assumptions  than 
originally  described  by  Lewis  and  Fay,  in  which  case  sub- 
stantially longer  distances  result.  In  the  author's  opinion, 
the  model  of  Germeles  and  Drake  provides  a more  plausible 
estimate  of  the  LNG  dispersion  process  following  a large 
rapid  spill  than  the  Fay  or  CHRIS  models,  since  the  model 
incorporates  a rational,  if  simplified,  description  of  an 
anticipated  gravity  spread  phase.  Further  effort  to 
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improve  this  type  model  as  an  alternative  to  a more 
complex  numerical  procedure  has  merit,  particularly  for 
routine  usage  where  time  and  expe  ise  constraints  are 
important . 

5.  The  estimate,  using  Feldbauer's  model,  of  5.2  miles 
for  the  downwind  distance  to  the  5%  concentration  level 
following  a 25,000  cubic  meter  spill  can  be  attributed 
to  the  predicted  dilution  and  corresponding  extreme 
width  (~2  miles)  of  the  cloud  at  the  end  of  the  gravity 
spread  phase.  Feldbauer's  allowance  for  air  entrainment 
during  the  gravity  spread,  which  involves  the  assumption 
of  a constant  cloud  depth,  is  based  on  observations  of 
small  spills  (10  M^)  and  the  extension  to  very  large  spills 
appears  uncertain.  Further,  representation  of  the  cloud 

at  the  end  of  the  gravity  spread  phase  as  a series  of 
dispersed  point  sources  on  a line  perpendicular  to  the 
direction  of  cloud  travel  is  not  realistic  in  view  of  the 
resulting  prediction  of  shorter  distances  with  increasing 
atmospheric  stability. 

6.  The  primary  reason  for  the  even  shorter  downwind 
distance  (~1  mile)  to  the  5%  concentration  level  predicted 
by  SAI  for  an  even  larger  (37,500  cubic  meters)  spill 
appears  to  be  the  predicted  highly  turbulent  motion  and 
associated  air  entrainment  induced  during  the  gravity 
spread  phase  of  the  cloud. 

7.  In  the  author's  opinion,  the  predicted  maximum  distances 
of  about  5 miles  by  Feldbauer  and  about  1 mile  by  SAI  for 
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flammable  cloud  travel  following  instantaneous  release 
of  25,000  cubic  meters  of  LNG  onto  water  cannot  be 
rationalized  on  the  basis  of  any  argument  thus  far  advanced 
except  that  of  gravity  spread/air  entrainment  effects,  and 
experimental  verification  of  these  effects  has  not  been 
adequately  demonstrated. 

8.  It  was  not  possible  within  the  limits  imposed  by  this 
review  to  evaluate  the  accuracy  of  the  predictions  published 
by  SAI . Rather,  the  author  has  reviewed  the  methodology 
described  by  SAI  and  believes  that  such  techniques  hold 
the  most  promise  for  accurate  prediction  of  vapor  disper- 
sion from  catastrophic  spills  on  water.  A program  designed 
to  evaluate  the  accuracy  of  the  SAI  model  or  other  models 
of  similar  generality  should  now  be  considered  i.igh  priority 
The  Recommendations  section  of  this  report  addresses  this 
need . 

The  following  recommendations  are  made: 

1.  The  Science  Applications,  Inc.  (SAI)  LNG  vapor  disper- 
sion model  should  be  more  thoroughly  evaluated.  This  will 
require  the  cooperation  of  SAI  due  to  the  proprietary 
nature  of  their  computer  programs  which  are  required  for 
solution  of  the  model  equations.  Further  evaluation  of 
the  SAI  model,  or  other  similar  models  based  on  simultaneous 
solution  of  the  mass,  momentum  and  energy  balance  equations 
which  may  become  available,  should  address  the  following 
requirements : 


a.  The  methodology  of  describing  turbulent  mass, 
momentum  and  energy  transfer  should  be  critically 
evaluated.  A literature  search  should  be  conducted 
to  identify  and  evaluate  experimental  data  supporting 
the  assumption  of  first-order  (eddy  diffusivity, 
thermal  conductivity  and  viscosity)  turbulent  transfer 
phenomenological  relationships  for  describing  turbulent 
transfer  in  the  lower  atmosphere. 

b.  An  error  analysis  should  be  done  to  provide  some 
means  for  estimating  the  confidence  level  in  the 
technique  used  to  assign  numerical  values  to  the 
turbulent  transfer  coefficients. 

c.  Sufficient  calculations  should  be  made  with  the 
model  to  determine  the  sensitivity  of  the  results 
predicted  by  the  model  to  uncertainties  in  the  transfer 
coefficients  identified  in  b.  above. 

d.  An  analysis  should  also  be  made  of  the  liquid 
spread,  vapor  generation,  and  heat  transfer  models  used 
in  the  specification  of  the  boundary  conditions  to 
determine  the  sensitivity  of  the  model  predictions. 

e.  The  numerical  stability  and  accuracy  of  the 
algorithm  used  for  computer  solution  of  the  equations 
should  be  criticallv  evaluated. 


2.  A series  of  computations  should  be  made,  using  the 
SAI  model,  of  the  downwind  distance  to  the  time  average 
5%  concentration  level  for  "instantaneous"  LNG  spills  as 
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a function  of  spill  size.  The  range  of  spill  sizes  should 
be  from  10  cubic  meters  to  25,000  cubic  meters  with 
sufficient  points  between  to  adequately  characterize  the 
predicted  relationship  between  flammable  cloud  travel  and 
spill  size. 

3.  The  result  of  2 above  should  be  compared  with  a 
similarly  prepared  relationship  between  flammable  cloud 
travel  and  spill  size  predicted  using  the  Germeles  and 
Drake  model  described  in  this  report.  It  is  anticipated 
that  the  results  will  be  in  substantial  agreement  for  very 
small  spill  sizes  but  the  comparison  should  indicate  the 
smallest  spill  sizes  for  which  significant  differences 
appear  in  predicted  downwind  distance.  Such  a comparison 
should  also  provide  guidance  for  determining  a lower  bound 
on  the  size  of  experimental  spills  which  may  be  required 
to  assess  large  spill  behavior. 

4.  In  anticipation  of  experimental  spills  which  may  be 
required  to  provide  confidence  in  predictions  of  large  spill 
behavior,  an  evaluation  should  be  made  of  the  experimental 
data  requirements  associated  with  verification  of  model 
predictions . 

5.  Additional  experimental  spills  should  be  performed 
only  after  completion  of  the  program  outlined  above,  and 
such  spills  should  be  performed  for  the  purpose  of  model 
evaluation.  Large  "demonstration  spills"  have  been 
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suggested  recently,  largely  as  a result  of  the  variation 
in  predictions  which  has  been  the  subject  of  this  report. 
It  is  the  opinion  of  this  author  that  validation  of  models 
should  still  be  the  primary  goal  of  further  test  programs; 
"demonstration"  of  the  effects  of  large  spills  without 
heavy  reliance  on  models  should  be  avoided. 
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IV.  INTRODUCTION 


A number  of  predictions  of  LNG  vapor  cloud  formation  and 
dispersion  which  might  result  from  a catastrophic  LNG  spill 
onto  water  have  been  published.  Order  of  magnitude  differences 
in  these  predictions  of  the  area  adjacent  to  the  spill  iich  could 
be  exposed  to  flammable  LNG  vapor/air  mixtures  are  significant 
in  the  overall  assessment  of  the  potential  risk  of  marine  trans- 
portation of  LNG. 

With  respect  to  LNG  spills  onto  water,  the  predictions  of 
the  following  groups  have  been  repeatedly  cited  in  the  literature 
related  to  safety  of  marine  LNG  transportation: 

1.  U.  S.  Bureau  of  Mines  - Burgess  et  al . (1,  2) 

2.  American  Petroleum  Institute  - Feldbauer  et  al.  (3) 

3.  Cabot  Corporation  - Germeles  and  Drake  (4) 

4.  U.  S.  Coast  Guard  CHRIS  (Chemical  Hazard  Response 
Information  System  ) - Arthur  D.  Little,  Inc.  (5) 

5.  Professor  James  Fay,  Massachusetts  Institute  of 
Technology  (6) 

6.  Federal  Power  Commission  (7) 

7.  Science  Applications,  Inc.  (8) 

Other  groups  have  published  information  related  to  vapor 
dispersion  from  LNG  spills  onto  \fater  (9,  10,  11).  However, 
these  studies  have  not  resulted  in  predictions  of  downwind 
travel  of  flammable  gas  mixtures  to  be  expected  in  large  acci- 
dent scenarios  and  were  therefore  not  reviewed  in  this  report. 
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I 

In  addition,  numerous  studies  have  been  made  concerning  the 
dispersion  of  vapor  clouds  from  LNG  spills  on  land.  For  an 
extensive  citation  of  such  work,  the  reader  is  referred  to 
U.  S.  Coast  Guard  document  CG-478,  "Liquefied  Natural  Gas  - 
Views  and  Practices  - Policy  and  Safety",  1 February  1976, 
available  from  the  Cargo  and  Hazardous  Materials  Division 
(G-MHM/83) , U.  S.  Coast  Guard,  Washington,  D.C.,  20590. 

Table  IV-1  shows  the  maximum  downwind  distance  to 
the  time-average  5%  vapor  concentration  level  following  an 
instantaneous  25,000  cubic  meter  spill  as  predicted  by  the 
models  suggested  by  the  seven  groups  above.  The  distances, 
with  the  exception  of  the  estimate  attributed  to  SAI,  were 
computed  by  the  author  using  the  procedure  suggested  by  the 
investigating  groups  cited.  The  corresponding  distance  for 
Sai's  model  could  not  be  computed  due  to  the  proprietary 
nature  of  the  SAI  computer  model.  Table  IV-1  therefore 
includes,  for  comparison,  the  distance  predicted  by  SAI  for 
a 37,000  cubic  meter  spill  as  described  in  their  risk  assess- 
ment study  prepared  for  the  Western  LNG  Terminal  Company  (8) . 

In  reviewing  Table  IV-1,  it  should  be  noted  that  the 
meteorological  conditions  suggested  as  applicable  by  the 
groups  are  not  necessarily  the  worst  that  might  have  been 
assumed.  Specifically, 

1.  The  0.75  mile  distance  obtained  with  the  FPC  model 
reflects  the  assumption  of  neutral  atmospheric  stability 
values  (Pasquill  D) , as  recommended  by  the  FPC  staff. 
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TABLE  IV-1 


MAXIMUM  DOWNWIND  DISTANCE  TO  TIME-AVERAGE  5%  CONCPNTRATION 
LEVEL  FOLLOWING  25,000  M3  INSTANTANEOUS  SPILL  OF  LNG  ONTO 
WATER  (Assumes  5 MPH  Wind  except  as  noted  and  Meteorological 
Conditions  Considered  Applicable  by  Investigating  Groups) 


MODEL 

U.  S.  Bureau  of  Mines  (1,  2) 
American  Petroleum  Institute  (3) 
Cabot  Corporation  (4) 

U.  S.  Coast  Guard  CHRIS  (5) 
Professor  James  Fay  (6) 

Federal  Power  Commission  (7) 
Science  Applications,  Inc.  (8) 


DISTANCE  (MILESj 
25.2  - 50.3* 
5.2 
11.5 
16.3** 

17.4** 

0.75 

1.2*** 


* A range  was  presented  to  indicate  uncertainty  in  vapor 
evolution  rate 

**  Wind  velocity  not  considered  explicitly  in  model 

***  For  37,500  cubic  meter  instantaneous  release,  wind 
velocity  « 6.7  MPH 
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2.  The  5.2  mile  distance  obtained  with  the  American 
Petroleum  Institute  Model  assumes  vertical  dispersion 
characterized  as  Singer  and  Smith  D and  horizontal 
dispersion  characterized  as  Pasquill  C,  following  the 
procedure  suggested  by  Feldbauer  (3) . 

3.  The  SAI  prediction  assumes  neutral  atmospheric 
stability;  however  the  result  is  claimed  to  be  essentially 
the  same  when  stable  atmospheric  conditions  are  assumed. 
The  SAI  model  also  gives  longer  downwind  travel  distances 
when  higher  wind  velocities  are  assumed  in  contrast  to 
the  other  models. 

4.  The  predictions  obtained  using  the  Bureau  of  Mines 
(Burgess)  model,  the  CHRIS  model.  Fay's  model,  and  the 
Cabot  (Germeles  and  Drake)  model  assume  stable  atmos- 
pheric conditions. 

The  results  shown  in  Table  IV-1  are  specifically  for  a 
25,000  cubic  meter  instantaneous  spill  (except  as  noted  for 
Sciences  Applications,  Inc.).  These  predictions  do  not 
consider  the  possibility  that  flammable  concentrations  of 
vapor  might  exist  at  greater  distances,  since  the  5%  level 
used  for  the  calculation  must  be  considered  a time-average 
concentration.  Nevertheless,  the  variation  shown  reasonably 
characterizes  the  extreme  range  of  predicted  results  which 
is  the  basis  for  the  present  controversy  regarding  the  assess- 
ment of  the  vapor  cloud  hazard  from  LNG  spills. 
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The  variation  in  these  predictions  is  significant  in 
assessing  the  potential  hazard  associated  with  a large 
accidental  release  of  LNG  on  water  if  the  release  should 
occur  without  immediate  ignition  of  the  flammable  vapor 
mixture  at  the  accident  site.  The  sudden  release  of  large 
amounts  of  LNG  onto  water  is  practically  realizable  only  as 
a result  of  a high  energy  collision.  Immediate  ignition  is 
considered  extremely  likely  if  such  a collision  should  occur/ 
because  of  the  frictional  heating  anticipated  with  such  a 
collision  and  ignition  sources  which  would  result  from 
damaged  equipment.  If  ignition  of  the  vapors  does  not  occur 
at  the  spill  site/  formation  of  a large  vapor  cloud  also 
presupposes  the  virtual  absence  of  ignition  sources  in  thd 
area  close  to  where  the  cloud  is  being  formed.  For  these 
reasons,  an  accident  scenario  which  assumes  formation  of  a 
vapor  cloud  extending  over  large  populated  areas  before 
ianition  is  extremely  unlikely,  even  if  formation  of  such 
clouds  might  occur  in  the  absence  of  ignition. 

However,  the  predictability  of  dispersion  of  vapors  from 
accidental,  catastrophic  release  of  LNG  does  have  a bearing  on 
the  safety  related  management  of  LNG  vessel  traffic,  as  it  does 
on  the  management  of  other  hazardous  cargoes.  This  is  true 
because  the  zone  around  an  accident  which  might  be  subjected 
to  flammable  vapor  concentrations  resulting  fron  non-ignited 
spills  (however  remote  the  probability)  has  an  effect  on  the 
routing  and  traffic  control  of  vessels  and  would  influence 
emergency  response  procedures.  Furthermore,  although  it  appears 
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extremely  unlikely  that  large  populated  areas  could  be  exposed 
to  a flammable  vapor  cloud,  since  ignition  is  likely  when  the 
advancing  front  of  the  cloud  reaches  such  areas,  the  ability 
to  predict  dispersion  is  required  to  assess  the  damage  which 
would  result  from  an  early  ignition.  This  is  true  for  two 
reasons.  First,  the  burning  of  a flammable  cloud  cannot  be 
adequately  predicted  without  knowledge  of  the  composition 
of  the  cloud.  Second,  an  ability  to  predict  vapor  dispersion 
is  required  to  assess  potential  damage  which  might  result  from 
detonations  of  vapor-air  mixtures.  Although  the  experience  to 
date  indicates  that  detonation  of  unconfined  I.NG  vapor-air 
mixtures  is  not  likely,  a good  method  of  vapor  dispersion  pre- 
diction would  be  valuable  in  attempts  to  understand  the 
circumstances  under  which  LNG  vapor-air  (or  a variety  of  other 
fuels  and  chemicals)  detonations  might  be  expected,  such  as 
partial  confinement  and  high  energy  initiation. 

The  purpose  of  this  study  is  fivefold: 

1 . To  provide  a detailed  description  of  the  mathematical 
models  upon  which  published  predictions  of  LNG  vapor  travel 
downwind  of  catastrophic  LNG  spills  onto  water  have  been 
based . 

2.  To  estimate,  using  these  models,  the  maximum  downwind 
travel  of  flammable  LNG  vapojr-air  mixtures  for  a "standard" 
spill  scenario,  so  that  a valid  comparison  can  be  made  of 
the  results  obtained  when  the  different  models  are  used  to 
describe  the  same  event. 
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3.  To  identify  the  reason  for  differences  in  predictions 
which  occur  when  the  different  models  are  used  to  describe 
the  same  event,  and  to  assess  the  technical  credibility 

of  the  methodology  which  results  in  such  differences. 

4.  To  define  the  present  "state  of  the  art"  in  LNG  vapor 
dispersion  modeling,  with  particular  emphasis  on  the  extent 
to  which  the  present  state  of  the  art  justifies  reliance 

on  existing  published  predictions  in  formulating  LNG  safety 
management  programs . 

5.  To  provide  recommendations  for  further  work  which  would 
increase  confidence  in  the  predictability  of  vapor  dis- 
persion from  LNG  (and  other  volatile  chemicals)  accidentally 
spilled  on  water. 

The  "standard  scenario"  LNG  spill  which  is  assumed  in  this 
report  for  purposes  of  comparison  of  the  above  models  is  an 
instantaneous  release  of  25,000  cubic  meters  of  LNG  onto  water. 
25,000  cubic  meters  is  representative  of  the  largest  single- tank 
capacity  of  ships  constructed  to  date  or  on  order.  Although 
as  many  as  six  tanks  may  be  incorporated  into  an  LNG  ship, 
an  accident  resulting  in  simultaneous  rupture  of  more  than  two 
tanks  is  not  considered  credible.  In  the  event  of  simultaneous 
rupture  of  two  tanks,  instantaneous  release  from  both  tanks  is  not 
considered  credible.  The  vapor  travel  following  instantaneous 
release  of  25,000  cubic  meters  would  be  expected  to  be  even  more 
extensive  than  would  be  expected  from  the  actual  release  of  LNG 
following  simultaneous  rupture  of  two  25,000  cubic  meter  tanks. 
Thus,  the  instantaneous  release  of  25,000  cubic  meters  of  LNG 
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on  water  provides  a conservative  upper  limit  on  the  size  of  a 
spill  which  might  conceivably  occur,  even  though  such  a spill 
is  considered  extremely  unlikely. 


V.  BASIS  FOR  LNG  VAPOR  DISPERSION  MODELS 


A number  of  different  predictions  of  LNG  vapor  cloud 
formation  and  dispersion  resulting  from  an  accidental  LNG 
spill  onto  water  have  been  published.  Although  the  pre- 
dictions reflect  wide  disagreement  by  the  parties  involved 
as  to  the  extent  of  hazard  associated  with  downwind  travel 
of  flammable  gas-air  mixtures,  it  is  important  to  realize 
that  all  of  the  mathematical  models  that  have  been  used  to 
make  such  predictions  h-*ve  a common  basis.  It  is  therefore 
expedient  to  provide  the  necessary  physical  and  mathematical 
basis  which  is  common  to  all  models  to  be  discussed. 


V-A.  PHYSICAL  PROCESSES  INVOLVED  IN  LNG  VAPOR  DISPERSION 

LNG  vapor  dispersion  in  the  atmosphere  involves  the 
mixing  with  air  of  a gas  which  is  much  colder  and  denser  than 
air.  A valid  description  of  the  process  should  account  for 
the  following  processes,  which  may  occur  simultaneously: 

1.  Heat  transfer  effects  due  to  mixing  the  cold  gas, 
formed  from  the  boiling  LNG,  with  warmer  air  (which  may  con- 
tain water  vapor) , and  heat  transfer  from  beneath  the  cloud 
(ground  or  water) . 

2.  Gravity-induced  spreading  effects  resulting  from 
non-uniform  density. 

3.  Dispersion  (dilution  with  air)  of  the  gas  due  to 
turbulent  fluid  motion. 


The  heat  transfer  and  gravity  spreading  effects  can  be 
described  in  mathematical  equation  form  by  application  of 
the  principles  of  accountability  of  energy  and  acco  mtability 
of  momentum.  The  concentration  variations  in  the  LNG 
vapor/air  mixture  can  be  described  in  mathematical  equation 
form  by  application  of  the  principle  of  accountability  of 
mass.  The  general  equations  of  accountability  of  energy, 
momentum,  and  mass  are,  respectively} 

Accountability  of  Energy 

-+  _v  i 

3pH  = -VpHv-V*q  + DP  - T : Vv  (V-l) 

3 1 Dt 

1 1 i 1 i l 


I 


II 


III 


Accountability 

of 

Momentum* 

3pv 

3t 

= - v • p vv 

l - 

-V 

-+■ 

-+ 

•t  + pg 

j i i 

(V-2) 

1 i 

II 

Accountability 

of 

Mass** 

3p 

3 1 

-V 

= -v*pv 

(total  mass) 

(V-3a) 

3 C 
at 

• — _l 

= -v*c$ 

i 1 

(methane  or 

LNG  component) 

(V-3b) 

I 

II 

• 

where 

P 

* density  of  air-gas  mixture 

K 

= enthalpy  (energy  content) 
air-gas  mixture 

of 

v = velocity  vector,  decomposable  into 

x,  y,  z components  u,  v,  w,  respectively 

**► 

q =*  heat  transfer  vector,  decomposible 
into  components  qx,  qy,  q2 
P = pressure 

t « time 

* Coriolis  Forces  have  been  neglected 

**  Molecular  Diffusion  has  been  neglected 
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g = gravity  force  vector,  decomposible 
into  3 components  gx“0,  gy=0.  9z 

C = concentration  of  gas-air  mixture 

t = stress  tensor,  decomposible  into  9 
components  txx,  tXy*  ^ i *^yx*  "^yy ^ 

Tyz'  Tzx'  Tzy'  Tzz 

Each  of  the  above  equations  can  be  understood  as  being  a state- 
ment of  the  general  principle  of  accountability: 

"The  rate  of  accumulation  of  a quantity  (energy, 
momentum,  or  mass}  at  a given  location  is  equal  to 
the  net  rate  at  which  the  quantity  is  transferred 
into  (or  out  of)  that  location  from  its  surroundings, 
plus  the  rate  at  which  that  quantity  is  being  produced 
at  that  location."  ’ - - , 

* * 

The  groups  of  terms  labelled  I,  II,  III  respectively  in  ’ * * 
Equations  V-l,  2,  3 correspond  to  the  accumulation,  transfer, 
and  production  mentioned  above  and  are  further  explained  in 
Table  V-l. 

Equations  V-l,  2,  3 are  differential  equations.  Use  of 
these  equations  to  describe  LNG  vapor  dispersion  requires 
specification  of  initial  conditions  and  boundary  conditions. 
Initial  conditions  include  a description  of  the  water 
and  atmospheric  conditions  at  the  site  and  a description 
of  the  initiation  of  the  spill.  Boundary  conditions 
may  also  include  water  and  atmospheric  conditions  but 
in  addition  include  description  of  momentum,  energy , and  mass 
transfers  at  the  boundary  of  the  region  being  modeled.  In 
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IDENTIFICATION  OF  TERMS  IN  ENERGY,  MOMENTUM,  AND  MASS  BALANCE  EQUATIONS  (V- 
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theory,  the  differential  equations,  along  with  presc  ibed 
initial  and  boundary  conditions  can  be  solved  to  obtain  a 
complete  description  of  the  cloud  behavior.  However,  a con- 
siderable quantity  of  additional  information  must  be  available 
as  input  to  the  solution  of  these  equations,  and  solution 
for  the  general  three-dimensional  case  is  not  practically 
realizable . 

V-B.  SIMPLIFICATION  OF  GENERAL  MODEL  FOR  PRACTICAL 

APPLICATION 

In  practice,  simplifying  assumptions  are  made  in  the 
general  mathematical  models  in  order  to  arrive  at.  a model  for 
which  a relatively  few  input  data  are  required  and  which  do  not 
require  excessive  solution  time  ai.d  expense.  It  is  in  the 
simplification  of  the  general  model  for  the  purpose  of  predicting 
LNG  vapor  dispersion  that  the  differences  between  various  in- 
vestigators' predictions  result. 

Classical  prediction  of  pollutant  dispersion  in  the  atmos- 
phere, which  in  theory  should  be  described  with  Equations  V-1,2, 

3,  has  focused  primarily  on  dispersion  of  relatively  small 
quantities  of  material  such  as  smoke,  radioactive  isotopes, 
chemical  fumes  and  dusts.  In  such  situations  it  is  commonly 
assumed  that  rhe  pollutant  material  is  present  in  sufficiently 
small  quantities  that  it  does  not  directly  affect  the  motion  of 
the  air  into  which  it  is  placed.  Rather,  the  notion  is  that  the 
pollutant  particles  simply  follow  the  (already  established)  motion 
of  the  air.  In  such  cases,  material  being  dispersed  can  be 
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viewed  as  "tracers"  of  the  air  motion.  Implicit  in  this 
approach  is  the  assumption  that  there  are  no  heat  transfer 

effects  between  the  pollutant  and  the  air  and  that  there  is 
no  momentum  exchange  between  the  two.  in  such  a case,  the 
requirement  for  the  energy  and  momentum  balances  is  eliminated. 
With  the  additional  assumption  that  the  material  has  the  same 
density  as  air  (i.e.  the  mixture  is  "neutrally  buoyant"), 
the  system  of  Equations  V-l,  2,  3 reduces  to  one  equation  for 
the  conservation  of  LNG  vapor: 

!£-  -V-C v (V-4 ) 


The  quantity  , where  C is  the  concentration  of  the  pollutant 
in  mass  per  unit  volume  and  v is  the  local  velocity,  has  the 
physical  units  of  a mass  flux  term,  i.e.  mass/  area  - time. 

The  quantity  Cv  must  be  related  to  the  concentration  profile  in 
order  to  get  an  equation  which  can  be  solved.  To  this  end 
it  is  commonly  assumed  that  the  mass  flux  is  proportional  to 
the  concentration  gradient.  In  this  case, 

Cv  = -kvC  (V-5) 

where  vc  - the  local  concentration  gradient 

k **  the  "Ficks  Law"  diffusion  coefficient 
Substituting  Equation  V-5  into  Equation  V - 4 we  have 

-J.£  - V • k C 

at  (v-6) 
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V-B-I.  Instantaneous  Release  (Puff)  Model 


If  the  pollutant  is  assumed  to  enter  the  atmosphere 
instantaneously,  in  amount  Q,  from  a single  point,  the  solution 
to  Eg  (V  - 6}  is 


Q 

C (x,y , z ,t)  = exp 

1 r*2  ♦ y2  + Z2' 

(4irt)  3/2  (KxKyK2)  1/2 

«t  [Kx  Ky  Kz 

where  Kv,  K , K_  are  constant  diffusion  coefficients 
a y z 

for  diffusion  in  the  x,  y,  z directions,  respectively 


Furthermore,  arguments  based  on  statistical  analyses  of  random 
turbulence  (13)  indicate  that  for  large  diffusion  times  (how 
large  is  large  enough  depends  on  the  particular  application, 
and  in  any  case  is  not  readily  determined) , the  mean  square 
diffusion  distance  is  given  by 


x2  (t) 


°x  = 2Kxt 


(V-8) 


y2(t) 


0y2  = 2IC  t 

J y 


Z2(t) 


0z2 


= 2K  t 
z 


where  oXf  oy,  az 


are  the  standard  deviations 
of  the  concentration  distributions 
in  the  yz,  xz,  and  xy  planes,  respectively 


36 


Combination  of  Equations  V - 7 and  V - 8 gives 


C(x,y,z,t)  = 


0(2*) 


r3/2 


ax0yaz 


exp 


[•-Hi 


2 + xl  + *L 

2 a y2  cfZ2 


(V~9 ) 


If  the  diffusion  process  is  assumed  to  be  superimposed  oh 
(translated  with)  a mean  wind  in  the  x direction  having 
velocity  u , a coordinate  transformation  gives 


C (x,y ,z ,t) 


Q(2irf  3/2 

exp 

°xlayl0zl 


’(x-ut)  2 


Jxl2 


Jyl2 


(V-10) 


Equation  V-10  predicts  the  concentration  of  the  gas  at  a 
position  x,  y,  z (relative  to  the  release  point)  at  time  t, 
given  an  instantaneous  point  source  release  of  the  gas  of 
magnitude  Q,  mean  wind  speed  u , and  diffusion  coefficients 
CTxI  ' ayi  '°zl  where  the  subscript  I has  been  used  to  denote 
association  with  an  instantaneous  release.  Equation  V-10 
excludes  all  heat  transfer  effects,  momentum  transfer  effects, 
and  gravity  effects  associated  with  materials  having  density 
different  from  air,  and  can  be  viewed  as  describing  the  growth 
of  a puff  or  cloud  as  it  is  carried  downwind  with  the  mean  wind 
velocity  u , as  shown  in  Figure  V-l.  Figure  V-1A  depicts  the 
position  of  the  outer  limit  of  the  cloud  which  continues  to 
increase  in  size.  Figure  V-lB  depicts  the  position  of  the  5% 
isopleth  (line  of  constant  concentration)  as  the  cloud  moves 
downwind.  The  portion  of  the  cloud  at  ground  level  with 
concentration  above  5%  increases  in  size  at  first  due  to  spreading 
of  the  vapors  and  then  shrinks  in  size  due  to  further  dilution 
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size  increases 


with  air.  At  some  distance  downwind,  labeled  in  Figure  V-1B, 

the  entire  cloud  is  below  the  5%  flammable  limit.  The 
distance  is  presumed  to  be  the  downwind  limit  of  the 

hazardous  zone  resulting  from  the  spill. 

Since  the  use  of  the  method  is  dependent  on  the  availability 
of  the  dispersion  coefficients  ( o^,  °yi»  czl^  which  are 
determined  in  practice  from  the  average  behavior  of  a number  of 
puffs,  consideration  must  be  given  to  the  probability  of  single 
puffs  having  downwind  distances  (to  the  5%  limit)  greater 

than  XMAX* 

The  maximum  concentration  predit  ted  by  Equation  V - 10 
occurs  at  the  puff  center  (x=ut , y=o  , z=o) . The  maximum 
concentration  at  a given  downwind  distance  is  therefore  given 
by  the  equation 

3 -V2 

c (x,y=0 , z-0,  t-Vu)  . ° (2"  (V-U) 

°xl°ylazl 

In  Equation  V-ll  the  right  hand  side  has  been  multiplied  by 
2 to  account  for  the  presence  of  the  water  surface. 

Equation  V-ll  is  the  basic  equation  used  by  Germeles 
and  Drake  (4) , the  U.S.  Coast  Guard  CHRIS  model  (5)  and  Fay 
and  Lewis  (6)  for  the  prediction  of  atmospheric  dispersion  of 
LNG  vapor  from  large  rapid  spills.  Fay  does  modify  the 
equation  to  force  asymptotic  behavior  which  he  considers  more 
applicable  to  the  cloud  development.  The  differences  in 
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process  is  rapid  due  to  high  heat  transfer  rates,  the  gas 
release  rate  from  the  spreading  LNG  liquid  pool  is  highly 
transient,  with  the  general  characteristics  shown  in  Fig.  V-2. 
Point  A refers  to  the  instant  when  a quantity  of  MG  is  spilled 
on  the  water.  The  segment  AB  corresponds  to  the  period  of  time 
when  the  liquid  pool  is  spreading,  the  increase  in  release  rate 
during  this  period  is  primarily  due  to  the  increase  in  LK  water 
contact  area  since  evaporation  rate  per  unit  area  of  surface 
(for  large  spills)  is  thought  to  be  relatively  constant.  At 
point  B the  liquid  pool  stops  growing  and  the  gas  release  rate 
remains  constant  until  the  pool  begins  to  break  up  at  point  C. 
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RATE  OF  GAS  RELEASE,  MASS/TIME 


The  decrease  in  release  rate  along  CD  corresponds  to  the  boil 
off  from  the  broken  patches  of  LNG.  At  point  D the  release  is 
complete. 

Several  LNG  evaporation  models  (1,  2,  3,  14,  15,  16,  17, 

18,  19)  have  been  proposed  for  quantitative  prediction  of  the 
pattern  shown  in  Figure  V-2.  Some  of  these  models  have  been 
used  in  predicting  LNG  vapor  dispersion  on  water  and  will  be 
discussed  later  in  this  report.  However,  the  important  point 
to  be  made  here  is  that  the  gas  release  cannot  be  instantaneous 
(the  use  of  the  terminology  "instantaneous  spill",  which  implies 
instantaneous  release  of  a quantity  of  LNG  onto  water  has  often 
been  confused  with  the  terminology  of  "instantaneous  release" 
of  gas  as  implied  by  Equation  V - 10) . A second  feature  of 
importance  illustrated  by  Figure  V-2  is  the  highly  transient 
nature  of  the  gas  release  rate.  This  is  important  because  an 
alternative  approach  (to  the  instantaneous  release  model)  is  to 
model  the  gas  release  rate  as  being  constant  in  time. 

V-B-2.  Steady  Release  (Plume)  Model 

Figure  V-3  illustrates  the  fact  that  a continuous 
release  of  material  can  be  viewed  as  the  rapid  successive 
release  of  (instantaneous)  puffs.  The  concentration  at  a 
given  point  downwind  resulting  from  a rapid  succession  of 
puffs  is  obtained  by  adding  the  contributions  from  all  the 
puffs  to  the  point  in  question.  This  corresponds  to  integration 
of  Equation  V-10  from  time  zero  to  time  infinity.  This 
integration  is  not  straight-forward  since  the  dispersion 
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TIME  DEVELOPMENT  OF  SINGLE  PUFF 


ILLUSTRATION  OF  SUCCESSIVE  PUFFS  AT  AN  INSTANT 


RAPID  SUCCESSION  OF  PUFFS  FORMS  CONTINUOUS  PLUME 


Figure  V-3.  ILLUSTRATION  OF  INSTANTANEOUS  RELEASE  OR  PUFF, 
RELEASE  ~>F  SUCCESSIVE  PUFFS,  AND  A CONTINUOUS 
RELEA?  STEADY  PLUME) 
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coefficients  cxI,  are  functions  of  time  and  distance. 

To  simplify  the  model  development/  it  is  therefore  commonly 
assumed  tha . dispersion  of  each  puff  in  the  downwind  direction 
is  negligible  in  comparison  with  its  movement  associated  with 
the  mean  /ind  velocity  u . The  result  of  the  time  integration 
of  the  instantaneous  release  equation  (Equation  V - 10)  is  the 


equation  widely  used  for  predicting  the  concentration  of  a gas 
or  particulate  material  dispersed  from  a ground-level  point 

source,  in  a wind  with  mean  velocity  u in  the  X direction,  at  a 
constant  rate  Q' : 


C (x,y ,z)  = 


TTO'  a u 

y z 


exp 


5 • i W 


(V-12) 


In  Equation  V-12  the  right  hand  side  has  been  multiplied  by 
2 to  account  for  the  presence  of  the  water  surface. 

Since  this  method  is  also  dependent  on  dispersion 
coefficients  ( oy  , a2  ) which  are  determined  in  practice  from 
the  time  average  behavior  of  a plume,  consideration  must  be 
given  to  the  probability  of  existence  of  gas  pockets  having 
downwind  distances  (to  the  5%  average  value)  greater  than 

The  maximum  concentration  predicted  by  Equation  V-12 
occurs  at  the  plume  center  line  at  ground  level  ( y = 0,  z = 0) . 
The  maximum  concentration  at  a given  downwind  distance  is  then 
given  by  the  equation 

Q' 

C(x,y=0,z=0)  = (V-13) 

TTOyOz  U 
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Equation  V - 13  is  the  basic  equation  used  by  Burgess  et  al. 

(1,  2),  Feldbauer  et  a1.  (3),  and  the  Federal  Power  Commission 
Staff  (7)  to  model  the  dispersion  of  LNG  vapors  from  spills  onto 
water.  The  differences  in  results  obtained  by  these  investigators 
result  from  differences  in  treatment  of  the  rate  of  gas  addition 
to  the  atmosphere,  in  allowance  for  the  effect  of  the  "area" 
source,  in  the  use  of  different  values  for  the  dispersion 
coefficients,  and  in  modifications  intended  to  describe  heat 
transfer  and  non-uniform  density  effects.  These  differences  will 
also  be  considered  in  detail  later  in  this  report. 

Table  V-2  shows  a summary  of  the  input  parameters  which 
must  be  specified  by  the  user  to  predict  downwind  concentrations 
of  vapor  using  Equations  V - 11  and  V - 13.  As. will  be  demon- 
strated in  the  detailed  analyses  of  the  various  predictions 
that  have  been  made  using  Equation  V - 11  and  Equation  V - 13, 
all  of  the  variation  in  the  reported  results  can  be  attributed 
to  the  following  factors. 

1.  The  maximum  downwind  distance  to  the  lower  flammable 
limit  concentration  is  strongly  dependent  on  the  amount  of 
material  released,  Q,  in  Equation  V - 11  or  on  the  rate  cf 
addition  of  LNG  vapor  to  the  atmosphere,  Q' , in  Equation  V - 13. 

2.  The  maximum  downwind  distance  to  the  lower  flammable 
limit  concentration  is  strongly  dependent  on  the  numerical 
values  of  the  dispersion  coefficients,  oxJ,  CTyi'°zi  used  in 
Equation  V - 11  and  a ,a_  used  in  Equation  V-13.  These  dis- 

y z 

persion  coefficients  in  turn  are  strongly  dependent  on  the 
atmospheric  conditions  at  the  spill  site.  Further,  the 
specification  of  these  dispersion  coefficients  in  the  scientific 
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TABLE  V - 2 


INPUT  PARAMETERS  REQUIRED  FOR  VAPOR  DISPERSION  PREDICTIONS 


MODE  OF  RELEASE 


INPUT  PARAMETERS  REQUIRED 


Instantaneous , 

Point  Source 

(Puff) , Equation  V-ll 


Steady,  Point  Source 
(Plume) , Equation  V-13 


xl, 


Q = amount  released 

o T o T = Horizontal  and  vertical 
yl , zl 

dispersion  coefficients 
for  instantaneous  re- 
lease 

Q'  = rate  of  release  of  material 
into  atmosphere 

„ _ Horizontal  and  vertical 


dispersion  coefficients 
for  steady  Plume 
u = mean  wind  velocity 
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literature  is  the  result  of  actual  experimental  measurements 
(see  Appendix  I)  from  releases  of  material  with  essentially 
neutral  buoyancy.  Therefore,  effects  such  as  those  associated 
with  the  low  temperatures  and  high  densities  of  LNG  vapor  are 
not  included  in  literature  compilations  of  dispersion  coefficients. 
Various  attempts,  which  are  empirical  in  nature,  to  account  for 
this  by  ''doctoring"  the  coefficients  obtained  from  neutrally 
buoyant  dispersion  measurements  are  responsible  for  much  of  the 
variation  in  predicted  results  based  on  Equation  V - 11  and  V ~ 13. 

3.  Equations  V - 11  and  V - 13  include  no  provisions  for 
the  LNG  vapor  puff  or  plume  to  spread  due  to  gravity  effects  as 
might  be  expected  due  to  the  density  of  the  cold  LNG  vapors. 
Treatment  of  effects  resulting  from  gravity  spreading  of  the 
vapors  resulting  from  latge  spills  has  varied  widely,  with 
correspondingly  varying  results. 

4.  Equations  V - 11  and  V - 13  assume  entry  into  the 
atmosphere  from  a point  source,  while  an  LNG  spill  onto  water 
is  an  area  source  of  LNG  vapor.  Attempts  to  estimate  the 
effect  of  the  area  source,  while  utilizing  Equations  V - 11 
and  V - 13,  are  responsible  for  some  of  the  variation  in 
reported  predictions  of  downwind  vapor  travel. 

V - B - 3 Combined  Mass,  Momentum,  and  Energy  Balance  Models 

A significantly  different  approach  to  the  prediction  of  LNG 
vapor  dispersion  following  an  accidental  spill,  which  involves 
solution  of  the  system  of  Equations  V-1,2,3  with  less  restrictive 
simplifying  assumptions,  has  been  published  by  Science  Applications, 
Incorporated  (8).  This  approach  results  in  estimates  of  maximum 


47 


downwind  distance  to  the  lower  flammable  limit  that  are  an 
order  of  magnitude  shorter  than  some  of  the  earlier  estimates 
which  were  based  on  the  use  of  Equations  V-ll  and  V-13.  A 
later  section  of  this  report  describes  in  detail  the  methodology 
associated  with  SAX's  predictions.  >■ 


VI.  SURVEY  OF  VAPOR  DISPERSION  PREDICTIONS  ASSUMING 
INSTANTANEOUS  RELEASE  OF  VAPOR  FROM  INSTANTANEOUS  SPILL 
OF  25,000  M3  OF  LNG  ONTO  WATER  - CLASSICAL  PUFF  MODELS 

Germs les  and  Drake  (4)  , the  Coast  Guard  (5)  and  Fay 
and  Lewis  (6)  have  published  predictions  for  a "worst  ca.se", 
instantaneous,  release  of  25,000  cubic  meters  of  LNG  onto 
water.  Table  VI  - 1 is  a summary  of  the  vapor  dispersion 
predictions  obtained  using  the  models  suggested  by  Fay  and 
Lewis,  Germeles  and  Drake  and  CHRIS  for  a 25,000  M3 
instantaneous  spill  onto  water  during  stable  atmospheric 
conditions.  Table  VI  - 2 presents  results  predicted  for  neutral 
weather  conditions.  All  three  groups  assume  applicability  of 
Equation  V - 10,  the  puff  model,  to  the  dispersion  of  the  vapor 
following  a spill.  Germeles  and  Drake  and  the  Coast  Guard 
CHRIS  method  assume  the  vaporized  LNG  initially  forms  a 
cylindrical  pancake  of  radius  re  and  height  he  where  re  is 
equal  to  the  radius  of  the  liquid  pool  at  the  end  of  the 
evaporation  period.  The  equations  used  by  Germeles  and  Drake 
and  the  CHRIS  model  for  predicting  the  radius  of  the  pool  at 
the  end  of  the  evaporation  period  (maximum  pool  radius)  are 
shown  in  Table  VI  - 3.  Table  VI  - 3 also  includes,  for 
[ i comparison,  other  models  for  maximum  pool  radius  which  have 

appeared  in  the  literature.  Table  VI  - 4 shows  the  maximum 
pool  radius,  evaporation  time,  and  height  (assuming  a 
cylindrical  cloud  of  pure  LNG  vapor  at  its  boiling  point) 

1 calculated  for  a 25,000  M3  spill  using  the  equations  shown 

» in  Table  VI  - 1. 
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(2)  Germeles  and  Drake  Model,  using  Gif ford-Pesquill  "D-Neutral"  Plume  Dispersion 
Coefficients  (Appendix  1) 

(3)  CHRIS  Model,  using  Gif ford-Pasguill  "D-Neutral"  Plume  Dispersion  Coefficients  (Appendix  1) 
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PUBLISHED  MODELS  FOR  PREDICTING  EVAPORATION 
TIME  AND  MAXIMUM  POOL  SIZE  FOR  INSTANTANEOUS 
SPILLS  OF  LNG  ONTO  WATER 


EQUATION  FOR 
MAXIMUM  RADIUS 

EQUATION  FOR 
EVAPORATION  TIME 

SOURCE 

(1) 

re 

7.4  V3/8 

v4 

t - 8.8  V 

Raj/Kalelkar  (15) 

h1^ 

hV2 

(used  by  Germeles- 
Drake  and  CHRIS) 

(2) 

re  " 

5/12 

4.7  V 

te  = 3.3  V1/3 

Fay  (14) 

(3) 

re  = 

10.4  v"/l2 

te  = 14.5  VV3 

Hoult  (16) 

(4) 

r = 
e 

7.3  V3//g 

t = 7.9  V1/4 
e 

Hoult  (17) 

(5) 

re  = 

7.6  V3/8 
h1/8 

1 / M 

tc  * 12.4  V 
h1'* 

Otterraan  (18) 

(6) 

r = 
e 

9.07  V3/8 
h1" 

t.  = 10.56  V1/4 
e 

hX/2 

Muscari  (19) 

where  V 

= Volume  of  Spill,  ft3  LNG 

re 

= Maximum  Pool  Radius,  ft 

*-e 

= Evaporation  Time,  sec 

h 

= Liquid  Regression  Rate, 

in/min 

52 


. » *:  mix  - 


TABLE  VI  - 4 


PREDICTION  OF  INITIAL  LNG  VAPOR  CLOUD  SIZE 
FOLLOWING  INSTANTANEOUS  SPILL  OF  25,000  M3  ON 
WATER  USING  MODELS  OF  TABLE  VI-I 


SOURCE 

VAPOR  CLOUD  RADIUS 
(FT) 

EVAPORATION 
TIME  (SEC) 

VAPOR  CLOUD 
HEIGHT  (FT) 

Raj/Kalelkar  (15) 
(used  by  Germeles- 
Drake  and  CHRIS) 

1255 

270 

43 

Fay  (14) 

1417 

316 

34 

Hoult  (16) 

3136 

1390 

7 

Hoult  (17) 

1239 

242 

44 

Otterman  (18) 

1289 

380 

41 

Muscari  (19) 

1539 

324 

29 

* Vapor  Cloud  Haight  = 241  vxiq^Tre^ 

where  241  ■ Gas  Specific  Volume  at  Boiling  Point 

Liquid  Specific  Volume  at  Boiling  Point 


Fay  and  Lewis  and  Germeles  and  Drake  assume  the  vapor 
generated,  being  heavier  than  air,  will  spread  laterally 
across  the  water  surface.  Fay  assumes  the  cloud  spreads 
without  appreciable  mixing  with  air,  while  Germeles  and  Drake 
allow  for  air  entrainment  at  the  top  of  the  cloud.  The 

termination  of  the  so-called  "gravity  spread"  phase  was 
considered  by  Fay  and  Lewis  to  be  the  point  where  the  cloud 
becomes  neutrally  buoyant  due  to  heat  transfer  from  the  water 
below  the  cloud.  Germeles  and  Drake  terminated  this  phase 
of  their  model  at  the  point  where  the  cloud  becomes  buoyant 
under  no  wind  conditions,  or  when  wind  is  present,  at  the  point 
where  the  gravity  spread  velocity  of  the  cloud  equals  the  mean 
wind  velocity.  Fay  models  the  "warming  up"  process  of  the 
cloud  as  resulting  only  from  convective  heat  transfer  between 
the  water  surface  and  the  cloud,  while  Germeles  and  Drake 
consider  heat  effects  due  to  convection  and  mixing  with  entrained 
air,  including  the  latent  heat  effects  of  condensation  and 
freezing  of  water  vapor.  The  vapor  cloud  at  the  end  of  the 
gravity  spread  phase,  as  described  by  Fay  and  Lewis  and 
Germeles  and  Drake  in  Items  II  and  III  in  Table  VI  - 1 , are 
used  as  starting  points  for  their  models  of  the  atmospheric 
dispersion  phase.  All  three  use  vapor  dispersion  models  based 
on  Equation  V - 1G,  restated: 


-3/2 

Q ( 2 TT  ) 
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Fay  notes  that  the  maximum  concentration  at  any  point 
downwind  must  occur  at  ground  level  at  the  cloud  (puff)  center. 


or  x=ut,  y=0,  z~0,  and  he  assumes  = ffxi*  Both  assumptions 
are  widely  practiced  and  appear  to  be  justified  within  the  error 
of  existing  experimental  data  when  other  assumptions  of  the 
model  (neutral  buoyancy,  dilute  pollutant/air  concentrations) 
are  valid. 

Equation  VI  -1  then  becomes 


'tn 


(VI  - 2) 


(2v  3) 


L/2. 


Jyi 
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where  the  subscript  m denotes  "maximum" 


In  Equation  VI  - 2 the  right  hand  side  has  been  multiplied  by 
2 to  account  for  the  presence  of  the  water  surface. 

Since  Q is  the  amount  of  LNG  vapor  added  (instantaneously, 
according  to  the  development  of  Equation  V - 1)  and  the  cloud 
is  assumed  to  be  pure  at  the  end  of  the  gravity  spread  phase. 
Fay  substitutes 


Q = w r_  h„ 
vm  v 


where  r^  and  hv  are  the  radius  and  height  of  the  pure  vapor 
cloud  to  get 


o r 


. Emmii ■ 


Equation  VI  - 3 includes  the  effect  of  horizontal  and 
vertical  mixing.  Fay  argues  that  horizontal  mixing  would  be 
suppressed  near  the  spill  due  to  the  shallow  depth  of  the 
cloud.  He  assumes  that  near  the  spill  only  vertical  mixing 
would  occur  and  that  the  resulting  vertical  distribution  would 
be  Gaussian.  Hence  at  these  ’’intermediate"  distances  Fay 
argues  that  the  maximum  ground  level  concentration  would  be 


Finally,  Fay  argues  that  at  the  location  of  the  spill  the  concentra- 
tion must  be  unity  (i.e.  1 FT3  of  LNG  vapor  per  FT3  of  space).  Based 
on  consideration  of  Equations  VI-3  and  VI-4  and  the  requirement 
of  unity  concentration  at  the  source.  Fay  proposed  the 
following  modified  form  of  Equation  VI  - 1,  which  asymptotically 
yields  Cm  = 1 at  the  spill  location.  Equation  VI  - 4 at 

intermediate  distances (where  a <<r  and  o >>h  ) and  Equation 

yl  vm  zl  v 

VI  - 3 at  large  distances  from  the  spill  source  (where  ay.I>>rvin  and 


Since  Fay  wanted  to  compare  his  model  prediction  with  the 
experimental  data  reported  by  Feldbauex  (3) , he  assumes  that 
at  large  distances  from  the  spill  the  time  average  concentration. 
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C,  of  a passing  cloud  is  given  by 


C = cm 
2 

The  equation  for  C corresponding  to  Equation  VI  - 5 then 
becomes 


r 

2 

r h 

c 

= 

vm 

V 

rvm  + 2 ayj 

(VI  - 6) 


Fay  compares  results  obtained  from  his  analysis  of  Feldbauer's 
data  with  the  prediction  of  Equation  (VI  - 6)  and  obtains  rough 
agreement. 

In  order  to  define  the  maximum  downwind  flammable  extent 
of  the  cloud,  we  are  primarily  interested  in  the  prediction  o£.~' 
the  distance  at  which  C - 0.0  5 (or  some  fraction  theregf*,* 

ill  m- 

K 

depending  on  assumptions  of  peak-to-average  concentration  ratio) 

/ 

as  predicted  by  Equation  VI  - 5.  In  order  to  solve  for  this 

distance,  the  dispersion  coefficients  o and  o T must  be 

yi 

specified  as  functions  of  the  downwind  distance.  Fay  assumes 
Oyj  and  azl  to  be  given  by  the  following  equations  (see 
Appendix  I) . 


ayl 

azl 


Neutral  Stability 


0.06 


0.92 

x 


0.15 


0.70 

x 


Very  Stable 


0.02 

0.05 


x0-89 

0.61 

x 


(VI  - 


7) 


where  a , o , x are  in  meters 
yi  zl 
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Equations  VI  - 7 are  estimated  correlations  for  instantaneous  - 


source  values  of  o T and  <r  T given  by  Slade  (13)  and  are  based 

y I z 

largely  on  the  data  of  Hogstrom  (20) . The  correlations  for 
the  "very  stable"  condition  denote  the  approximate  limit  of  the 


most  stable  data  of  Hogstrom. 

The  solution  of  Equation  VI  - 5 for  Cm  as  a function  of 
x,  the  downwind  distance,  assuming  "very  stable"  coefficients, 
is  shown  in  Figure  VI  - 1 for  a 25,000  instantaneous  LNG 
spill.  From  Figure  VI  - 1,  the  distance  to  a maximum  concentra- 
tion of  5%  is  predicted  to  be  25,000  meters.  The  distance  to 
Cm  = 0.025  (incorporating  a 2 to  1 peak-to-average  ratio)  is 
50,000  meters.  The  distance  to  Cm  = 0.01  (incorporating  a 
5 to  1 peak-to-average  ratio)  is  100,000  meters. 

Values  for  r^  and  hv  of  816  meters  and  2.9  meters 
respectively  were  used  in  the  calculation  of  downwind  distances 
shown  in  Figure  VI  - 1 and  Tables  VI  - 1 and  VI  - 2 . These 
values  were  taken  from  Lewis'  thesis  (40)  . The  values  of  2.9  meters 
for  hv  and  816  meters  for  r^  correspond  to  a pure  vapor  cloud 
volume  at  the  LNG  boiling  point,  approximately  240  times  the  spilled 
liquxd  volume.  Although  Fay  and  Lewis'  paper  indicates  (on  the  last 
line  of  page  491  in  Reference  6)  that  a pancake  neutrally 
buoyant  pure  vapor  cloud  of  radius  r^  and  height  hv  forms  over 
the  spill,  the  results  which  appear  in  Lewis'  thesis  and  which 
correspond  to  Figure  5 of  Fay  and  Lewis'  paper  apparently  are 
based  on  values  of  r^  and  hv  of  816  meters  and  2.9  meters 
respectively.  If  the  height  of  the  cloud  is  determined  from 
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MAXIMUM  DOWNWIND  CONCENTRATION,  VOL.  %/100 


the  volume  at  the  temperature  corresponding  to  neutral 
buoyancy  (v-l55F)  the  height  is  increased  by  a factor  of 
about  1.45  and  the  predicted  downwind  distance  with  very 
stable  weather  increases  to  approximately  23  miles.  In  a 
recent  communication , Dr.  Fay  suggests  that  the  hv  should 
be  determined  from  the  volume  of  the  pure  gas  cloud  at 
0°C.  In  this  case  hv  is  estimated  to  be  7.1  meters  and  the 
downwind  distance  to  the  5%  level  with  stable  weather 
conditions  is  calculated  to  be  28.0  miles.  (See  Appendix  II) 
Fay's  rationale  for  the  development  of  Equation  VI  - 5 
included  uhe  requirement  that  it  agree,  at  long  distances, 
with  Equation  VI  - 3,  restated 
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Recalling  that 


where  V = volume  of  gas  released 

r^  - radius  of  pure  gas  cloud  at  end 
of  gravity  spread  phase 

hv  = heiyht  of  pure  gas  cloud  at  end 
of  gravity  spread  phase, 

Equation  VI  - 8 is  equivalent  to 


(VI 


m 


V 


( 2ir ' 


V; 


2 

°yl  a 


zl 


(VI 


8) 


9) 


Note  that  Equation  VI  - 9 is  just  the  equation  for  the  maximum 
downwind  concentration  (ground  level,  z = 0,  and  cloud  center. 
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60 


y = 0}  for  a puff  or  instantaneous  release  of  vapor  volume  V. 
The  solution  of  Equation  VI  - 9 is  plotted  in  Figure  VI  - 1 
for  comparison  with  Equation  VI  - 5,  using  a , ozI  suggested 
by  Slade  for  "very  stable"  atmospheres.  Although  Fay's  model 
(Equation  VI  - 5)  approaches  the  solution  to  Equation  VI  - 9 
for  very  long  distances,  the  distances  predicted  for  5%  con- 
cer‘  ition  (lower  flammable  limit  for  methane)  for  a 25,000 
cubic  meter  spill  is  significantly  different  for  the  two 
equations  (by  a factor  of  3.7).  The  important  point  to  be 
made  is  that  Fay's  model  can  be  viewed  as  a model  for  the 
int  source  instantaneous  release  of  25,000  LNG  as 


vapor,  modified  to  give  a finite  concentration  (Cm  = 1)  at 
the  source. 


VI-B . PREDICTIONS  USING  GERMELES  AND  DRAKE'S  MODEL 
Although  the  final  prediction  of  downwind  distance  to  a 
given  concentration  by  Germeles  and  Drake  is  also  based  on  the 
use  of  Equation  VI  - 1 , the  classical  diffusion  model  for  the 
dispersion  of  a "puff"  (instantaneous  release  of  vapor) , other 
procedures  in  their  model  differ  significantly  from  those  of 


1.  Germeles  and  Drake  allow  for  entrainment  of  air  by 
the  LNG  cloud  as  it  spreads  across  the  water  surface 
immediately  following  the  release  (which  is  treated  as 
if  the  vapor  release  is  instantaneous) . This  results 
in  a cloud  which  is  to  be  used  as  the  start  of  the 
dispersion  phase  (to  be  described  by  Equation  VI  - 1) 
that  is  already  diluted  with  air. 


2.  Germeles  and  Drake  terminate  the  initial  gravity 

spread  phase  of  the  cloud  (during  which  time  air  is 

/ 

entrained)  at  the  point  where  it  reaches  neutral  density 
under  no-wind  conditions  or  when  the  velocity  of  the 
edge  of  the  spreading  cloud  falls  to  the  mean  wind 
velocity. 

3.  Germeles  and  Drake  argue  that  an  analysis  of  Hogstrom's 
data  for  dispersion  coefficients  for  instantaneous  release 

I 

do  not  justify  Slades  estimated  correlation,  particularly 
for  "very  stable"  weather: 


0.02  x 


0.89 


0.05  x 


0.61 


Instead,  they  recommend  the  use  of  the  Pasquill  F 
stability  "plume"  dispersion  coefficients  for  stable 
weather  and  Pasquill  D stability  coefficients  for  neutral 
weather  conditions. 

4.  Since  Drake  and  Germeles  assume  a "starting  point" 
cloud  of  22%  vapor  (specific  to  the  case  being  considered) 
for  the  atmospheric  dispersion  prediction,  they  correct 
the  result  predicted  by  Equation  VI  -1  by  subtracting 
the  distance  required  for  a 22%  concentration  to  occur 
downwind  of  a point  source  instantaneous  release.  This 
method,  usually  referred  to  as  the  specification  of  a 
"virtual  source",  is  a common  practice  for  allowing  for 
the  effect  of  an  area  source.  The  method  is  illustrated 


in  Figure  VI  - 2 , where  X is  the  distance  computed  for 

C =0.05  using  Equation  VI  - 1 and  X,,  is  the  distance 
m v 

computed  for  = 0.22.  The  downwind  distance  to 
Cm  = 0 . C 5 from  the  actual  (area)  source  is  then 


X’  = X - xv 


(VI  - 10) 


As  shown  in  Table  VI  - 1,  the  size  of  the  initially 
formed  pure  LNG  vapor  cloud  over  the  spill  has  a radius  of 
383  meters  and  a height  of  13  meters.  Germeles  and  Drake 
assume  that  during  the  gravity  induced  spread  the  cloud  can 
be  represented  by  its  average  spatial  thermodynamic  state. 

That  is,  the  cloud  at  any  instant  is  assumed  uniform  in 
temperature  and  composition. 

The  equation  used  to  predict  the  gravity  spread  of  the 
vapor  cloud  was  proposed  by  Yih  (22)  to  describe  density  intrusion 
weather  phenomena  such  as  the  movement  of  cold  fronts : 


dR 

dt 


where 


H = cloud  height,  ft 

t = time,  sec 

3 

p = cloud  density,  lb/ ft 

k - 2 


(VI  - 11) 
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The  cloud  density,  p,  varies  due  to  air  entrainment  and  heat 
transfer  between  the  cloud  and  its  surroundings. 

Germeles  and  Drake  assume  that  air  is  entrained  at  the 
upper  surface  of  the  spreading  LNG  cloud  as  the  clouds  spread 
laterally.  If  the  volume  of  air  entrained,  d Qe,  by  an  annular 
area,  2rrdr,  of  the  top  surface  of  the  cloud  is 


dQe  = « Uc  2rrdr 


(VI  - 12) 


where  = local  velocity  of  the  cloud  surface 

r r dRi 

assumed  = _ — 

L R dt  _J 

= entrainment  coefficient 


then 
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From  the  principle  of  mass  conservation. 


dM 

dt 


= p Q 
Ma  e 


where  M = mass  of  the  "mixed"  cloud 


(VI  - 13) 


(VI  - 14) 


(VI  - 15) 
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From  the  energy  conservation  principle. 


d(MCT) 


dt 


CaVeTS  + 6V  + <K, 


(VI  - 16) 


where 
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m — 
x — 

a = 
0v“ 

«„  - 


heat  capacity  of  mixed  cloud 

temperature  of  mixed  cloud 

refers  to  air  only 

heat  of  condensation  and  freezin 
of  water  in  cloud 

heat  transferred  by  convection, 

natural  (Q  ) or  forced  (£)*)  , 
n *• 


Qn 

Qf 
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Solution  of  the  four  simultaneous  equations; 
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(VI  - 21) 


gives  the  concentration  and  temperature  development  of  the 
cloud  during  the  gravity  spreading  phase. 

Figures  VI  - 3 and  VI  - 4 show  the  development  of  the 
cloud  radius  (R) , temperature  (T) , height  (H)  and  concentration 
(C)  for  a 25,000  M3  instantaneous  spill,  using  Germeles  and 
Drakes'  gravity  spread  model.  The  initial  temperature  and 
humidity  of  the  air  and  the  values  used  for  “ , the  entrainment 
coefficient,  and  f,  the  friction  factor  in  Equation  VI  - 17, 
are  also  shown  in  Figures  VI  - 3 and  VI  - 4. 

Solution  of  Equation  VI  - 1 for  y = z = 0 and  x = ut 

gives 


c = (VI  - 22) 

m 17 

->  1/2  2 
( 2 rr  ) °yl0zl 

Q in  Equation  VI  - 22  was  assumed  by  Germeles  and  Drake  to  be 
the  volume  of  the  pure  vapor  cloud  formed  at  ambient  conditions 
(70°F,  1 atm)  or  approximately  630  times  the  spilled  liquid 
volume.  Solving  Equation  VI  - 22  by  trial  and  error,  using 
_ -the  - Clifford.  JPaLS.qu.i.1  ^cor  r e 1 at  i on  s for  cy,az  vs  X (Appendix  I) 
for  "D-Neutral"  weather  for  C_  = 0.22  gr tffes~X^=  0 0 0 meters. 

Solvina  Equation  VI  - 22  similarly  for  Cm  = 0.05  gives  X 2"4(^)0 
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900  L-  INITIAL  COMblTIONS  FOR  GRAVITY  SPREAD  CALCULATION 


INITIAL  CONDITIONS  FOR  GRAVITY  SPREAD  CALCULATION 
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E FROM  COMPLETE  EVAPORATION,  SECONDS 


meters.  The  downwind  distance  to  5%  concentration  for  neutral 

•> 

stability  conditions  is  then  9800  - 5000  = 4800  meters 
(3.0  miles)  as  shown  in  Table  VI  - 2.  A similar  calculation 
shows  the  distance  to  the  5%  concentration  for  moderately 
stable  (F)  conditions  to  be  about  10,000  meters  (^11.5  miles) 
as  shown  in  Table  VI  - 1. 

Calculations  were  made  to  determine  the  effect  of 
variation  in  a , the  entrainment  coefficient,  on  the  downwind 
distance  to  the  LFL  as  predicted  by  Germeles  and  Drake. 

Figure  VI  - 5 shows  the  average  concentration  of  the  cloud 
during  the  gravity  spread  phase  for  values  of  “ of  0.01, 

0.05,  0.1,  0.20,  and  0.50.  The  first  vertical  hash-mark  on 
each  « curve  on  Figure  VI- 5 denotes  the  time  (and  concentra- 
tion) when  the  gravity  spread  phase  would  be  terminated  for  a 
10  MPH  wind.  The  second  vertical  hash-mark  on  each  « curve 
denotes  the  time  (and  concentration)  when  the  gravity  spread 
phase  would  be  terminated  for  a 5 MPH  wind.  For  “ = 0.5,  the 
downwind  concentration  drops  below  5%  before  the  cloud  edge 
velocity  decreases  to  5 MPH  and  before  the  cloud  becomes 
neutrally  buoyant. 

VI  - C.  PREDICTIONS  USING  U.S.  COAST  GUARD  (CHRIS)  MODEL  (5) 

The  U.S.  Coast  Guard  has  published  methods  for  estimating 
downwind  dispersion  of  vapors  from  spills  of  LNG  or  other 
cryogenic  liquids  in  its  "CHRIS”  - Chemical  Hazards  Reponse 
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RAINMEIWT  COEFFICIENT  ON  GRAVITY  SPREAD 
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System.  These  methods  were  developed  by  Arthur  D.  Little,  Inc. 
under  contract  to  the  Coast  Guard.  ADL's  model  for  vapor 
dispersion  from  an  instantaneous  LNG  spill,  as  incorporated  in 
CHRIS,  is  also  based  on  Equation  VI  - 1,  the  classical 
diffusion  model  for  the  dispersion  of  a "puff"  (instantaneous 
release  of  vapor) . 

To  determine  the  downwind  distance  to  the  5%  (average) 
concentration.  Equation  VI  - 1 was  simplified  for  the  ground 
level,  centerline  case  (x  = ut , y = 0,  z = 0) 
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m 


( 2 it 3 > 2 , 

yl  zi 


(VI  - 23) 


In  Equation  VI  - 2 3 the  right  side  has  been  multiplied  by  2 to 
account  for  the  presence  of  the  water  surface. 

For  a 5 MPH  wind  and  stable  weather  conditions,  values 
of  the  downwind  distance  are  assumed  until  Equation  VI  - 23 
predicts  5%  concentration.  The  dispersion  coefficients  for 
stable  weather  conditions  are  taken  from  the.Pasquill  plume 
dispersion  coefficient  charts  shown  in  Appendix  I.  It  was 
recognized  by  ADL  (5)  that  the  application  of  these  coefficients 
t5''>fcla^  dispersion  of  a puff  (instantaneously  released  vapor) 
is  debatable,  I5"ut  such  use  was  recommended  until  more  experi- 
mental data  are  available.  Using  the  coefficients  representing 
Pasquill  F stability,  the  downwind  distance  to  the  5%  (average) 
concentration  is  determined  (by  trial  and  error)  to  be  approxi- 
•-juately  30  ,000  meters. 
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The  CHRIS  model  accounts  for  the  area  source  nature  of 


r 


an  evaporating  LNG  P°°l  by  locating  a virtual  source  at  a 
distance  5 pool  diameters  upwind  of  the  center  of  the  pool 
as  shown  in  Figure  VI  - 6.  The  liquid  pool  diameter  is 
estimated  using  the  maximum  pool  radius  model  proposed  by 
Raj  and  Kalelkar  (15)  shown  in  Table  VI  ~ 3/  restated: 


r 

e 


7.4 


(VI  - 24) 


where  V = volume  of  spill,  ft^  LNG 

r = Maximum  Pool  Radius , ft 
s 

h = Liquid  Regression  Rate,  in  /min 
For  a 25,000  cubic  meter  instantaneous  spill.  Equation  VI  - 24 
gives  a maximum  pool  diameter  (2re)  of  766  meters.  The 
distance  between  the  pool  center  and  the  virtual  source  is 
then  5 pool  diameters,  or  3830  meters.  Subtracting  this 
distance  from  the  result  given  by  Equation  VI  - 23,  the  downwind 
distance  is  26,200  meters  or  16.3  miles,  as  shown  in  Tabie  VI  - 1 . 


The  CHRIS  model  described  above  assumes  the  instantaneously 
formed  cloud  to  be  at  ambient  temperature  and  pressure  (70°F, 
1 atm) , and  there  is  no  provision  for  gravity  spreading  or 
heat  transfer  effects. 
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Figure  VI-6.  SCHEMA  t 1C  OF  LNG  POOL  AND  LOCATION  OF  VIRTUAL 
POINT  SOURCE  FOR  25,000  m3  INSTANTANEOUS  SPILL 
AS  SUGGESTED  BY  ADL,  INC.  IN  CHRIS  MODEL 


VI -D  COMPARISON  OF  RESULTS  BASED  ON  INSTANTANEOUS 


VAPOR  RELEASE  MODELS 

The  predictions  shown  in  Table  VI  - 1 and  VI  - 2 by 
Germeles  and  Drake,  Fay,  and  CHRIS  of  maximum  downwind  distance 
to  the  5%  and  2 1/2%  time-average  vapor  concentration  following 
instantaneous  release  of  25,000  M3  of  LNG  as  vapor  during  neutral 
and  stable  atmospheric  conditions  appear  to  indicate  fair  agree- 
ment. The  maximum  variation  is  about  25%  from  the  mean  value  for 
the  downwind  distance  to  the  5%  vapor  concentration  during  stable 
weather  conditions.  However,  this  "agreement"  is  due  to  com- 
pensating differences  in  the  approaches. 

The  gravity  spread  portion  of  the  Germeles  and  Drake 
model  determines  the  concentration  which  is  assumed  to  represent 
the  starting  point  for  dispersion  resulting  from  atmospheric 
turbulence.  This  estimated  concentration  (22%  for  the  conditions 
chosen  for  illustration)  directly  determines  the  virtual  source 
correction  as  indicated  in  Figure  VI  - 2.  The  virtual  source 
distance  Oty.  of  Figure  VI  - 2)  for  a 25,000  cubic  meter  instan- 
taneous release  during  stable  weather  co  :tions,  using  the 
Germeles  and  Drake  model,  is  approximately  x4,000  meters  or  8,5 
miles.  The  CHRIS  model,  however,  estimates  the  virtual  source 
distance  to  be  five  pool  diameters  or  approximately  3,800  meters 
(~2.a  miles).  Since  the  estimation  techniques  do  not  difrer 
except  in  estimating  the  virtual  source  distance,  the  difference 
in  predictions  by  CHRIS  arl  Germeles-Drake  can  be  directly 
attributed  to  the  greater  virtual  source  correction  resulting 
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from  the  gravity  spread  effects  included  in  the  Germeles 
Drake  model. 

A "comparison"  of  the  predictions  of  the  Germeles  and 
Drake  and  Fay  models  is  more  difficult.  Four  factors  affect- 
ing the  predictions  of  these  models  must  be  recognized. 

1.  Fay's  modification  of  the  classical  dispersion 
equation  to  force  a unity  concentration  at  the  source 
tends  to  shorten  his  predicted  distances  in  comparison 
to  the  results  obtained  with  simple  application  of  the 
puff  model  (as  shown  in  Figure  VI  - 1)  and  the  model 
of  Germeles  and  Drake. 

2.  Fay's  model  has  been  used  in  this  report  assuming 
the  total  vapor  volume  released  from  the  spill  to  be 
the  saturated  vapor  volume  of  LNG  at  1 atmosphere 
pressure,  or  approximately  240  times  the  liquid  volume. 

The  total  volume  of  vapor  released  from  the  spill  is 
assumed  in  the  Germeles  and  Drake  predictions  to  be  the 
volume  of  methane  at  standard  conditions  (70°F,  1 atm) 
or  approximately  630  times  the  liquid  volume.  If  the 
larger  volume  is  used  in  Fay's  model,  as  suggested  by 
Fay  in  a recent  communication  to  this  author  (see 
Appendix  2,  Fay's  comments),  a much  longer  distance 
(^28  miles)  results. 

3.  Fay  uses  the  "very  stable"  category  puff  dispersion 
coefficients  presented  by  Slade.  Germeles  and  Drake  argue 
in  their  paper  that  the  very  stable  puff  dispersion 
coefficients  correlation  suggested  by  Slade  is  not  suffi- 
ciently justified  from  an  analysis  of  the  original  data. 
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and  that  the  Pasquill  F stability  coefficients  which 
represent  "plume"  dispersion  data  are  more  applicable 
in  their  analysis.  This  choice,  however,  considerably 
shortens  the  downwind  distance  to  the  5%  level  wl  .n 
using  the  Germeles  and  Drake  model.  If  the  very  stable 
puff  dispersion  coefficients  of  Slade  are  used  in 
Germeles  and  Drake's  model  the  calculated  distance  to 
the  5%  level  shown  in  Table  VI  - 1 would  be  approximately 
40  miles.  Conversely,  if  the  Pasquill  F stability 
coefficients  are  used  in  Fay's  model  instead  of  the  very 
stable  puff  coefficients  cited  by  Slade,  the  predicted 
distance  is  cut  roughly  in  half  as  shown  in  Figure  VI  - 7. 

4.  Fay's  model  does  not  address  the  possibility  of  air 
entrainment  during  the  gravity  spread.  This  factor 
considered  alone  would  tend  to  give  a longer  distance 
using  the  Fay  model  than  the  Germeles  and  Drake  model. 

In  view  of  these  important  differences  in  the  three  models 
particularly  the  differences  between  the  Germeles  and  Drake 
and  Fay  models,  the  "agreement"  indicated  in  Table  VI  - 1 and 
Table  VI  - 2 must  be  considered  fortuitous. 
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VII.  SURVEY  OF  VAPOR  DISPERSION  PREDICTIONS 
ASSUMING  STEADY  RELEASE  OF  VAPOR  FROM 
INSTANTANEOUS  SPILL  OF  25,000  M3  LNG 
ONTO  WATER-CLASSICAL  PLUME  MODELS 


Burgess  (1,  2),  Feldbauer  (3),  and  the  Federal  Power 
Commission  (7)  have  published  predictions  for  a "worst  case. 
Instantaneous,  release  of  25,000  cubic  meters  of  LNG  onto 
water.  All  three  assume  the  applicability  of  the  classical 
steady  release  or  "plume"  model  for  atmospheric  dispersion, 
Equation  V-12,  restated: 


C(x,y,z)  = 


trcr  a u 

y z 


exp 


[- ; fy’  ■ i w 


(VII-l) 


Table  VII-l  is  a summary  of  the  vapor  dispersion  predictions 
obtained  using  the  models  suggested  by  Burgess,  Feldbauer,  and 
the  FPC  for  a 25,000  cubic  meter  spill  of  LNG  onto  water. 

The  differences  in  downwind  distances  to  the  5%  concentra- 
tion level  shown  in  Table  VII-l  can  be  attributed  to  four  factors: 

1.  The  value  of  Q',  the  rate  of  vapor  flow  into  the 
atmosphere,  has  been  estimated  by  different  methods, 
with  widely  varying  results.  In  all  cases,  however, 
the  predictions  reflect  the  assumption  of  a steady 
vapor  flow  rate  from  the  spill  site.  (This  condition 
is  implicit  in  Equation  VII-l) . 
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TABLE  VII  - 1 

LNG  VAPOR  DISPFRSION  PREDICTIONS  EOR  25,000 
INSTANTANEOUS  SPILL  ON  WATER  BASED  ON  CLASSICAL  PL 


2.  Some  groups  have  included  effects  due  to  gravity 
spreading  of  the  cold  LNG  vapors;  others  did  not. 

Where  included  (FPC,  Feldbauer) , the  modeling  processes 
were  dissimilar. 

3.  The  values  used  for  Oy  and  az,  the  horizontal  and 
vertical  dispersion  coefficients,  were  not  always  the 
same.  Different  sources  of  these  data  have  been  used 
and  "adjustments"  have  been  made  to  these  data  in  an 
effort  to  more  accurately  reflect  the  expected  LNG 
cloud  behavior.  Finally,  the  predictions  made  have  not 
always  assumed  the  same  meterological  stability  conditions, 
e.g.,  "neutral"  vs.  "stable". 

4.  Modifications  have  been  made  to  account  for  the 
area  nature  of  the  source  (Equation  VII  - 1 describes 
the  dispersion  from  a point  source)  and  the  modeling 
processes  were  dissimilar. 

A description  of  methods  used  by  each  of  the  four  groups 
to  obtain  the  predictions  in  Table  VII  - 1 follows.  The 
calculation  of  vapor  flow  rate,  allowances  for  gravity  spreading, 
selection  and  modification  of  dispersion  coefficients  to  "fit" 
LNG  behavior,  and  allowances  for  the  effect  of  area  sources  are 
described  in  detail.  A description  of  sources  of  dispersion 
coefficient  data  from  which  ail  of  the  groups  selected  some  data 
is  shown  in  Appendix  I. 
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VII-A . PREDICTIONS  USING  BURGESS  * MODEL  (1,  2) 


Burgess 1 model  for  LNG  vapor  dispersion  is  the  classical 
plume  model,  Equation  V - 12,  restated: 


C(x,y,z) 


y z “ 


exp 


hw-u-y] 


(VII  - 2) 


where  Q'  = 


a 


y ' 


0 


z 


u - 


rate  of  LNG  vapor  flow  rate 
downwind 

horizontal  and  vertical  coefficients 
of  dispersion,  respectively 

mean  wind  velocity  in  the  X - direction 


Calculation  of  Vapor  Flow  Rate 

Based  on  data  obtained  from  approximately  steady  spills 
of  LNG  at  rates  of  the  order  of  one  cubic  meter  per  minute, 
Burgess  found  the  maximum  diameter  of  the  LNG  pools  to  be 
given  by 

D - 6.3  W1/3  (VII  - 3) 

where  D = maximum  pool  diameter,  feet 
W = weight  of  LNG  spilled,  lbs 

The  corresponding  evaporation  time  was  found  to  be 

T = 2.5  WX/3  (VII  - 4) 

where  T * evaporation  time,  sec 

W = weight  of  LNG  spilled,  lbs 
Twenty  five  thousand  cubic  meters  of  LNG  weighs  23.4  x 106 
lbs.  Equations  VII  - 3 and  VII  - 4 therefore  give  a maximum 
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pool  diameter  of  1800  feet  and  an  evaporation  time  of 

12  minutes.  This  corresponds  to  an  average  vapor  production 

rate  of  approximately  750,000  ft"*  per  second  (at  70°F,  1 atm). 

The  peak  evaporation  rate  occurs  when  the  pool  size  is 

maximum.  Burgess  used  a steady  LNG  boil  off  rate  of  0.037 
2 

lb/ft  - sec  based  on  his  experimental  results.  The  maximum 
evaporation  rate  was  then  estimated  to  be  about  2,000,000  ft^ 
per  second  (at  70°F,  1 atm) . Burgess  then  treats  the  problem 
as  a steady  release  with  750,000  ftvsec  and  2,000,000  ft  /sec 
as  lower  and  upper  limits  on  the  vapor  flow  rate  Q'. 

Source  of  Dispersion  Coefficients 

Burgess  used  dispersion  coefficient  correlations  proposed 
by  Singer  and  Smith  (Appendix  I).  Singer  and  Smith's  correla- 
tion of  tjy  and  a ^ with  downwind  distance  X can  be  represented 
by  the  equations  shown  below. 

ATMOSPHERIC  DISPERSION  COEFFICIENTS 
FROM  SINGER  AND  SMITH  (21) 

USED  IN  BURGESS'  MODEL 

Gustiness  Classification 

(Meteorological  Stability)  Plume  Dimensions  (ft) 

0.78  0.78 

C (Neutral)  0.42  X 0.29  X 

0.71  n 7i 

D (Stable)  0.44  X 0.087  X 

Burgess  found  that  in  order  to  fit  his  data  from  small  spills 
using  Equation  V. 1 - 2,  the  pronounced  layering  (gravity  spreading) 
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of  the  vapor  which  he  observed  in  his  experiments  had  to  be 
accounted  for.  He  found  that  reasonable  agreement  between 
the  model  predictions  and  his  small  scale  experimental  data 


r t 


was  obtained  when  the  correlations  for  oz 
replaced  by  0.2c  . 

y 


above  were 


Provision  for  Area  Source 

Burgess  makes  no  provision  for  the  area  source  nature 
of  the  spill.  The  predictions  are  made  with  Equation  VII  - 2 
which  assumes  the  vapor  is  released  from  a point  source. 

Burgess’  predictions  of  downwind  distance  to  the  average 
5%  concentration  level  following  instantaneous  release  of 
25,000  M3  LNG  in  a 5 MPH  wind  under  stable  weather  (Singer 
and  Smith  D category)  conditions  are  shown  in  Table  VII  - 1. 


i 
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VII-B . PREDICTIONS  USING  FBLDBAUER  MODEL  (3) 


Feldbauer  (3)  has  published  results  of  spill  tests 
ranging  in  size  from  250  to  2700  gallons  (approximately  1 to 
10  cubic  meters) . Spill  times  varied  from  3 seconds  for  the 
smallest  spills  to  30  seconds  for  the  largest.  The  basic 
model  used  by  Feldbauer  to  describe  atmospheric  dispersion 
is  the  classical  plume  model,  Equation  V - 12,  restated: 


C(x,y ,z) 


exp 


no  a u 
y z 


; ftf  - ; w ] 


(VII  - 5) 
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Calculation  of  Vapor  Flow  Rate 


Downwind  vapor  concentrations  were  monitored  by  hydro- 
carbon detectors  set  in  lines  at  right  angles  to  the  wind 
d ection.  From  concentration  vs.  time  measurements  at  all 
locations  in  a line,  and  from  the  wind  velocity,  the  total 
vapor  flow  rate  past  a line  of  sensors  as  a function  of  time 
was  calculated.  These  data  were  used  to  predict  the  maximum 
vapor  flow  rate  from  the  spill  area.  The  maximum  vapor  flow 
rate  was  then  used  for  Q'  in  Equation  VII  - 5. 

Figure  VII  - 1 shows  Feldbauer's  suggested  correlation 
for  the  maximum  LNG  vaporization  rate  from  an  instantaneous 
release  of  LNG  onto  water.  Figure.  VII  - 2 shows  Feldbauer's 
suggested  correlation  of  maximum  downwind  vapor  flow  vs. 
maximum  LNG  vaporization  rate. 

For  a 25,000  M3  spill  (6,600,000  gal.),  from  Figure 
VII  - 1,  i « 130,400  lb/sec  (3.1  x 10  ft3/sec  at  70°F,  1 atm), 
and  from  Figure  VII  - 2 for  a 5 MPH  wind,  q/W  = 0.2.  There- 
fore, the  maximum  downwind  vapor  flow  rate  from  a 25,000  M3 
instantaneous  spill  is  estimated  to  be 

q = 0.2  (130,400)  ■ 26,080  lb/sec 

= 6.2  x 105  ft3/sec  (at  70°F, 

1 atm) 

It  should  be  noted  that  the  rationale  for  the  downwind  vapor 
flow  rate  being  lower  than  the  evaporation  rate  is  the 
accumulation  of  dense  vapor  over  the  spill  site. 
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MAXIMUM  VAPORIZATION  RATE  w FOR  INSTANTANEOUS  SPILLS, 
FELDBAUER  et  al.  (3) 


- 


WIND  SPEED  MPH 


Figure  VII  2.  EFFECT  OF  WIND  SPEED  FOR 
INSTANTANEOUS  SPILLS, 
FELDBAUER  et  al.  (3) 
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Allowar.  ;e  For  Gravity  Spreading  Effects 

Feldbauer's  gravity  spread  analysis  is  based  on  the 
following  equation  for  the  plume  width  (during  gravity 
spreading)  as  a function  of  downwind  distance  from  the  spill 
point  suggested  by  Fannelop  and  Waldman  (41) . 


(VII 


L •=  2.2 

• — < 

AP  ghLu 

v3  2/3 

X 

1 

P 

u 

where  L = 

plume  width, 

ft 

p = plume  density/  lb/ft J 

A p = difference  between  plume  and  air 

densities 

g = gravitational  acceleration,  ft/sec'1 

h = plume  height,  ft 

u * plume  speed,  ft/sec 

x = distance  downwind,  ft 


By  taking  the  derivative  of  Equation  VII  - 6 the  following 
equation  for  the  rate  of^lateral  (radial)  spread  with  respect 
to  downwind  distance  traveled  is  obtained. 


dL 

dx 


2.2 


“ ** 

V3 

1 

2/3 

L 

AP  gh 

P 

* • 

u 

X 

(VII 


Equation  VII  - 7 was  used  to  predict  the  gravity  spread  of 
the  cloud  as  follows. 
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The  spreading  plume  is  assumed  to  be  uniform  in 
concentration  and  density  and  approximately  rectangular  in 
cross  section.  At  any  cross  section  of  the  plume  the 
total  mass  {vapor  plus  air)  flow  rate  M (lbs  per  second) 
is  given  by 


M 


100  Q' 


hLup 


(VII  - 8) 


where  Q’  = vapor  flow  rate,  lb/sec 

C = vapor  concentration,  volume  % 

h = height  of  plume,  ft 

L = width  of  plume,  ft 

p = density  of  plume,  lb/ft3 

u = plume  velocity 

Solving  Equation  VII  - 8 for  L gives 


100  Q' 

L = (VII  - 9) 

Chup 


In  Equation  VII  - 9 Q'  has  already  been  specified. 

C,  h,  u,  p must  be  determined.  A relation  between  C and  p 
is  developed  assuming  adiabatic  mixing  of  air(70°F,  70% 
relative  humidity)  with  LNG  vapor  at  its  boiling  point  as 
shown  in  the  first  two  columns  of  Table  VII  - 2.  Based  on 
temperature  measurements  made  during  the  tests,  corrections 
were  made  to  the  density  to  reflect  the  addition  >of  heat  to 
the  cloud  due  to  heat  transfer  from  its  surroundings.  These 


89 


■ '•fL^Kr^nr  w-inw^- v 


corrections  result  in  the  density  values  given  in  column  3 
of  Table  VII  - 2. 


TABLE  VII  ~ 2 

CALCULATED  PLUME  DENSITY  AS  A 
FUNCTION  OF  PLUME  CONCENTRATION 
FROM  FELDBAUER  (3) 


adiabatic 

CORRECTED* 

han6  Mol€i% 

I*3 

p x 10 

P X 103 

100 

115.18 

75 

92.76 

92.38 

50 

81.45 

81.15 

30 

75.86 

75.43 

20 

74.28 

73.87 

0 

74.13 

74.13 

* corrected  for  heat  transferred  to  cloud  from  surroundings 
using  experimental  cloud  temperature  data. 
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The  plume  height  h is  estimated  by  calculating  the 
amount  of  vapor  accumulation  over  the  spill,  assuming  its 
shape  to  be  cylindrical,  and  solving  for  h from  the 
relation 

4V 

h = (VII  - 10) 

irD2 

where  V = volume  of  vapor  accumulation  over 
spill 

D = diameter  of  spill 

Based  on  coi  ^lations  derived  from  their  own  data,  the 

3 

diameter  of  a 25,000  M spill  was  determined  to  be  2036 
feet  and  the  volume  of  vapor  accumulated  was  calculated 

8 "5 

to  be  2.1  x 10  ft  (at  LNG  boiling  temperature  and  1 atm). 

Solving  for  h from  Equation  VII  - 10  gives  h “ 66.2  ft. 

Feldbauer  et  al.  then  suggests  multiplying  this  value  by 

0.6  to  account  for  "diffusion  effects".  Thus,  the  initial 

3 

value  of  h for  a 25,000  M spill  is  66.2  x 0.6  * 40  feet. 
This  value  of  h is  assumed  to  remain  constant  throughout 
the  gravity  spread. 

Finally,  u,  the  plume  velocity  is  estimated  by  assuming 
a linear  relation  between  vapor  weight  percent  of  the  cloud 
and  the  percent  of  the  wind  speed  attained.  The  resulting 
non-linear  relation  between  volume  % vapor  and  percent  of 
wind  speed  attained  is  shown  in  Figure  VII  - 3. 
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Equation  VII  - 9 is  then  used  to  calculate  the  plume  width 
(L)  as  a function  of  the  vapor  concentration,  as  shown  in  Table 
VII  - 3.  The  table  is  terminated  at  a vapor  concentration  of 
22.3%  since  the  plume  density  approaches  that  of  the  air  at  that 
concentration,  i.e.,  the  plume  becomes  neutrally  buoyant' at  that 
point  and  the  gravity  spread  phase  of  the  calculation  is  terminated. 


TABLE  VII  - 3 

PLUME  WIDTH  VS.  CONCENTRATION  DURING 
GRAVITY  SPREAD  (FROM  EQ.  VII-9) 


Vapor  Concentration,  Mole  % 
100 
75 
50 
40 
30 

22.3 


Plume  Width,  Ft. 


4036 

4655 

5793 

7610 

10,180 


Equation  VII-7  is  then  used  (Feldbauer  multiplied  this 
equation  by  2/3  in  order  "to  fit  their  data")  to  calculate  the 
relation  between  plume  concentration,  plume  width,  and  downwind 
distance  traveled  during  the  gravity  spread  phase  as  shown  in 
Table  VII-4.  The  gravity  spread  calculation  was  terminated  when 
the  plume  reached  neutral  buoyancy,  where  the  plume  is  predicted 
to  be  10,180  feet  wide  and  40  feet  high. 
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TABLE  VII-4.  RESULTS  OF  INTEGRATION  OF 
EQUATION  VII  - 7 to  DESCRIBE  GRAVITY  SPREAD  OF 
VAPOR  CLOUD  FROM  25,000  M3  LNG  SPILL 


Downwind  Distance, 


Methane  Mole  % 

Cloud  Width,  L,  ft 

X,  ft 

100 

75 

4036 

90.3 

50 

4655 

128 

40 

5793 

280 

30 

7610 

596 

22.3 

10,180 

1200 

Source  of  Dispersion 

Coefficients 

Feldbauer  suggested  the  use  of  the  following  atmospheric 
stability  classifications  for  describing  the  conditions  present 
during  their  test.  Their  report  implies,  but  does  not 
explicitly  state,  that  they  consider  these  conditions  to  be 
generally  representative  of  stability  to  be  expected  over  water. 

Dispersion  Coefficients  Atmospheric-Stability  Coefficients 
Horizontal  Coefficient ,Oy  Gifford  Pasquill  - "C" 

Vertical  Coefficient , oz  Singer  and  Smith  - "D" 

The  API  approach  is  unique  in  that  all  other  predictions  based 
on  classical  plume  models  to  date  have  utilized  the  same  atmos- 
pheric stability  category  for  estimating  horizontal  and  vertical 
dispersion  coefficients  (although  Burgess  did  modify  the 
vertical  dispersion  coefficients  to  fit  his  spill  data) . 
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Provision  for  Area  Source 


Feldbauer's  model  considers  the  atmospheric  dispersion 
of  the  vapor  cloud  to  begin  with  a cloud  _ ,200  feet  wide. 
They  suggest  that  this  "source"  for  the  classical  model. 
Equation  VII-5,  is  too  large  to  be  represented  by  a point 
source.  They  assume  the  source  to  be  represented  by  a number 
of  point  sources  spread  equidistant  along  a line  equal  in 
length  to  the  width  of  the  cloud  resulting  from  the  gravity 
spread  calculation.  Following  this  method,  the  dispersion 
in  this  analysis  was  assumed  to  be  represented  by  11  point 
sources,  separated  by  equal  distances  of  1000  feet.  Each 
point  source  was  assumed  to  emit  the  total  vapor  flow  rate 
obtained  for  the  25,000  M3  spill  (Q ' = 26,000  lb/sec  or 
6.3  x 105  ft3/sec  at  70°F,  1 atm)  divided  by  11.  a schematic 
of  the  arrangement  is  shown  in  Figure  VII  - 4.  The  downwind 
concentration  is  a maximum  on  the  centerline  of  the  center 
source.  This  maximum  downwind  concentration  is  computed  by 
adding  the  contribution  of  all  eleven  point  source  plumes  to 
the  concentration  at  the  given  distance  on  the  centerline  of 
the  center  plume.  The  concentration  on  the  centerline  of  the 
center  plume  at  any  distance  downwind  is  obtained  from  the 
equation 


C(x,y,z=0) 


11  5.73  x 104 


l uoyo2u 


(VII 
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-POINT  SOURCES 


PLUME  CONCENTRATIONS  ARE  ADDITIVE  AT  ANY  POINT. 
MAXIMUM  CONCENTRATION  OCCURvS  ON  CENTERLINE 
OF  CENTER  PLUME. 


Figure  VII-4.  MULTIPLE  POINT  (LINE)  SOURCE 
REPRESENTATION  OF  VAPOR 
SOURCE  FOR  25,000  M3  INSTAN- 
TANEOUS LNG  SPILL 


where  y = 0 for  the  center  plume  and 

y = 1000,  2000,  3000,  4000,  & 5000  feet 

for  the  plumes  on  each  side  of  the 
center  plume 

To  calculate  the  distance  downwind  to  the  (average)  5% 
concentration,  a distance  is  assumed,  c and  o are  read 

y Z 

from  the  Pasquill  "C"  and  Singer  and  Smith  "D"  dispersion 
coefficient  charts  respectively  (Appendix  I)  and  Equation 
VII-11  is  solved  for  C.  This  process  is  repeated,  by  trial 
and  error,  until  the  calculated  downwind  concentration  is 
5%.  The  result  for  a 25,000  spill,  as  shown  in  Table 
VII  - 1,  is  5.2  miles.  The  downwind  distance  to  the (time 
average)  2.5%  concentration  level,  calculated  using  the 
same  procedure,  is  9.5  miles. 


VH-C.  PREDICTIONS  USING  FPC  MODEL  (7) 

The  FPC  predictions  of  LNG  vapor  dispersion  are  also 
based  on  the  classical  plume  dispersion  model.  Equation 
(V  - 2)  restated: 


C(x,y,z) 


cr 


vo  a 

y z 


u 


(VII  - 12) 


Since  the  method  used  by  FPC  to  estimate  the  vapor  flow  rate, 
Q',  depends  on  the  extent  of  gravity  spread,  it  is  expedient 
to  describe  their  handling  of  the  gravity  spread  process  first. 
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Allowances  for  Gravity  Spread  EffectB 

To  calculate  the  pool  size  and  evaporation  time  for  a 
25,000  instantaneous  spill,  the  FPC  staff  (7)  uses  the 
gravity  spread  relations  proposed  by  Raj  and  Kalelkar  (15) : 


EQUATION  FOR  POOL 
RADIUS  (MAXIMUM) 


EQUATION  FOR 
EVAPORATION  TIME 


r 


e 


t 


e 


7.4 


h 


V4 


8.8  Vc 


(VII  - 13) 


(VII  - 14) 


where  re 

Vo 

h 


maximum  pool  radius,  ft 
evaporation  time.  Bee 

3 

volume  of  spill,  ft  LNG 
liquid  regression  rate,  in/min 


The  regression  rate  is  assumed  to  be  one  inch  per  minute, 

2 

■which  is  equivalent  to  a vapor  flux  of  0.037  lb/ft  - sec  or 

2 

a constant  heat  transfer  rate  >f  approximately  30,000  BTU/hr-ft  . 
This  is  consistent  with  evaporation  rates  used  by  other 

3 

investigating  groups.  For  VQ  « 25,000  M and  h ■ 1 in/min 
the  following  values  were  obtained  using  Equations  VII  - 13 
and  VII  ~ 14. 


MAXIMUM  POOL  DIAMETER  - 2511  feet 
EVAPORATION  TIME  - 270  sec  (4 .5  min) 

The  LNG  vapor  from  the  liquid  pool  is  assumed  to  "pile  up" 
in  a cylindrical  volume  over  the  spill.  The  diameter  of  the 
pure  vapor  cylindrical  volume  is  assumed  equal  to  the  maximum 


A 
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liquid  pool  diameter.  The  pure  vapor  is  assumed  to  be  at 
the  LNG  boiling  point,  112  K,  at  atmospheric  pressure.  At 
this  condition  the  specific  volume  of  the  vapor  is  approxi- 
mately 250  times  that  of  the  liquid.  The  height  of  the  pure 


where  h£  - initial  height  of  pure  vapor  cloud 

3 

V0  = volume  of  LNG  spilled,  ft 

re  = radius  of  pure  vapor  cloud,  assumed 
equal  to  the  maximum  liquid  pool 
radius,  ft 

*) 

For  a 25,000  M instantaneous  spill  the  height  of  the  pure 
vapor  croud  initially  formed  is  determined  from  Equation 
VII  - 15  to  be  45  feet. 

The  FPC  staff  assumes  the  pure  vapor  cloud  formed  over 
the  spill  site,  as  described  above,  spreads  out  laterally  due 
to  gravitational  forces.  The  spread  of  this  pure  cloud,  which 
is  assumed  to  remain  pure  during  the  gravity  spread  process, 
is  calculated  using  the  following  equation. 


where  r = cloud  radius , ft 
t *=  time,  sec 

K = constant,  (K  - 2 assumed) 


t 
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g ••  acceleration  of  gravity,  ft/eec 

p ■ density  of  cloud } 

I arbitrary,  but  consistent 
PA  ■ density  of  air  J units 

H » height  of  cloud,  ft 

Equation  VII  - 16  was  proposed  by  Yih  (22)  as  a model  f'**  the 
density  intrusion  phenomenon,  such  as  the  movement  of  a layer  of 
cold  dense  air  Into  warmer  air  (the  movement  of  weather  "cold 
fronts"  is  an  example) . It  was  used  later  by  Fay  to  describe 
the  spread  of  oil  slicks  (23)  and  LNG  (24)  on  water.  It  can 
be  derived  from  physical  first  principles  if  it  is  assumed 
that  the  only  forces  involved  in  the  spread  are  gravitational 
and  inertial  forces,  i.e„  that  surface  tension  and  friction 

2 

forces  are  neglected.  Substitution  of  the  relation  H « V/*r 
into  Equation  VII  - 16  and  integration  with  respect  to  time 
(assuming  V to  be  constant)  gives  a relation  for  cloud  radius 
as  a function  of  time. 


r * 

4Kg 

’ P-PA- 

/4  1/2  (VII  - 17) 

V t 

« 

. P . 

However,  the  total  volume  of  the  cloud  is  assumed  to  be 
increasing  due  to  heat  transfer  from  the  surroundings.  It  is 
assumed  that  the  entire  process  is  to  be  followed  until  th '? 
cloud  density  decreases  to  that  of  the  surrounding  air, 
after  which  time  the  gravity  spread  modeling  process  is  tar  * 
minated  and  atmospheric  dispersion  (Associated  with  atmospheric 
turbulence)  is  assumed  to  dominate.  The  temperature  at  which 
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pure  LNG  vapor  equals  the  density  of  air  is  assumed  to  be 
151  K.  Using  the  ideal  gas  law  it  is  assumed  that 

dT  = dV 

T ~V~  (VII  - 18) 


during  the  expansion  process.  It  follows  that  the  final 
volume  of  LNG  vapor  (at  151  K)  is  related  to  the  initial 
volume  of  LNG  (at  115  K)  by  the  relation 

V = 337  V 

VN  o 

where  V„  ■>  volume  of  pure  vapor  cloud  at  151  K 
(neutrally  buoyant) 

V = volume  of  pure  liquid  at  115  K 

° (boiling  temperature  of  LNG) 

Assuming  that  the  value  for  V in  Equation  VII  - 17  can  be 
reasonably  represented  by 


Vo  + VN 


and  that  the  cloud  density,  which  is  also  changing,  is 
represented  by  the  log  mean  value,  the  solution  of  Equation 

3 

VII  - 17  for  a 25,000  K instantaneous  spill  gives  the 
following  relation  for  the  radius  of  the  spreading  cloud  as 
a function  of  time. 

r = 5550  t 

where  r = radius  of  cloud,  ft  (VII  - 19 

t « time,  sec 
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The  spreading  process  described  by  Equation  VII  - 19  is 
terminated  at  the  time  the  cloud  becomes  neutrally  buoyant 
(T  - 151  K) . The  time  required  for  the  cloud  to  reach 
T « 151  K is  estimated  by  calculating  the  amount  of  heat 
required  to  raise  its  temperature  from  112  K to  151  K and 
dividing  that  amount  by  the  rate  of  heat  transfer  to  the 
cloud  from  the  water  surface  and  the  air  around  the  cloud. 

The  heat  absorption  required  to  raise  the  cloud  to 
neutral  buoyancy  is 

% * Cp  AT 

* (0.5  BTU/lb  R)  (151  K-112  K)  (1.8  K/R) 

= 35.1  BTU/lb  (VII  - : 

The  rate  of  heat  transfer  to  the  cloud  is  estimated  as  the 
sum  of  the  heat  transfer  rates  from  the  water  and  the 
surrounding  air,  Qw  and  respectively. 


® * k,  + Qa 

* KAATW  + hAATa  (VII  - 


where  K = thermal  conductivity,  water 
(3.13  x 10"4  BTU/meter  sec) 

A » area,  cloud  - water  interface, 
5549  r t 


r 1 

1 at2 

ATW  “ 

ata  - AT2  - atJ  j 

' In  

L J / 

ATi 

AT2  * 

273  - 112  - 161  K 

AT.  ■ 

273  - 151  - 122  K 

a = thermal  diffusivity,  water  (1.42  x 10  7 
meter2/sec) 

For  a 25,000  M instantaneous  spill  the  total  rat  of  heat 
transfer.  Equation  VII  - 21,  reduces  to 

Q = 1.32  x 104  t + 2.06  x 105  t1/2 


where  Q = total  heat  transfer  rate  to  vapor 
cloud,  BTU/sec 

t = time,  sec 

The  total  heat  transferred  up  to  the  time  when  neutral 
buoyancy  occurs  (t^)  is 


(VII  - 22) 


/*n  r 4 

Q = Jo  Qdt  = J 1.32  x 104  t + 2 


= 6.6  x 10  t + 1.37  x 103  tN 

N 


06  x 10'  t 
3/2 


1/2] 


dt 

(VII  - 23) 


= 35.1  W (using  Equation  VII  - 20) 

where  Q = total  heat  transferred  to  cloud,  BTU 
W = total  mass  of  cloud,  lb 

Solving  Equation  VII  - 23  for  t gives  a time  to  neutral 

N 

buoyancy  of  60  sec.  From  Equation  VII  - 19  the  diameter  of 
the  cloud  at  the  time  when  it  becomes  neutrally  buoyant  is 
then  3785  feet. 

Summarizing,  the  condition  of  the  cloud  at  the  end  of 
the  gravity  spread  process  is  estimated  to  be  as  shown  in 


Table  VII  - 5. 


TABLE  VII  - 5.  VAPOR  CLOUD  DESCRIPTION  AT  END 
OF  GRAVITY  SPREAD  PROCESS  - 25,000  M3  SPILL  (FPC) 

CLOUD  DIAMETER  = 3785  ft 

CLOUD  HEIGHT  = 28  ft 


CLOUD  COMPOSITION  = 100%  LNG  vapor 

CLOUD  TEMPERATURE  = 151  K 


103 


J 


I 

! 


i 


Calculation  of  Vapor  Flow  Rate 

The  FPC's  method  of  calculating  the  value  of  the  vapor 
flow  rate  for  use  in  the  classical  plume  dispersion  equation 
is  unique.  They  assume  that  the  pure  vapor  cloud  which 
exists  at  the  end  of  the  gravity  spread  process  (see  Table 
VII  - 5)  will  release  vapor  from  its  upper  surface  at  a 
rate  determined  by  the  rate  at  which  heat  is  absorbed  by  the 
(now  neutral)  cloud  from  the  surrounding  air.  The  release 
' rate  is  calculated  from  the  following  relation: 


h A 

Q'  » 

C 


(VII  - 24) 


where  h * heat  transfer  coefficient,  air  to 
cloud  (2.99  x 10“3  BTU/m2  sec  F) 

• 

A - area  of  top  surface  of  neutrally 

buoyant  cloud  (*^2  * 1.13  x 10'  ft2) 

C * average  sensible  heat  capacity  of 
cloud  (0.5  BTU/lb  F) 

From  Equation  VII  - 24  the  vapor  flow  rate  is  calculated 


to  be  6250  lb/sec. 


Source  of  Dispersion  Coefficients 

The  FPC  staff  estimated  dispersion  coefficients  from 
the  correlations  presented  by  Gifford  and  Pasquill  (see 
Appendix  I) . The  data  presented  in  the  charts  of  Appen- 
dix I have  been  reduced  to  analytical  equation  form  and 
programmed  in  a computer  subroutine  by  Zimmerman  and 
Thompson  (25) . 
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The  values  of  the  horizontal  and  vertical  dispersion 
coefficients  are  determined  from  the  following  equations. 


0y  = 465.1  (X  + Xy)  tan  |c-[d  In  (X  + Xy) ]j  (VII  - 25) 

o = ax*  L V J 


(VII  - 26) 


where  a 

y 

a 

z 

X 


X 


v 


horizontal  dispersion  coefficient,  meters 
vertical  dispersion  coefficient,  meters 
downwind  distance,  meters 

upwind  distance  to  virtual  source,  meters 


a,  b,  c,  d * constants  derived  from  curve  fit  of  Gifford 
Pasquill  charts  (Appendix  I) 


The  values  for  c and  d are  functions  of  stability  class 
only.  The  values  of  a and  b are  functions  of  distance  as 
well  as  stability  class.  Values  of  a,  b,  c,  d are 
reproduced  from  Zimmerman  (25)  in  Tables  VII  - 6,  VII  - 7 
and  VII  - 8. 


TABLE  VII-6.  VALUES  OF  a AND  b USED  IN 
EQUATION  VII  - 26  FOR  D STABILITY  CLASS 


Downwind 
Distance  (km) 

a (meters) 

b (dimensionless) 

0.3-1 

32.093 

0.E1066 

1-3 

32.093 

0.64403 

3-10 

33.504 

0.60486 

10-30 

36.650 

0.56589 

>30 

44.053 

0.51179 
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TABLE  VII  - 7 . VALUES  OF  a AND  b 
FOR  USE  IN  EQUATION  VII- 26 
FOR  F STABILITY  CLASS 


Downwind 


Distance  (km) 

a (meters) 

b (dimensionless) 

0.2-0. 7 

14.457 

0.7841 

0. 7-1.0 

13.953 

0.6847 

1-2 

13.953 

0.6323 

2-3 

14.823 

0.5450 

3-7 

16.187 

0.4349 

7-15 

17.836 

0.4151 

15-30 

22.651 

0.3268 

30-60 

27.074 

0.2744 

>60 

34.219 

0.2172 

TABLE  VII  - 8. 

VALUES  OF  C AND 

d FOR  USE  IN  EQUATION  VII-25 

Stability  Class 

c (degrees) 

d (degrees) 

D - Neutral 

8.333 

0.72382 

F - Stable 

4.167 

0.36191 

It  should  be  noted  that  the  FPC’s  published  predictions 
of  downwind  distance  to  the  lower  flammable  limit  <5% 
average)  have  been  based  on  the  assumptions  of  D-Neutral 
stability  meteorological  conditions. 
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Allowance  for  Area  Source  Effects 


The  FPC  model  accounts  for  the  area  source  nature  of 

an  evaporating  LNG  pool.  The  area  source  is  treated  as 

a virtual  point  source  located  at  a distance  *upwind  of 

the  spill  which  corresponds  to  a horizontal  standard 

deviation  a given  by  the  relation: 

*o 

o = D'/4.3  (VII-27) 

yo 


where  av  = standard  deviation  at  spill  site 
* o 

equivalent  to  area  source  width 

D'  = width  (diameter)  of  area  source  at 
spill  site 

Equation  VII  - 27  effectively  treats  the  area  source  as 

a cross  wind  line  source  with  a normal  distribution,  and 

was  suggested  by  Holland  (26)  and  Turner  (27) . 

For  a 25,000  M3  instantaneous  spill,  the  FPC  estimate 

for  D'  is  1154  meters  (3790  ft).  Equation  VII  - 27  then 

gives  a value  of  av  = 268  meters  (880  ft).  From 

jo 

Equation  VII  - 25  the  virtual  distance,  Xy, is  determined 

to  be  X = 4.0  km  for  D-Neutral  conditions, 
v 

Applying  the  classical  plume  dispersion  equation 
(Equation  VII  - 12)  to  the  centerline  condition  (y=0) , at 
an  effective  emission  height  H, 


C (x,y=0 , z=H) 


Q' 


™yaz  U 


exp 


(VII-28 ) 
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The  effective  emission  height  is  determined  from  the  relation 


VAvg  / 1Tr  Avg 


* 10.1  meters 


(V1I-29) 


where  V = time  averaged  volume  of  pure  cloud 
Avg  during  gravity  spread, equals  7.39  X 
106  M3 

rA  = (577  + 383)/  2 - 480  M (average  of 

g initial  and  final  gravity  spread 
radii 

Substituting  values  for  H and  Q'  determined  above » with 
u = 5 MPH  (2.24  M/sec)  into  Equation  VII  - 28  and  utilizing 
the  relations  for  a and  oz  given  in  Equations  VII  - 25 
and  VII  - 26  and  Tables  VII-6,  VII-7,  and  VII-8  for  the  D 
stability  class,  the  following  relation  is  obtained: 


2,838,000 


1 r 4. 

2 L32- 


093  X °*644 


1 2 ’ 


(VII-30) 


ir  (465.1)  (x+4 .0)  tan  ^8 . 333-  (0 . 7238  In  (x+4 . 0 ))](32 . 09)  X°  * 644  (2 . 24 ) 
where  C is  in  gm/M3 

By  trial  and  error,  the  solution  of  Equation  VII-30  for  X, 

the  downwind  distance  to  the  average  5%  concentration  level 

3 3 

(36.6  gm/M  ) following  instantaneous  release  of  25,000  M of 

LNG  during  D-Neutral  weather  conditions,  is  found  to  be  1.2  km 

or  0.75  miles  as  shown  in  Table  VII-1. 


VII-D.  COMPARISON  OF  RESULTS 

BASED  ON  STEADY  RELEASE  MODELS 

The  downwind  distances  to  the  time  average  5%  concentration 

level  calculated  for  a 25,000  M3  instantaneous  spill  using  the 
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models  proposed  by  Burgess  (1,2),  Feldbauer  (3)  and  the  FPC 
staff  (7)  are  plotted  in  Figure  VII-5  as  a function  of  vapor 
release  rate  used  in  the  predictions.  The  largest  predicted 
distance,  50  miles,  obtained  by  Burgess  using  a vapor  flow 
rate  equal  to  his  predicted  peak  evaporation  rate  is  almost 
70  times  greater  than  the  distance  of  0,75  miles  predicted 
by  the  FPC  staff.  The  downwind  distances  calculated  using 
Burgess ' model  with  a vapor  flow  rate  equal  to  his  predicted 
average  evaporation  rate  and  with  Feldbauer' s model  lie  in 
between . 

Burgess'  model  does  not  account  for  area  source  effects 
or  effects  due  to  gravity  spreading  immediately  following  the 
spill.  Furthermore,  his  predictions  for  the  "worst  case" 

25,000  M3  instantaneous  spill  assume  very  stable  meteorological 
conditions.  Burgess  used  linger  and  Smith's  dispersion 
coefficients  for  the  D-gustiness  category  which  are  a close 
approximation  to  the  most  stable  weather  category  (F)  of 
Pasquill  (see  Appendix  1) . The  uppermost  line  in  Figure 
VII-5,  drawn  through  Burgess'  predicted  values,  therefore, 
represents  a "worst  case"  downwind  distance  to  the  5% 
concentration  level  as  predicted  by  the  classical  point  source 
steady  plume  dispersion  model.  The  extreme  effect  on  these 
predicted  distances  of  the  values  used  for  the  dispersion 
coefficients  a and  o is  seen  when  the  same  calculations  are 

y z 

carried  out  for  weather  conditions  described  by  Burgess  as 
B2~gustiness  classification  (representative  of  unstable 
meteorological  conditions) . The  lower  line  of  Figure  VII-5 
represents  the  downwind  distance  to  the  5%  concentration 


BURGESS 


© API 


BURGESS  B2  GUSTINESS  (UNSTABLE) 


level  calculated  using  the  Burgess'  model  for  B2“gustiness 
category  classification  with  the  values  of  «2  set  equal  to 
0.2  <r  co  better  describe  the  vertical  dispersion  of  the 
dense  LNG  vapors. 

All  of  the  predictions  of  the  downwind  distance  to 
the  5%  concentration  level  fall  inside  these  two  lines, 
and  the  reasons  for  the  different  values  predicted  are 
indicated  by  the  location  of  the  particular  prediction  in 
relation  to  these  two  "bounding  cases". 

The  prediction  by  Feldbauer  of  a downwind  distance 
of  5.2  miles  can  be  attributed  to  two  factors.  First, 
the  estimate  of  a much  lower  vapor  flow  rate  due  to  the 
assumption  of  accumulation  of  the  vaporized  LNG  over  the 
spill  site  leads  to  a shorter  distance.  Secondly,  the 
treatment  of  the  vapor  source  as  a line  source  almost 
2 miles  wide  markedly  reduces  the  downwind  distance  below 
that  whic  would  be  predicted  using  a point  source.  Since 
this  line  source  width  results  directly  from  their  treatment 
of  the  gravity  spread  phase,  the  API  allowance  for  gravity 
spread  is  a strong  factor  in  the  shorter  predicted  distance. 

It  might  be  expected  that  API's  use  of  atmospheric  stability 
category  Pasquill  - C in  the  horizontal  and  Singer  and  Smith 
D in  the  vertical  direction  would  result  in  a much  shorter 
distance  than  would  have  been  obtained  if  the  Singer  and 
Smith  horizontal  stability  category  D had  been  used.  However, 
this  .is  not  the  case.  Calculations  were  made  to  determine 
the  difference  in  downwind  distance  which  would  be  obtained 
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using  the  Singer  and  Smith  D stability  category  for  determining 
the  horizontal  as  well  as  vertical  dispersion  coefficients. 

The  distance  was  calculated  to  be  slightly  shorter  than  5.2 
miles.  This  surprising  result  can  be  explained  by  referring 
to  Figure  VII-6.  Figure  VII- 6~a  is  a schematic  representation 
of  the  additive  nature  of  the  point  sources  representing  the 
10,400  feet  line  source  previously  described,  using  the 
dispersion  coefficients  suggested  by  Feldbauer.  Figure  VII- 6-b 
is  a schematic  representation  of  the  additive  nature  of  the 
point  sources  representing  the  10,400  feet  line  source,  using 
horizontal  (as  well  as  vertical)  dispersion  coefficients 
representing  Singer  and  Smith  D stability  category.  Since 
the  horizontal  dispersion  of  the  individual  point  source 
plumes  is  reduced,  the  plumes  to  either  side  of  the  center 
plume  contribute  less  to  the  center  plume,  and  the  downwind 
distance  along  the  centerline  plume,  which  is  the  maximum, 
is  correspondingly  reduced.  Hence,  Feldbauer 's  predicted 
distance  of  5.2  miles  should  properly  be  attributed  to  the 
lower  vapor  production  rate  and  the  large  .gravity  spread 
effect. 

The  smallest  downwind  distance  to  the  5%  concentration 
level,  0.75  miles  predicted  by  the  FPC  staff,  can  be 
attributed  primarily  to  two  factors.  First,  the  low  value 
utilized  for  the  vapor  flow  rate,  143,000  ft'Vsec  (70°F,  1 atm), 
is  the  primary  reason  tor  the  short  distance  predicted. 
Secondly,  the  use  of  Pasquill  D stability  category  dispersion 
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A OFF  CENTER  PLUMES  CONTRIBUTE  TO  CENTERLINE 
CONCENTRATION  - PASQUILL  "C"  HORIZONTAL, 
"F"  VERTICAL  STABILITY  CATEGORY 


B - OFF  CENTER  PLUMES  CONI  RIBUTE  LESS  TO  CENTER 
LINE  CONCENTRATION-SINGER  AND  SMITH  "D" 
STABILITY  CATEGORY 


Figure  VII  6.  SCHEMATIC  REPRESENTATION  OF  ADDITIVE 
PROPERTY  OF  POINT  SOURCES  REPRESENT 
ING  A LINE  SOURCE 


coefficients  rather  than  the  "worst  case”  F coefficients 
Iso  contributes  to  the  shorter  distance.  The  correction  for 
the  area  nature  of  the  source  resulted  in  a less  important 
reduction  in  the  predicted  downwind  distance. 
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VIII.  VAPOR  DISPERSION  PREDICTIONS  BASED  ON  SOLUTION 
OF  COMBINED  ENERGY,  MOMENTUM,  AND  MASS  BALANCE  EQUATIONS  - 
SCIENCE  APPLICATIONS,  INCORPORATED  MODEL 


Science  Applications,  Inc.  (SAI)  has  made  predictions  of 
dispersion  of  LNG  vapor  from  large  LNG  spills  on  water  in  a 
series  of  risk  assessment  studies  done  for  Western  LNG  Terminal 
Company  (8).  SAI's  approach  involves  solution  of  the  system 
of  equations  representing  the  accountability  of  mass,  momentum, 
and  energy  associated  with  an  LNG  spill.  Equations  VIII-1, 
VIII-2,  and  VIII-3  are  balance  equations  for  mass,  momentum, 
and  energy  respectively,  restated  as  follows: 

Accountability  of  Mass 

£P  = - V ♦ pv  (VUI-la) 

3 1 


3C  = - V-Cv  (VUI-lb) 

3t 

Accountability  of  Momentum 

3Cv  = - V'pW  - V * t + pg  (VIII-2) 

3t 

Accountability  of  Energy 

9 pH  = - V * pH v - V*q  + DP  - t:Vv  (VIII-3) 

9t  Dt 


i i 


i j 
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.»■**> 


!'  t.» 


where  p 


density  of  gas-air  mixture 


H = enthalpy  (energy  content)  of  gas-air 
mixture 

v =■  velocity  vector,  decomposable  into 
components  u,  v,  w 

q = heat  transfer  vector,  decomposable 
into  components  qx,  qy,  q2 

$ 

t “ stress  tensor,  decomposable  into  9 
components 


T:cx  • Txy ' Txz  > TyX ' t yy  / xyz » Tzx  • T zy  * Tzz 
P = pressure 

t * time 

g - gravity  force  vector,  decomposable 

into  components  gx  = 0,  gy  *»  0,  gz  = 

32.2  ft/sec2 


J2_ 

Dt 


substantial  derivative  operator, 


3 „ 

— + v.  V 

at 


Solution  of  Equations  VIII-  1,  2,  3 with  appropriate 
boundary  conditions  describing  the  LNG  vapor  source,  the  air 
temperature  and  humidity,  and  the  heat  transfer  between  the 
gas-air  mixture  and  its  surroundings  should  provide  a complete 
description  of  the  vapor  cloud  development  and  dissipation. 

The  following  section  describes  SAI's  simplification 
of  Equations  VIII-  1,  2,  3,  assignment  of  boundary  conditions, 
and  specification  of  input  data. 
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Accountability  of  Mass 

Neglecting  molecular  diffusion  in  Equation  VUI-lb  and 
expanding: 


3 P = 

r 3 pu  + 

3 pv 

+ 3 pw  I 

(VIII-4a) 

It 

ay 

17  J 

3C  = 

lu3C  + v 3C  + w 

3cl 

- C [ 3u  + 3v  + 3w j 

It 

^ 3x 

3y 

3zJ 

[ 3x  3y  3z J 

(VIII-4b) 

Accountability  of  Momentum 


Equation  VIIX-2  is  expanded,  with  vertical  accelerations 
and  viscous  forces  neglected  in  the  equation  for  accountability 
of  vertical  (Z)  momentum: 


3pu 

“It 


9pv 

“It 


- I 3 puu  + 3 puv  + 3 puwl 

j lx  ly  Tz  J 

- f a * + _i  t + _2  tzx  . - 3p 

L 3x  17  y 3z  J 3x  (VIII-5a) 

- I 3 puv  + 3 pw  + 3 pvw  I 

[n  17  II  J 

” J*  Txy  + Tyy  4 Tzy  I ~ 

[lx  1 Iyyy  H yJ  17  (VIIl-5b) 
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- 3P  + pg 
3z 


(VIII-5C) 


Accountability  of  Energy 

The  energy  balance  is  simplified  by  neglecting  viscous 
dissipation  (heating  due  to  fluid  friction)  and  heat  trans- 
fer by  conduction  (heat  transfer  due  solely  to  temperature 

i -*■ 

gradients),  represented  by  the  terms  x :Vv  and  V?q 
respectively,  in  Equation  VIIX-3. 


3 pH  - - V • pHv  + DP  (VIII-6) 

3t  Dt 

The  fluid  motion  to  be  described  is  turbulent. 

Following  standard  practice,  the  variables  velocity, 
enthalpy,  concentration,  and  pressure  are  expressed  as 
the  sum  of  a mean,  or  time  averaged,  component  and  an 
instantaneous  deviation  from  that  mean  value,  as  follows: 

-*•  + -♦ . 
v *5  v + v 

u *=  u + u* 

V * V + v* 

w ■ w + w’ 

(VIII-7) 

H •=  H + H' 

C - C + C’ 

p = i>  + P’ 

Substituting  Equation  VIII-7  into  Equations  VIII-4,  5,  6, 
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taking  time  averages,  and  noting  that  v'  =u'  * v'  = w'  «H' 


3 pu 
"5x 


+ 3 pv  + 3 pw  | 

Jy  J 


- ju3C  + v3C+w3c|  - V-(C'v*; 
|_  3x  3y  3z  J 

- I 3 puu  + 3 puv  + 3 puw 

[a5  Sy  3i  J 


- I 3 pu ' u ' + 3 pu* v'  + 3 p'uTwT 

j^3x  Tiy  3z 

' [it  T*x  + 4 Tyx  + ^ Tzx] ' 


[3  puv  + 3 pvu  + 

15  Ty 


5 pvw 

75 


[5  p v 1 vi r + 3 pvYv'  + 3 pv ' w ' j 

7x  3y  75  J 

[3  T„v  + 3 Tyy  + 3 TZy  ” 3P 

75  xy  7y  yy  75  yJ  77 


(VIII-8a) 


(VIII-8b) 


(VIII-9a) 


i 

1 

t 


r 


0 = - 3_P  + pg 

3z 


(VIII-9C) 


3 pH  = - V-pHv  - VpH'v'  + DP  (VIII-10) 

3t  Dt 


where  the  subscript  1 in  Equations  VIII-9a,  9b, 
denotes  "laminar"  shear  stresses. 


The  ideal  gas  equation  of  state  is  used  to  relate  density 
and  temperature: 


P 

p = — r (VIII-11) 

nRT 


where  p 
P 
R 
T 
n 


* density  of  gas-air  mixture 
= pressure 

= ideal  gas  constant 
= temperature 

= moles  of  gas  mixture,  jl-C  + C | 

L.Ma 

* molecular  weight  of  air,  29 

= molecular  weight  of  methane,  16 


This  formulation  assumes  the  pressure  is  equal  to  the  sum 
of  the  partial  pressure  of  air  and  methane  and  neglects  any 
contribution  by  water  to  the  pressure. 


\ 

\ 

i 
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The  enthalpy  of  the  mixture  of  water  vapor  (or  ice) , 
methane  and  air  is  assumed  given  by  the  expression: 

H = jca(l-C)  + CmcJ  T + WL0  ( 1-C)  f (T)  (VIII-12) 

where  Ca  = heat  capacity  of  air,  0.24  cal/gm  C 

= heat  capacity  of  methane,  0.52  cal/gm  C 
W = mixing  ratio  of  water  vapor  in  air 

Lq  = latent  heat  of  condensation  and  freezing 
of  water  (675  cal/gm) , assumed  to  occur 
over  a temperature  range  of  -1  C to  1 C 

f (T)  = a linear  function  representing  the 
temperature  dependence  of  the  phase 
transition 


The  system  of  Equations  VTII-8,  9,  10,  11,  12  cannot 
be  solved  without  relating  the  terms  involving  the  velocity, 
concentration  and  enthalpy  deviations  (the  primed  quantities 
in  SquaA  ns  VIII-8,  9,  10)  to  the  mean  values  of  those 
quantiti  , This  is  known  as  the  "closure"  problem  of 
turbulence  modeling.  The  simplest  form  of  "closure"  which 
has  been  ore  osed  (the  so-called  First-Order  closure)  assumes 
that  the  product  of  the  deviation  variables  are  proportional 
to  the  gradient  of  the  associated  mean  values  of  the  same 
variables . This  method  is  used  by  SAI  in  their  model  for 
LNG  vapor  dispersion.  Specifically,  the  following  relations 
are  assumed: 


pu'u'  = Kxl  3u 
3x 

pu'vf  = Kx2  3u 

3y 


(VIII-13-a) 
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pu*  w1 

(VIIl-13-a) 

pv'u* 

= K-H 

pv'v' 

- KxS  |2 
3y 

p^' w' 

“ Kx6  |2 
9z 

pH ' U 1 

= pcT'u*  = 

- kx 

9T 

3x 

pH  ' V 1 

= pcT'v'  = 

- ky 

9T 

3y 

(VIII-13-b) 

pH  'wT 

= pcT'w'  = 

- kz 

3T 

u 1 

'C' 

— - K' 

ac 

X 

9x 

(VIII-13-c) 

V' 

'C* 

= - K ' 

9C 

y 

9y 

- 

w1 

rcr 

= - K' 

9C 

z 

9 z 
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Kx5  = KH 


(VIII-14-a) 


and  Kxl  = Kx2  » Kx4  - 


Kx3  = Kx6  = ^ 


^x 


k 


H 


k 

z 


= *V 


K' 

x 


K' 

z 


% 


(VIII-14-b) 


(VIII-14-c) 


9 


where 


K = "eddy  viscosity" 

k = "eddy  thermal  conductivity" 

K'=  "eddy  diffusivity" 

c = heat  capacity 

subscript  H denotes  horizontal 

subscr !.pt  v denotes  vertical 


Substituting  Equations  (VIII-13,  14)  into  Equations  (VIII-8, 
9,  10)  and  neglecting  laminar  shear  stresses,  the  following 
equations  result. 


9C  = 

at 


_a  pu  + _a_  pv  + _a  pw  | 

3x  3y  3z 

u ac  -i-  v ac  + w aC 
ax  3y  3zJ 

__a  k ' ac  + _a  ka  ac  + _a  k, 

3v  **  ^ vr  gy  gy 


(VIII-15-a) 


(VIII-15-b) 
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puu  + _3_  puv  + _3  puw 
3y  3z 


(VIII-16-a) 


[_3^  R 3u  + _3  K 3u  + 3 K-  3u 

3x  3x  3y  3y  3z  3zJ 

[3  pvu  + 3 pvv  + 3 pvw 

lx  ly  Tz  J 

[3  K 3v  + 5 K 3v  + 3 K 3v 

3x  H 3x  3y  H 3y  3z  v 3z  J 


(VIII-16-b) 


3P  = - pg 

3z 


(VIII-16-c) 


3H  = 
3t 


pHu  + _3  pHv  + __3  pHw 
3y  3z 


T 3 k 

[3x 


„ 3T  + _3  k 3T  + 
H 3x  3y  H 3y 


_i  kv 

3z  3zJ 


+ DP 
Dt 


(VIII-17) 


Equation  VIII-17  is  further  simplified  by  assuming 
V*v  = 0 (neglecting  compressibility  of  the  gas-air  mixture) 


Then 


V * pHv  = vVpH  + pHV*v 
= v’VpH 
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and  Equation  VIII-17  becomes 


9H  = 

at 


(VIII-19) 


where  r , the  adiabatic  lapse  rate  (vertical  temperature 
gradient  for  a neutrally  stable  atmosphere) , is 
included  to  insure  that  a "neutral"  atmosphere 
is  not  perturbed  by  the  turbulent  diffusion. 

Since  the  hydrostatic  approximation  (Equation  VIII-16-c) 
provides  a relation  between  the  pressure  and  altitude,  it  can 
be  used  to  transform  the  proceeding  equations  so  that  pressure 
is  an  independent  variable  and  altitude  is  a dependent 
variable.  furthermore,  a dimensionless  pressure,  ° , can  be 
defined  as  follows: 


a 


P-PT  _ P-Pt 


IT 


(VIII-20) 


where  o - dimensionless  pressure 
P = mean  local  pressure 

Ps  = pressure  at  earth's  surface,  may  depend 
on  position  and  time 


PT  = pressure  at  upper  boundary  of  atmospheric 
region  being  considered 
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Transformation  of  Equations  VIII-15,  16  and  19  to 
the  x,  y,  a , t coordinate  systems  gives  the  following  system 
of  equations  to  be  solved 


3tt  + _3_  (ttu)  + __3  (tfv)  + _3  (iro)  = 0 
3t  3x  3y  3o 


(VII1-21) 


Du  + 3_<fc  + a = g^p  8 8_u 

Dt  3x  p 3x  r2  3o  3o 


+ -J.  xh  in'*  + _2  k 3^1 

3x  3xJ  3y  3yJ 


Dv  + 3$  + o5tt  = 
Dt  3y  p 3y 


= 2?J±  Jl  pXt7  ix 

-rr 2 ^ [ 3oJ 

+ _3  KH  3v  + _3  K 3v 

3x  ' 3xJ  3y  [ 3yj 


(VIII-22) 


DH  = DP  + £p  _3_  [kv  [gp  3T  - rl 

Dt  Dt  it  3o  I 3o  ! 

+ _3  I k 3T I + _3  Ik  3T 

3x  [ H 8xJ  3y  [ H 3yJ 


(V1II-24) 
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dc  = ais.  _JL  [\\  1£  + _1 

Dt  9o  L v 3°J  3x 


(VIII-25) 


air  + PT 
P = 

_ (VIII-26) 

+— 1 
Mm  J 

(1-C)  + r C T + WL_ ( 1-C) f (T)  (VIII-27) 

a m j o 

where  4>  = gz 

and  D_  = + u + v _3  + (substantial  derivative 

Dt  3t  3x  3y  3a  operator  in  x,  y,  c , t 

coordinate  system) 

The  "mathematical'*  problem  of  LNG  vapor  dispersion 
consists,  therefore,  of  solution  of  the  set  of  Equations 
VIII-21  through  VIII-27  with  appropriate  boundary  conditions, 
using  finite  difference  (digital  computer)  methods.  A 
Circular  LNG  spill  shape  is  assumed.  A three  dimensional 
(x,  y,  a ) qrid  system  is  laid  out  to  enclose  the  volume  of 
the  atmosphere  into  which  the  LNG  is  evaporated  from  the 
spill.  Due  to  the  symmetry  of  the  assumed  spill  and  the 
resulting  symmetry  of  the  dispersion  process.,  only  half  of 
the  vapor  cloud  development  need  to  be  described.  A reflective 
boundary  condition  is  therefore  incorporated  at  a vertical 
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plane  through  the  center  of  the  spill  coincident  with  the 
wind  direction.  The  grid  system  is  illustrated  in  Figure 
VIII-1.  Figure  VIII-1  also  illustrates  the  type  of 
boundary  conditions  applied  to  the  boundaries  of  the  grid 
system  (28) . 

As  stated  previously,  SAI  assumes  that  turbulence 
associated  transfers  of  mass,  momentum  and  energy  (heat) 
are  proportional  to  the  local  gradients  in  mean  concentration, 
velocity,  and  temperature,  respectively,  in  the  flow  field. 

An  immediate  requirement  is  specification  of  the  eddy 
transfer  coefficients: 

= horizontal  "eddy  diffusivity" 

K'  = vertical  "eddy  diffusivity" 

Kh  = horizontal  "eddy  viscosity" 

Ky  = vertical  "eddy  viscosity" 

kH  = horizontal  "eddy  thermal  conductivity" 

= vertical  "eddy  thermal  conductivity" 

SPECIFICATION  OF  EDDY  TRANSFER  COEFFICIENTS 
The  key  process  in  SAI ' s specification  of  the  eddy  transfer 
coefficients  is  the  method  of  specifying  the  value  of  the 
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3u/3o  = 0 
8d/3 a - 0 
o = 0 


T = T» 


ALONG  THIS  BOUNDARY,  a 


0, 


ALONG  THIS  BOUNDARY,  FOR  x,y  INSIDE  POOL  AREA,  a = v = 0 

T = 204  R,  cr  = 1,  ENTHALPY  FLUX  = mH  VAPOR 
m = 0.04  Ib/ft2  SEC 

FOR  x.y  OUTSIDE  POOL  AREA,  u = i7  = <J  = 0 
T _ Twater,  a ~ 1,  q = hA(^T) 


ALONG  THIS 
BOUNDARY 

u = u (o  ) 
v = 0 

oj  - 0 

T = TH20 


SPILL 

POOL 


r(T) 


X ALONG  THIS  BOUNDARY 

3T/3y  = 0 
i.-C/Dy  = 0 
dx/dy  = 0 


AT  TIME  t - 0,  T = T0(cr),c  = 0,  a =u  (a),  v - a>  = 0 


Figure  Villi.  DESCRIPTION  OF  SAI  "SIGMET"  GRID 
SYSTEM  AND  BOUNDARY  CONDITIONS 
USED  FOR  LNG  VAPOR  DISPERSION 
PREDICTION 
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vertical  "eddy  viscosity  coefficient",  Kv.  All  of  the 
remaining  eddy  transfer  coefficients  are  determined  from 
the  values  assigned  to  Kv. 

The  prediction  of  the  vertical  eddy  viscosity  coeffi- 
cients Ky  by  SAI  is  based  on  a method  proposed  by  Hanna  (29) . 
This  method  assumes  that  the  vertical  mixing  efficiency  of 
the  atmosphere  (which  is  quantified  by  the  value  of  Kv)  is 
dependent  on  the  mean  eddy  sizes  and  the  amounts  of  turbulent 
energy  carried  by  the  eddies.  Since  the  eddy  sizes  and 
amounts  of  turbulent  transfer  associated  with  eddy  movements 
are  related  to  the  energy  spectrum  of  the  vertical  fluctua- 
tions of  the  wind  speed,  it  is  hypothesized  that  the  eddy 
viscosity  should  be  dependent  on  the  characteristics  of 
the  vertical  velocity  spectrum  of  the  atmosphere.  Hanna 
assumes  that  the  vertical  velocity  spectrum  can  be  completely 
determined  by  two  quantities;  the  standard  deviation  of  the 
vertical  fluctuations  of  the  wind  velocity , ffw,  and  the  wave 
number  at  which  the  amount  of  vertical  turbulent  energy  is  a 
maximum,  1^.  Based  on  these  arguments,  Hanna  proposed  the 
following  relation. 

Kv  = Cx  ow  km  (VIII-28) 

For  nearly  neutral  conditions  near  the  ground,  Lumley  and 
Panofsky  (30)  proposed  that 


K 


0.4  u*  Z 


(VIII-29) 


1.3  u* 


(VIII-30) 


k = 0.3 

m 

z 


where  u*  = friction  velocity 
z = vertical  distance 


Assuming  Equation  VII-29,  30,  31  along  with  Equation  VII-28, 
the  constant  in  Equation  VII-28  is  0.09: 

K = 0.09  ^ 

v (VIII-32) 

Taylor  et  al . (31)  have  reported  the  following 
correlation  between  atmospheric  turbulence  scale  length  L 
and  ^ derived  from  spectra  of  vertical  air  velocity  measured 
from  aircraft  at  heights  between  10  and  1300  meters. 

L km  = 0.216  (VIII-33) 

SAI  assumes  L km  = 0.20.  Incorporating  this  expression  into 


Equation  VII-32, 


- 0.45  ow  L 


= 0.45  u <rt  L 


where  u = local  mean  velocity 

at  = standard  deviation  of  wind  direction 


Using  the  data  of  Taylor  et  al.  (31)  SAI  proposed  the  correla- 
tion shown  in  Table  VIII-1  to  describe  the  dependence  of  scale 
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length  on  vertical  distance,  with  atmospheric  stability  category 
as  a parameter. 


TABLE  VII 1-1 

CORRELATION  OF  TURBULENCE  SCALE 
LENGTH,  L,  WITH  HEIGHT  AND  ATMOSPHERIC 
STABILITY  PROPOSED  BY  SAI 


HEIGHT  (meters) 


10 

A 

18 

B 

15 

20 

30 

25 

30 

41 

34 

50 

62 

52 

75 

84 

71 

100 

105 

85 

STABILITY  CATEGORY 


C 

D 

E 

F 

G 

12 

10 

8 

7 

6 

21 

18 

16 

14 

12 

29 

25 

22 

20 

17 

44 

39 

35 

31 

27 

60 

52 

48 

43 

37 

74 

64 

60 

54 

46 

SAI  proposed  the  correlation  shown  in  Table  VIII-2 
between  the  standard  deviation  of  the  wind  direction, ^ , and 
vertical  distance,  with  atmospheric  stability  category  as  a 
parameter. 
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TABLE  VI I 1-2 

CORRELATION  OF  STANDARD  DEVIATION  OF 
WIND  DIRECTION,  <Jt  , WITH  HEIGHT  AND  ATMOSPHERIC 
STABILITY  PROPOSED  BY  SAI 


STABILITY  CLASS 
A 
B 
C 
D 
E 
F 
G 


(10  M) 

0.200  (radians) 
0.185 
0.157 
0.117 
0.061 
0.028 
0.012 


(30  M and  100  M) 
0.262  (radians) 
0.237 
0.184 
0.119 
0.056 
0.023 
0.009 


The  correlation  shown  in  Table  VIII-2  was  developed  by  SAI  based 
on  data  presented  in  the  Shoreline  Diffusion  Program  by  Smith 
and  Niemann  (32) . The  vertical  eddy  viscosity  coefficient,  Kv, 
can  be  specified  using  Equation  VIII-34  and  Tables  VIII-1  and 
VIII-2  if  the  vertical  height,  atmospheric  stability  category 
and  local  mean  velocity  are  known. 

Therefore,  SAI  assigns  a value  of  Ky,  the  vertical  eddy 
viscosity,  at  each  grid  point  of  their  numerical  solution  based 
on  the  vertical  height  and  local  velocity  calculated  at  that 


point  and  the  atmospheric  stability  category  which  is  assume 
to  characterize  that  location.  In  order  to  assign  an  atmos- 
pheric stability  category  at  a given  location  at  a given  time, 
SAI  uses  the  method  proposed  by  Smith  and  Howard  (33)  in  which 
the  atmospheric  stability  category  is  correlated  with  the 
vertical  temperature  gradient. 
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SAI  applies  the  classical  Reynolds  analogy  to  equate  the 
turbulent  transfer  coefficients  for  momentum,  mass  and  energy. 
The  components  of  the  vertical  eddy  viscosity  coefficients 
(see  Equation  VIII-14a)  are  also  assumed  equal. 

Hence 

kx3  = kx6  = KV  . (VIII-35) 

and  k2  = K'  = Ky 


Finally,  SAI  assumes  equality  of  turbulent  transfer 
coefficients  for  momentum,  mass,  and  energy  in  the  horizontal 
plane  and  assumes  the  x and  y components  of  these  coefficients 
equal,  hence 


and 


Kxl  = Kx2 
kx  = ky 
Kx  " Ky 


kx4  kx5  " kh 

kH 

KH 

kh 


The  horizontal  transfer  coefficients  (kH  * Kjj  * K^)  are  then 
estimated  from  the  ratios  of  horizontal  to  vertical  transfer 
coefficients  shown  in  Table  VIII-3.  Table  VIII-3  indicates 
enhancement  of  vertical  "diffusive"  power. of  the  atmosphere 
relative  to  horizontal  diffusive  power  when  the  atmosphere  is 
unstable.  The  ratios  in  Table  VIII-3  are  based  on  proprietary 
field  data  obtained  by  SAI. 


TABLE  VIII-3.  RATIO  OF  HORIZONTAL  AND  VERTICAL 

DIFFUSIVITIES  VERSUS  STABILITY  CLASS 

RATIO  PASQUILL  STABILITY  CLASS 

D E F 

KH  ” 

1-0  10  25 
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SPECIFICATION  OF  VAPOR  RELEASE  RATE 


SAI  assumed  that  the  liquid  pool  resulting  from  an 
instantaneous  release  of  LNG  onto  water  is  circular  and  increases 
in  size  as  described  by  the  following  equation: 


(VXII-37) 


where  r 
t 


g 

pw,pl 

h 


= pool  radius,  ft 
= time,  sec 

2 

= gravitational  acceleration,  ft/sec 

= densities  of  water  and  LNG,  respectively, 

any  consistent  units 


= pool  depth,  ft 


Equation  VIII-37  is  used  to  describe  the  growth  of  the  spill 
pool  until  a minimum  pool  thickness  is  reached  at  which  time 
the  pool  is  assumed  to  break  up.  The  minimum  pool  thickness 
is  determined  from  the  relation  proposed  by  Feldbauer  et  al. 
(3)  based  on  API  sponsored  Matagorda  Bay  test  data: 


7jmin 


0.0017 


D0-56 


(VIII-38) 


where  Zm^n  = minimum  pool  thickness,  ft 
D - pool  diameter,  ft 

Following  pool  breakup,  the  evaporation  rate  is  assumed  to 
decrease  according  to  the  following  relation  also  proposed  by 
Feldbauer  (3) : 
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exp  [-  0 .04  (t- 

L Pzmin 

= evaporation  rate  at  time  t,  lb/sec  |j 

= evaporation  rate  at  time  of  pool  ! I 

breakup,  Ib/sec  ij 

= LNG  density,  lb/ft3  ] 

’ J 

= time  of  pool  breakup,  sec 

Equations  VIII-37,  38,  39  are  used  to  calculate  the  pool  ’ 

evaporation  rate,  assuming  a constant  boiling  rate  per  unit 

2 

area  of  0.04  lb/ft  sec. 

Independent  calculations  by  the  author  of  the  vapor 
dispersion  following  a 25,000  M3  instantaneous  spill  were  not 
possible  due  to  the  proprietary  nature  of  SAI's  computer 
programs.  SAI  has  not  published  calculated  results 
for  a 25,000  M spill.  Therefore,  SAI's  results  for  a 
37,500  instantaneous  spill  are  discussed  here  for  comparison 
with  the  previous  estimates.  Table  VIII-4  shows  SAI's 
predictions,  based  on  Equations  VIII-37,  38,  39,  for  liquid 
pool  size  and  evaporation  rate  for  a 37,500  M^  instantaneous 
spill  onto  water.  Table  VIII-4  shows  total  vapor  production 
rate  as  a function  of  time.  In  SAI's  computer  simulation,  the 
evaporating  pool  is  represented  as  a variable  area  source  by 
simulation  of  LNG  vapor  addition  to  the  atmosphere  at  the 
appropriate  grid  points  indicated  in  Figure  VIII-1. 

The  downwind  distance  to  the  time  average  5%  vapor 
concentration  for  a 37,500  cubic  meter  spill  in  a 3 m/sec 
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MAX 


where  W 

W 

MAX 

P 


MAX 
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(6.7  mph)  wind  calculated  using  SAI's  model  (as  reported  by 
SAI  )is  1.2  miles. 


TABLE  VIXI-4.  EVAPORATION  RATE 
AND  LIQUID  POOL  RADIUS  PREDICTED  BY  SAI 
FOR  37 / 500  M3  INSTANTANEOUS  LNG  SPILL  ONTO  WATER 


VAPOR 

PRODUCTION  RATE 

Ft^/sec  at  7 OF, 

Time , sec  Pool  Radius,  ft  lb/sec  1 atm  


50 

620 

4.9 

X 

104 

1.2 

X 

106 

100 

869 

8.7 

X 

104 

2.1 

X 

106 

150 

1050 

14.0 

X 

104 

3.3 

X 

io6 

200 

1184 

17.8 

X 

104 

4.2 

X 

106 

250 

1184 

11 . 6 

X 

104 

2.8 

X 

106 

300 

1184 

6.8 

X 

104 

1.6 

X 

106 

350 

1184 

3.9 

X 

104 

9.3 

X 

105 

450 

1184 

13.4 

X 

103 

3.2 

X 

105 

520 

1184 

6.3 

X 

103 

1.5 

X 

105 

137 


A 


! 

1 

i 

IX  - ASSESSMENT  OF  LNG  VAPOR  DISPERSION  PREDICTABILITY  FOR  j 

i 

CATASTROPHIC  SPILLS  ONTO  WATER  j 

I 

j 

Published  predictions  of  LNG  vapor  cloud  formation  and  ] 

dispersion  following  a catastrophic  spill  of  LNG  on  water  can 
be  categorized  as  follows: 

1.  Predictions  which  utilize  classical  air  pollutant 
dispersion  models  originally  developed  to  describe  relatively 
near-field  dispersion  of  neutrally  buoyant  materials.  These 
models  are  based  on  the  general  observation  that  the  concen- 
tration profiles  downwind  of  a pollutant  source  are  reasonably 
accurately  represented  by  a Gaussian  or  normal  distribution. 

This  model  type  is  further  subdivided  to  describe  two 
different  dispersion  phenomena? 

a.  Dispersion  of  an  essentially  instantaneous 
release  of  a given  amount  of  pollutant  into  the  atmosphere, 
the  dispersion  being  associated  with  the  growth  of  this 
instantaneously  released  "puff",  or  cloud,  as  it  is  being 
translated  by  the  wind. 

b.  Dispersion  of  material  which  is  being  emitted 
at  a continuous,  steady  rate  forming  a "plume"  downwind  of 
the  emission  source. 

2.  Predictions  based  on  solution  of  the  combined  mass, 
momentum  and  energy  balance  equations.  The  classical  air 
pollutant  dispersion  equations  of  category  1 above  are  a 


138 


special  case  where  energy  effects  and  momentum  effects  are 
not  explicitly  considered.  In  cases  where  the  material 
added  to  the  atmosphere  has  a substantially  different  tempera- 
ture and  density  than  that  of  the  atmosphere  consideration  of 
energy  and  momentum  effects  can  be  important. 

Comparison  of  published  predictions  of  the  downwind 
travel  of  flammable  gas-air  mixtures  following  the  instantaneous 
release  of  25,000  cubic  meters  of  LNG  onto  water  identifies  the 
sensitivity  of  such  predictions  to  the  following  factors. 

a.  Characterization  of  atmospheric  stability 

b.  Allowances  for  area  source  effects 

c.  Specification  of  vapor  release  rate 

d.  Allowances  for  gravity  spread/air  entrainment  effects 

The  choice  of  atmospheric  stability  category  assumed 

applicable  to  the  accident  scenario  strongly  affects  the 
downwind  distances  predicted  using  models  based  on  the  classical 
pollutant  dispersion  equations.  The  use  of  stability 
characterizations  other  than  those  representing  "inversion" 
or  very  stable  conditions  for  "worst  case"  evaluation  is 
difficult  to  justify,  in  the  author's  opinion,  since  the 
latter  may  occur  frequently. 

Allowances  for  area  source  effects  incorporated  with 
the  classical  pollutant  dispersion  equation  models  rely  on 
specification  of  a point  "virtual"  source  (CHRIS,  Germeles 
and  Drake,  FPC) , or  line  source  representation  (Feldbauer) , for 
the  predictions  shown  herein.  Incorporation  of  these 
techniques  affects  the  predicted  distances  more  for  unstable 
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weather  conditions  than  for  stable  conditions,  since  the 
correction  for  the  initial  spreading  effect,  is  a smaller 
percentage  of  the  total  distance  in  the  latter  case. 

This  review  shows  that  much  of  the  variation  in 
predicted  downwind  distances  is  due  to  differences  in 
estimation  of  the  rai  e of  vapor  flow  into  the  atmosphere. 

For  example,  the  she ctest  distance  predicted  for  an 
instantaneous  25,000  cubic  meter  spill  is  0.75  miles  by 
the  FPC  staff.  Th  s short  distance  can  be  viewed  as  resulting 
primarily  from  the  low  estimated  rate  of  vapor  flow  into 
the  atmosphere  and  to  a lesser  but  still  important  degree 
from  the  use  of  neutral  weather  stability  dispersion  coeffi- 
cients. It  should  be  noted  that  the  FPC  model  predicts  an 

3 

evaporation  time  of  only  4.5  minutes  for  a 25,0Gu  m 
instantaneous  spill  which  corresponds  to  an  average  vapor 
production  rate  approximately  the  same  as  predicted  by  the 
ADL  - CHRIS  model.  However,  the  FPC  staff  assumes  that  this 
vapor  "piles  up"  above  the  liquid  pool  in  pure  form  and  only 
begins  to  enter  the  atmosphere  after  the  cloud  becomes 
neutrally  buoyant  (i.e.  when  its  density  reaches  that  of  the 
air) , during  which  time  it  spreads  as  a pure  cloud  to  a 
diameter  of  3,785  feet.  The  FPC  then  assumes  the  rate  of 
vapor  "release"  from  this  pure  cloud  to  the  atmosphere  is 
limited  by  the  rate  of  heat  transfer  from  the  surrounding 
air  to  the  cloud's  upper  surface.  This  assumption  results  in 
a vapor  release  rate  of  6,250  lb/sec,  which  indicates  the 
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cloud  would  release  vapor  from  its  top  surface  at  this  steady 
rate  for  almost  18  minutes.  There  appears  to  be  no  technical 
justification  for  this  description  of  the  vapor  flow  rate  and 
it  is  considered  by  the  author  to  be  unacceptably  low. 

Furthermore,  as  long  as  the  classical  air  pollutant  dispersion 
models  are  used,  there  does  not  appear  to  be  any  valid  reason 

\ 

why  the  worst  case  atmospheric  stability  conditions  should  not 
be  used  to  pred5  ;t  the  maximum  downwind  distance.  For  thee- 
reasons,  the  short  distance  predicted  by  the  FPC  staff  cannot  be 
accepted  based  on  their  technical  arguments. 

LNG  vapor,  when  it  is  initially  formed  at  the  boiling 
pool  surface,  is  at  a temperature  of  about  -260°F  and  the 
vapors  at  this  temperature  are  almost  1 1/2  times  as  heavy 
as  air.  When  large  quantities  of  this  dense  vapor  are 
rapidly  released  into  the  atmosphere  the  cloud  formed  should 
tend  to  remain  close  to  the  water  surface,  i.e.,  its 
vertical  dispersion  should  be  suppressed.  The  experimental 
spills  which  have  been  made  on  water  to  date  (1,  2,  3,  10) 
confirm  this  behavior.  Fay  (6),  Germeles  and  Drake  (4), 

Burgess  (1,  2),  Feldbauer  (3)  and  the  FPC  (7)  have 
all  attempted  to  modify  or  augment  the  classical  pollutant 
dispersion  models  to  account  for  this  effect.  However,  the 
methods  used  for  this  purpose  by  these  groups  are  not 
similar,  and  the  predicted  effect  on  dispersion  directly 
attributable  to  gravity  spread  action  varies  from  slight 
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(Fay,  Burgess,  FPC)  to  very  large  (Feldbauer) . The  results 
obtained  from  the  Germeles  and  Drake  model  are  sensitive 
to  the  numerical  values  of  the  parameters  which  relate  to 
the.  gravity  spread  phase  and  its  associated  air  entrainment. 

An  important  pattern  can  be  recognized  in  the  techniques 
surveyed  that  are  based  on  classical  air  pollution  dispersion 
models.  Where  gravity  spread  has  been  considered  along  with 
air  entrainment  by  the  advancing  gravity  spreading  cloud, 
results  show  that  inclusion  of  both  effects  can  markedly 
reduce  tne  prediction  of  downwind  travel  to  the  lower  flam- 
mable limit  for  very  large,  rapid  spills.  Variation  of 
the  air  entrainment  parameter  in  Germeles  and  Drake's  model 
by  a factor  of  5 results  in  prediction  of  the  average  cloud 
concentration  dropping  below  the  lower  flammable  limit 
during  the  gravity  spread  phase  of  a 25,000  M3  instantaneous 
spill.  Although  variation  of  the  air  entrainment  parameter 
by  a factor  of  5 upward  (and  10  downward)  may  not  represent 
a physically  realizable  range,  it  does  show  the  sensitivity 
of  the  resulting  prediction  to  the  numerical  quantification 
of  the  air  entrainment.  In  view  of  the  suggested  sensitivity 
of  the  result  to  the  degree  of  air  entrainment  by  the  spreading 
cloud,  the  importance  of  correctly  modeling  the  dispersion 
of  the  cloud  associated  with  gravity  spreading  of  the  dense 
cloud  is  apparent. 

If  gravity  spreading  induced  effects  are  not  considered 
to  be  important,  classical  models  suggest  that  substantial 
downwind  travel  of  the  cloud  will  occur  before  the  concentration 
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decreases  to  the  non-flammable  range.  Furthermore,  the 
predicted  downwind  distance  to  the  lower  flammable  limit 
following  very  large  spills  of  LNG  appears  to  depend  strongly 
on  the  degree  of  dispersion  attributed  to  the  initial  gravity 
spread  phase. 

In  the  author's  opinion  future  attention  should  be 
centered  on  the  development  and  verification  of  models  which 
include  some  explicit  procedure  for  describing  the  early 
development  of  the  cloud,  including  a method  for  quantifying 
the  air  entrainment  which  may  be  associated  with  gravity 
spread  induced  turbulence.  The  Germeles  and  Drake  model  and 
the  SAI  model  both  address  this  need;  Burgess'  model.  Fay's 
model,  the  FPC  model  and  the  CHRIS  model  do  not.  The  model 
suggested  by  Feldbauer  provides  for  mixing  of  air  and  vapor 
during  gravity  spread  by  assuming  the  cloud  depth  to  remain 
constant  during  the  spreading  process.  This  approach  appears 
to  be  based  on  the  experimental  observations  of  Feldbauer, 
and  the  validity  of  extension  to  very  large  spills  is  uncertain. 

It  is  also  the  opinion  of  the  author  that  any  model  to 
be  used  for  predicting  the  dispersion  of  vapor  from  very 
large  spills  should  take  into  account  energy  effects 
associated  with  mixing  LNG  vapor  and  air.  Furthermore, 
since  the  most  important  question  concerning  the 
validity  of  previously  used  models  concerns  the  degree  of 
dispersion  which  may  result  due  to  the  action  of  the  cloud 

itself  (i.e.  by  gravity  spreading  and  associated  air  entrain- 
ment) , future  attention  should  be  centered  on  development  of 


models  which  are  capable  of  accounting  for  simultaneous  effects  of 
energy  transfer,  gravity  induced  spreading,  and  turbulent 
diffusion.  The  Germeles  and  Drake  model  accounts  for  gravity 
spread  effects,  energy  effects  associated  with  mixing  LNG 
vapor  and  air,  and  air  entrainment,  utilizing  a lumped  para- 
meter approach  which  assumes  the  developing  cloud  to  be 
spatially  uniform  (but  changing  with  time)  during  the  gravity 
spread  phase.  The  Germeles  and  Drake  model  provides  a 
framework  for  inclusion  of  important  physical  effects,  even  if  in 
a simplified  form.  The  SAI  model  accounts  for  gravity  spread, 
energy  effects  associated  with  mixing  LNG  vapor  and  air,  and 
air  entrainment  by  solving  a less  simplified  form  of  the  mass, 
energy  and  momentum  balances.  In  this  regard,  the  SAI 
technique  provides  several  advantages  as  follows: 

1 . The  technique  allows  for  a representative  descrip- 
tion of  the  true  transient  nature  of  the  spill  phenomena. 

For  example,  the  rate  of  vapor  production  from  the  spill  can 
be  represented  in  a much  more  realistic  time  varying  form. 

2.  Inclusion  of  the  energy  balance  equations  allows 
description  of  the  temperature  development  of  the  cloud  in  a 
more  realistic  way.  In  the  SAI  method,  the  temperatures  and 
concentrations  in  the  cloud  are  considered  to  be  functions  of 
both  time  and  location,  whereas  even  the  most  sophisticated 
previous  models  (Germeles  and  Drake)  assume  the  cloud  tempera- 
ture and  concentration  during  the  initial  phase  of  development 
to  be  uniform  while  varying  with  time. 

3.  Phenomenological  relationships,  particularly  the 
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I 

. 

I 

coefficients  of  turbulent  diffusion,  can  be  specified  as  a ' 

function  of  both  time  and  position.  This  is  significant 

! 

since  the  turbulent  diffusion  properties  of  the  cloud  would  i 

i 

be  expected  to  vary  in  both  time  and  space  due  to  the  ] 

j 

progressive  mixing  of  the  cold  vapor  with  air.  The  simpler  j 

j 

classical  models  assume  implicity  that  the  turbulent  ’ 

diffusion  of  the  vapors  occurs  without  affecting  the  pre-existing 
turbulence  patterns  in  the  atmosphere. 

The  primary  reason  for  the  much  shorter  downwind 
distances  to  the  5%  concentration  level  predicted  by  SAI  for 
a catastrophic  spill  appears  to  be  the  predicted  highly 
turbulent  motion  associated  with  the  gravity  spread  phase. 

This  high  degree  of  predicted  turbulence  at  the  spreading 
cloud-air  interface  is  responsible  for  significant  air  entrain- 
ment by  the  cloud.  Since  the  predicted  turbulence  is  primarily 
induced  by  the  spreading  action  of  the  cloud,  this  provides 
an  explanation  for  why  the  turbulence  properties  assigned  to 
the  surrounding  atmosphere  at  the  time  of  the  spill  (i.e.  neutral 
vs.  stable)  do  not  markedly  affect  SAI's  predicted  results. 

The  results  of  the  model  predictions  indicate,  that  the  principal 
dispersion  of  the  vapor  to  the  point  where  the  concentration 
is  below  51  is  associated  with  effects  caused  by  the  cloud 
spread  itself,  rather  than  the  prevailing  atmospheric  conditions. 

It  is  interesting  to  note  that  the  gravity  spread  analyses 
proposed  by  Germeles  and  Drake  (4)  and  Feldbauer  (3)  lend 
support  to  this  idea. 
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However,  it  is  also  the  author’s  opinion  that  certain 
questions  must  be  answered  concerning  the  predictions  of  the 
SAI  model  before  the  results  cited  herein  can  be  confidently 
accepted.  The  specification  of  the  turbulent  diffusion 
coefficients,  e.g.,  the  "eddy  viscosity"  coefficient,  must  be 
more  carefully  evaluated.  As  has  been  described  in  Section 
VIII,  the  local  specification  of  these  transfer  coefficients 
is  a rather  complex  process  involving  several  assumptions. 

The  assumption  of  equality  of  coefficients  representing  mass, 
momentum,  and  energy  transfer  requires  careful  scrutiny.  It 
should  be  noted  that  Hanna  (29)  , whose  work  provides  the  basis 
for  SAI’s  estimation  of  turbulent  diffusion  coefficients,  did 
not  generally  support  this  assumption  for  determination  of  the 
coefficient  of  energy  transport.  Smith  and  Niemann  (32)  have 
raised  questions  about  the  general  validity  of  the  basic 
relationships  between  the  energy  spectrum  and  the  wind  speed 
and  direction  variation  assumed  by  Hanna  (20)  and  used  by  SAI. 
Furthermore,  it  is  difficult  to  assess  the  uncertainty  in  the 
ultimate  specification  of  the  coefficients  introduced  by  the 
relations  proposed  by  Taylor  (31)  for  turbulence  scale  length. 
For  example,  the  data  of  Taylor  et  al. , on  which  Equation 
VIII-33  is  based,  is  very  scattered.  Pasquill  (12)  has 
questioned  the  validity  of  such  precise  correlations  of  the 
turbulence  scale  length  with  energy  spectrum  parameters. 
Finally,  the  method  used  by  SAI  still  involves  the  requirement 
to  assign,  locally,  stability  categories  of  the  classical 
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type  (e.g.,  Pasquill  A-F)  to  the  developing  cloud.  The 
method  used  is  based  on  the  work  of  Smith  and  Howard  (33) 
and  considers  the  stability  category  to  be  a function  of 
the  local  temperature  gradient.  This  correlation  is  based 
on  measurements  of  atmospheric  turbulence  under  relatively 
stationary  (in  the  statistical  sense)  conditions  and  the 
assumption  that  the  same  correlation  applies  in  a methane 
rich  cloud  which  is  in  a highly  nonstationary  state  is  not 
obvious.  It  must  be  emphasized  that  the  advantages  which 
obtain  from  the  use  of  a complex  model  such  as  that  proposed 
by  SAT  can  be  easily  vitiated  by  the  incorporation  of  techniques 
for  descriptions  of  turbulence  which  are  not  easily  verified. 
Until  further  studies  validate  this  part  of  the  overall 
approach,  we  may  only  be  trading  uncertainty  in  the  classical 
models  for  a new,  but  no  less  important,  uncertainty  in  more 
complex  models . 

There  remains  the  problem  of  verification  of  the  numerical 
procedures  used  in  the  computer  solution  of  the  SAI  model. 

This  study  did  not  address  the  need  for  a thorough,  independent, 
evaluation  of  the  computer  program  to  verify  .the  numerical 
accuracy  and  stability  of  the  solution  technique. 

There  are  other  techniques  which  might  be  applied  to 
the  vapor  dispersion  problem.  The  obvious  one  which  might  be 
suggested  is  to  use  a turbulence  "closure"  model  of  higher 
order.  These  methods  are  proposed  when  the  assumption  of 
proportionality  between  the  mean  gradients  of  concentration. 
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velocity,  and  temperature  and  their  respective  turbulent 
transfers  is  not  considered  applicable.  A large  body  of 
literature  has  developed  in  this  area  (34,  35,  36)  but  this 
study  has  not  addressed  it  in  detail.  However,  until  the  basic 
questions  posed  above  concerning  the  SAI  model  are  answered, 
there  is  little  justification  for  pursuing  a more  sophisticated 
model.  The  adage  that  more  complexity  does  not  insure  more 
validity  applies  directly  to  this  problem. 

There  are  other  important  questions  related  to  the 
predictability  of  vapor  dispersion  from  catastrophic  spills 
of  LNG  on  water  which  have  not  been  addressed  in  this  report. 
All  of  the  predictions  of  downwind  distance  which  have  been 
surveyed  in  this  report  have  been  compared  at  the  time  average 
5%  concentration  level.  For  comparison,  the  downwind  distance 
to  lower  time  average  concentrations  has  been  shown  for  some 
of  the  models.  There  is  still  disagreement  as  to  the  magnitude 
of  the  peak-to-average  concentration  ratios  that  would 
characterize  a vapor  cloud  resulting  from  a catastrophic  spill 
and  this  affects  the  choice  of  time  average  concentration  which 
limits  the  flammable  region  of  the  cloud.  in  the  author's 
opinion,  this  uncertainty  does  not  affect  the  comparative 
assessment  of  the  models  discussed  in  this  report.  Unless  a 
model  can  be  developed  which  provides  accurate  time  average 
concentrations,  the  accurate  prediction  of  peak-to-average 
concentration  rates  effects  cannot  be  anticipated. 

An  additional  facet  of  LNG  vapor  cloud  dispersion  which 
is  important  to  the  assessment  of  potential  hazard  is  the  width 


of  the  flammable  zone.  Analysis  of  the  models  surveyed 
in  this  report  indicates  a marked  variation  in  the  width  of 
potentially  hazardous  zones.  Since  the  area  exposed  to 
potential  ignition  sources  and/or  the  burning  cloud  determine 
the  exposure  to  such  an  accident,  an  accurate  estimate  of  the 
shape  of  the  cloud  is  required.  However,  the  author  believes 
that  the  comparison  of  the  models  described  in  this  report  based 
only  on  predicted  downwind  extent  of  the  flammable  zone  is 
sufficient  to  justify  the  assessment  made  and  the  recommendations 
for  further  evaluation  which  are  offered. 

It  is  important  to  re-emphasize  that  this  report  is 
intended  to  deal  only  with  the  predictability  of  catastrophic 
LNG  spills  on  water.  The  conclusions  to  be  drawn  are  not 
necessarily  appropriate  for  the  consideration  of  the  predicta- 
bility of  vapor  dispersion  from  small  LNG  spills  on  water  or 
land.  For  small  spills,  the  allowance  for  heat  transfer 
effects  and  momentum  transfer  effects  in  the  prediction 
of  the  dispersion  appears  much  less  important.  Experimental 
evidence  from  LNG  spills  on  the  order  of  10  cubic  meters  and 
smaller  support  this  contention.  Figure  IX-1  is  a comparison 
of  downwind,  ground-level  concentrations  predicted  using  Burgess' 
model  and  the  SAI  model,  as  described  in  this  report,  with 
experimental  data  from  an  American  Gas  Association  experimental 
program  (33) . The  spill  described  was  a rapid  release  of 
14,000  gallons  of  LNG  on  land.  The  spill  was  confined  by  an 
80  feet  diameter,  1.5  feet  high  dike.  Maximum  vapor  production 
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AGA  MEASUREMENTS 
FROM  TEST  044 


Figure  IX-I.  COMPARISON  OF  GROUND  LEVEL  CONCENTRATIONS  FROM  AGA 

TEST  044  WITH  PREDICTIONS  BY  SAI  MODEL  AND  BURGESS'  MODEL 
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rates  were  measured  in  the  experiments  and  reported  as  0.09 

o 

cubic  meters  per  second  at  -260  F (0.72  inches  LNG  per  minute). 
The  atmospheric  stability  conditions  reported  for  the  experi- 
ment (A GA  044  in  the  test  series)  were  reported  as  Pasquill 
"C",  with  a wind  velocity  of  12  miles  per  hour.  The  vertical 
hash  marks  represent  the  range  of  concentrations  measured  at 
downwind  positions.  The  solid  predicted  curve  is  taken  from 
SAI's  published  risk  assessment  study  for  Oxnard,  California  (8) 
The  dashed  line  was  calculated  by  the  author  using  Burgess' 
model  with  Singer  and  Smith  stability  coefficients. 

Following  Burgess,  the  vertical  dispersion  coefficient  ct 

z 

was  equated  to  0.2 o-y  and  the  vapor  source  was  assumed 
concentrated  at  the  pool  center.  The  maximum  experimentally 
measured  vapor  production  rate  was  used  in  both  models. 

Note  that  the  downwind  distance  to  the  5%  level  is 
essentially  the  same  for  both  models.  This  is  in  contrast 
to  the  difference  in  downwind  distances  to  the  5%  concentra- 
tion predicted  using  these  models  for  a 25,000  cubic  meter 
spill  as  shown  in  Table  IV-1  of  this  report.  Two  things  are 
immediately  apparent  from  this  comparison.  First,  sufficient 
accuracy  may  be  obtainable  from  both  classical  dispersion 
models  and  the  SAI  model  for  the  prediction  of  LNG  vapor 
dispersion  from  small  spills  on  land  or  water.  Second, 
experimental  data  from  small  spills  validates  several  models 
for  prediction  of  that  type  of  phenomena,  while  shedding  no 
light  on  the  question  of  validity  of  the  models  for  predicting 
vapor  dispersion  from  very  large  spills. 
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X.  CONCLUSIONS 


1.  This  review  and  comparison  of  published  predictions 
of  the  downwind  travel  of  flammable  gas-air  mixtures 
following  the  instantaneous  release  of  25,000  cubic  meters 
of  LNG  onto  water  identifies  the  sensitivity  of  such 
predictions  to  the  following  factors. 

a.  Characterization  of  atmospheric  turbulence 
(stability) 

b.  Allowances  for  area-source  effects 

c.  Specification  of  vapor  release  rate 

d.  Allowances  for  gravity  spread/air  entrainment  effects 


2.  The  shortest  distance  to  the  time  average  5%  concentra- 
tion level  for  a 25,000  cubic  meter  instantaneous  spill 
predicted  by  the  models  reviewed  herein  is  0.75  miles.  This 
distance,  predicted  by  the  FPC  model,  results  primarily  from 
the  use  of  an  unrealistically  low  vapor  release  rate  and  the 
use  of  neutral  atmosphere  stability  characteristics.  The 
FPC  estimate,  in  the  author's  opinion,  is  not  justified. 


3.  At  the  other  extreme,  distances  of  the  order  of  tens  of 
miles  are  predicted  for  a 25,000  cubic  meter  instantaneous 
spill  under  stable  weather  conditions  using  continuous  plume 
models  (Burgess)  which  do  no/  account  for  any  heat  transfer 


cr  momentum  transfer  effects.  Such  estimates  are  not  justified 
in  this  author's  opinion. 

4.  Intermediate  distances  to  the  5%  concentration  level 
are  predicted  for  a 25,000  cubic  meter  spill  during  stable 
weather  conditions  by  Germeles  and  Drake  (11.5  miles). 

Fay  (17.4  miles)  and  the  CHRIS  model  (16.3  miles). 

The  difference  in  predicted  downwind  distances  obtained 
with  the  CHRIS  and  Germeles-Drake  models  can  be  attributed 
primarily  to  the  inclusion  of  gravity  spread/air  entrainment 
effects  in  the  Germeles  and  Drake  Model.  The  rough  agreement 
of  these  two  predictions  with  the  value  of  17.4  miles 
predicted  by  Fay  must  be  considered  fortuitous  since  the 
modeling  process  assumed  by  Fay  is  quite  different  from 
the  other  two.  Professor  Fay  now  believes  (42)  that  his 
model  should  be  used  with  different  assumptions  than 
originally  described  by  Lewis  and  Fay,  in  which  case  sub- 
stantially longer  distances  result.  In  the  author’s  opinion 
the  model  of  Germeles  and  Drake  provides  a more  plausible 
estimate  of  the  LNG  dispersion  process  following  a large 
rapid  spill  than  the  Fay  or  CHRIS  models,  since  the  model 
incorporates  a rational,  if  simplified,  description  of  an 
anticipated  gravity  spread  phase.  Further  effort  to 
improve  this  type  model  as  an  alternative  to  a more  complex 
numerical  procedure  has  merit,  particularly  for  routine 
usage  where  time  and  expense  constraints  are  important. 
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5.  The  estimate  using  Feldbauer's  model  of  5.2  miles  for 
the  downwind  distance  to  the  5%  concentration  level  following 
a 25,000  cubic  meter  spill  can  be  attributed  to  the  predicted 
dilution  and  corresponding  extreme  width  (—2  miles)  of  the 
cloud  at  the  end  of  the  gravity  spread  phase.  Feldbauer's 
allowance  for  air  entrainment  during  the  gravity  spread,  which 
involves  the  assumption  of  constant  cloud  depth,  is  based  on 
observations  of  small  spills  (10  M^)  and  the  extension  to  very 
large  spills  appears  uncertain.  Further,  representation  of 
the  cloud  at  the  end  of  the  gravity  spread  phase  as  a series 
of  dispersed  point  sources  on  a line  perpendicular  to  the 
direction  of  cloud  travel  is  not  realistic  in  view  of  the 
resulting  prediction  of  shorter  distances  with  increasing  atmos- 
pheric stability. 

6.  The  primary  reason  for  the  even  shorter  downwind  distance 
(~1  mile)  to  the  5%  concentration  level  predicted  by  SAI  for 
an  even  larger  (37,500  cubic  meters)  spill  appears  to  be  the 
predicted  highly  turbulent  motion  and  associated  air  entrain- 
ment induced  during  the  gravity  spread  phase  of  the  cloud. 

7.  In  the  author's  opinion,  the  predicted  maximum  distances 

of  about  5 miles  by  Feldbauer  and  about  1 mile  by  SAI  for  flam- 
mable cloud  travel  following  instantaneous  release  of  25,000 
cubic  meters  of  LNG  onto  water  cannot  be  rationalized  on  the 
basis  of  any  argument  thus  far  advanced  except  that  of  gravity 
spread/air  entrainment  effects,  and  experimental  verification 
of  these  effects  has  not  been  adequately  demonstrated. 

8.  It  was  not  possible  within  the  limits  imposed  by  this 
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review  to  evaluate  the  accuracy  of  the  predictions  published 
by  SAI.  Rather,  the  author  has  reviewed  the  methodology 
described  by  SAI  and  believes  that  such  techniques  hold 
the  most  promise  for  accurate  prediction  of  vapor  disper- 
sion from  catastrophic  spills  on  water.  A program  designed 
to  evaluate  the  accuracy  of  the  SAI  model  or  other  models 
of  similar  generality  should  now  be  considered  high  priority. 
The  Recommendations  section  of  this  report  addresses  this 
need . 
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XI.  RECOMMENDATIONS 

1.  The  Science  Applications,  Inc.  (SAI)  LNG  vapor  dispersion 
model  should  be  more  thoroughly  evaluated.  This  will  require 
the  cooperation  of  SAI  due  to  the  proprietary  nature  of  their 
computer  programs  which  are  required  for  solution  of  the  model 
equations.  Further  evaluation  of  the  SAI  models  or  other 
similar  models  based  on  simultaneous  solution  of  the  mass, 
momentum  and  energy  balance  equations  which  may  become 
available,  should  address  the  following  requirements: 
a.  The  methodology  of  describing  turbulent  mass, 
momentum  and  energy  transfer  should  be  critically 
evaluated.  A literature  search  should  be  conducted 
to  identify  and  evaluate  experimental  data  supporting 
the  assumption  of  first-order  (eddy  diffusivity, 
thermal  conductivity  and  viscosity)  turbulent  transfer 
phenomenological  relationships  for  describing  turbulent 
transfer  in  the  lower  atmosphere. 

b An  error  analysis  should  be  done  to  provide  some 
means  for  estimating  the  confidence  level  in  the  technique 
used  to  assign  numerical  values  to  the  turbulent  transfer 
coefficients . 

c.  Sufficient  calculations  should  be  made  with  the  model 
to  determine  the  sensitivity  of  the  results  predicted  by 
the  model  to  uncertainties  in  the  transfer  coefficients 
identified  in  b.  above. 

d.  An  analysis  should  also  be  made  of  the  liquid  spread, 
vapor  generation,  and  heat  transfer  models  used  in  the 
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specification  of  the  boundary  conditions  to  determine 
the  sensitivity  of  the  model  predictions, 
e.  The  numerical  stability  and  accuracy  of  the 
algorithm  used  for  computer  solution  of  the  equations 
should  be  critically  evaluated. 

2.  A series  of  computations  should  be  made,  using  the  SAI 
model,  of  the  downwind  distance  to  the  time  average  5% 
concentration  level  for  "instantaneous"  LNG  spills  as  a 
function  of  spill  size.  The  range  of  spill  sizes  should  be 
from  10  cubic  meters  to  25,000  cubic  meters  with  sufficient 
points  between  to  adequately  characterize  the  predicted 
relationship  between  flammable  cloud  travel  and  spill  size. 

3.  The  result  of  2 above  should  be  compared  with  a 
similarly  prepared  relationship  between  flammable  cloud 
travel  and  spill  size  predicted  using  the  Germeles  and 
Drake  model  described  in  this  report.  It  is  anticipated 
that  the  results  will  be  in  substantial  agreement  for  very 
small  spill  sizes  but  the  comparison  should  indicate  the 
smallest  spill  sizes  for  which  significant  differences  appear 
in  predicted  downwind  distance.  Such  e.  comparison  should 
also  provide  guidance  for  determining  a lower  bound  on  the 
size  of  experimental  spills  which  may  be  required  to  assess 
large  spill  behavior. 


4.  In  ancitipation  of  experimental  skills  which  may  be 
required  to  provide  confidence  in  predictions  of  large  spill 
behavior,  an  evaluation  should  be  made  of  the  experimental 
data  requirements  associated  with  verification  of  model 
predictions . 

5.  Additional  experimental  spills  should  be  performed  only 
after  completion  of  the  program  outlined  above,  and  such 
spills  should  be  performed  for  the  purpose  of  model  evaluation. 
Large  "demonstration  spills"  have  been  suggested  recently, 
largely  as  a result  of  the  variation  in  predictions  which 

has  been  the  subject  of  this  report.  It  is  the  opinion  of 
this  author  that  validation  of  models  should  still  be  the 
primary  goal  of  further  test  prograins;  "demonstration"  of 
the  effects  of  large  spills  without  heavy  reliance  on  models 
should  be  avoided. 
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APPENDIX  I 

DISPERSION  COEFFICIENT  DATA  FOR 
USE  IN  CLASSICAL  AIR  POLLUTANT 
DISPERSION  EQUATIONS 

Dispersion  coefficient  data  are  of  two  types: 

1.  Data  representing  the  cloud  width  (or  specified 
fraction  thereof}  as  a function  of  distance  traveled 
by  an  instantaneous  release  of  material  from  a point 
source.  Data  of  this  type  is  relatively  limited.  A 
survey  of  data  of  this  type  has  been  made  by  Slade  (13) 
from  which  the  suggested  correlations  shown  in  Tables 
A-I-l  were  proposed. 


TABLE  A-I-l 

SUGGESTED  ESTIMATES  FOR  oxrT  AND  c , SLADE  (13) 

yj.  zl 


Parameter 

Conditions 

x=100  meters 

x=4000  meters 

Approximate 

Correlation 

c , meters 

yi 

Unstable 

10 

300 

0.14 

x0.92 

Neutral 

4 

120 

0.06 

X°  *92 

Very  Stable 

1.3 

35 

0.02 

jO  .89 

czl,  meters 

Unstable 

15 

220 

0.53 

x0.73 

Neutral 

3.8 

50 

0.15 

x0-70 

Very  Stable 

0.75 

7 

0.05 

x0.61 

It  should  be  noted  almost  no  data  were  reported  by  Slade  for 

distances  beyond  about  4000  meters,  and  the  approximate  correla- 
tions for  dispersion  coefficients  as  a function  of  distance  were 
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in  the  .range  100  to  4000  meters. 

2.  Dispersion  coefficient  data  representing  the  standard 
deviation  of  the  horizontal  concentration  distribution,^, 
and  the  standard  deviation  of  the  vertical  concentration 
distribution,  az,  as  functions  of  travel  distance  from  a 
steady,  continuously  emitting  point  source.  The  horizon- 
tal and  vertical  coefficients,  „ , and  a respectively, 

y z 

used  by  all  investigators  have  been  obtained  from  two 
primary  sources. 

DISPERSION  COEFFICIENTS  PUBLISHED  BY  PASQUILL  <12 ) 

Pasquill  and  others  published,  around  1960,  estimation 

methods  for  and  a which  were  based  on  measurements  of 
Y z 

wind-direction  fluctuation.  Due  to  the  need  for  estimates  of 
dispersion  coefficients  when  wind  fluctuation  measurements  are 
not  available,  Pasquill  suggested  values  for  0y  and  <*z  based 
on  the  degree  of  atmospheric  stability.  He  further  suggested 
that  stability  be  estimated  from  wind  speed  and  insolation. 

The  correlations  proposed  by  Pasquill,  along  with  the  guidelines 
for  estimating  atmospheric  stability,  are  shown  in  Figures 
A-I-l  and  A-I-2 . 

DISPERSION  COEFFICIENTS  PUBLISHED  BY  SINGER  AND  SMITH  (21) 

Singer  and  Smith  published  estimation  methods  for  a y and  <jz 
derived  from  measurements  of  dispersion  of  oil  fog,  radioactive 
isotopes  and  uranine  dye  at  the  Brookhaven  National  Laboratory. 


166 


TRAVEL  DISTANCE  (m) 

RELATION  OF  PASQUILL  TURBULENCE  TYPES  TO  WEATHER  CONDITIONS 

A — Extremely  unstable  conditions  D — Neutral  conditions* 

B — Moderately  unstable  conditions  E — Slightly  stable  conditions 

C — Slightly  unstable  conditions  F — Moderately  stable  conditions 


Nighttime  conditions 


Surface  wind 
speed,  m/sec 

Daytime  insolation 
Strong  Moderate  Slight 

Thin  overcast 
or  > 4/8  cloudi- 
ness t 

< 3/8  cloudi- 
ness 

<2 

A 

A-B 

B 

2 

A-B 

B 

C 

E 

F 

4 

B 

B-C 

C 

D 

E 

6 

C 

C-D 

D 

D 

D 

>6 

C 

D 

D 

D 

, D 

* Applicable  to  heavy  overcast,  day  or  night 
t The  degree  of  cloudiness  is  defined  as  that  fraction  of  the  sky 
above  the  local  apparent  horizon  that  is  covered  by  clouds. 

Figure  A 1-1.  HORIZONTAL  DISPERSION  (COEFFICIENTS  BY  PASQUILL 
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Nighttime  conditions 

Surface  wind 

Daytime  insolation 

Thin  overcast 
or  > 4/8  cloudi- 
ness t 

^ 3/g  cloudi- 
. ness 

speed,  m/sec 

Strong 

Moderate 

Slight 

C 2 

A 

A-B 

B 

J 

2 

A-B 

B 

C 

E 

F 

4 

B 

B-C 

C 

D 

E 

6 

C 

CD 

D 

D 

D 

> 6 

C 

D 

D 

D 

D 

* Applicable  to  heavy  overcast  day  or  night 

t The  degree  of  cloudiness  is  defined  as  that  fraction  of  the  sky  above  the  local  apparent 
horizon  that  is  covered  by  clouds. 


Figure  A-l-2.  VERTICAL  DISPERSION  COEFFICIENTS  BY  PASQUILL 

I 
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The  most  important  sources  of  data  were  oil  fog  release 
experiments,  which  involved  emission  of  oil  fog  droplets  from 
a single  source  at  an  elevation  of  108  meters.  In  addition 
measurements  of  radioactive  isotope  (A^)  emission  from  the 
Brookhaven  Graphite  Research  Reactor  were  used  to  make 
"qualitative  or  at  best  crude  quantitative"  estimates  of  plume 
position  and  dimension.  The  source  height  in  this  case  was 
also  108  meters.  These  experiments  were  apparently  the  basis 
for  estimation  of  horizontal  and  vertical  dispersion  coefficients 
at  great  distances  from  the  release  point,  and  all  information 
at  distances  50  km  or  more  was  obtained  from  the  isotope  measure- 
ments. Uranine  dye  releases  from  a height  of  2 meters  also 
provided  a small  amount  of  data.  In  all  cases,  concentration 
data  were  mean  values  obtained  over  periods  ranging  from  30  to 
90  minutes.  The  atmospheric  stability  was  taken  into  account 
by  defining  5 "gustiness  classifications"  based  on  horizontal 
wind  direction  fluctuations  measured  at  the  release  site  with 
a Bendix  Friez  Aerovane  located  350  ft  (107  meters)  above  ground. 

The  definition  of  these  "gustiness  classifications"  is  shown 
be  low . 

Gustiness  Classification  Horizontal  Wind  Direction  Fluctuation 

A Fluctuations  of  wind  direction  > 90° 

B2  Fluctuations  ranging  from  40°  to  90° 

Fluctuations  ranging  from  14°  to  45° 

C Fluctuations  ranging  from  0 tc  15° 

D Essentually  no  fluctuation,  short 

term  fluctuation  do  not  exceed  15° 
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The  correlation  for  the  horizontal  dispersion  coefficient, 
o , developed  from  the  data  is  snown  in  Figure  A-I-3.  No 
correlation  was  developed  for  Type  A Gustiness  Classification 
since  the  condition  is  characterized  by  the  absence  of  organized 
horizontal  wind  flow  and  is  describabie  only  in  qualitative 
terms . 

The  vertical  structure  of  the  plumes  from  test  releases 
was  not  measured  directly;  the  vertical  dispersion  coefficients 
were  calculated  from  Equation  V-12  rearranged: 

Q ' 

a — 

z 

TT  o U C _ n 
y x,y=0,z=0 

The  correlation  for  the  vertical  dispersion  coefficient,  o , 

4M 

proposed  by  Singer  and  Smith  are  shown  in  Figure  A-I-4.  In 
contrast  to  the  correlation  proposed  by  Pasquill,  et.  al,  the 
Brc  . .'on  vertical  dispersion  coefficient  vs.  distance  is 
ted  as  being  of  the  form  ci=ax^,  similrr  to  the 
e J on  for  the  horizontal  dispersion  coefficient . Smith 
am1  Pas^- ill  presented  plots  of  typical  field  concentrations 
against  distance  which  show  reasonable  agreement  with  predictions 
using  their  dispersion  coefficient  correlations  out  to  about 
6000  meters.  Singer  and  Smith  emphasized  the  lack  of  precision 
in  the  definition  and  specification  of  the  vertical  dispersion 
coefficient,  which  are  tied  by  the  method  of  determination  using 
Equation  V-12  to  the  assumption  of  constant  wind  speed,  as  well 
as  the  assumed  correctness  of  the  model. 
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For  comparison  the  correlations  for  a and  proposed 
by  Pasquill  and  Singer  and  Smith  are  plotted  together  on  Figures 
A-I-5  and  A-I-6. 


172 


HORIZONTAL  DISPERSION  COEFFICIENT,  <’v,  METERS 


VERTICAL  DISPERSION  COEFFICIENT,  oZi  METERS 
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DEPARTMENT  OF  TRANSPORTATION 
UNITED  STATES  COAST  GUARD 


MAILING  ADDRESS: 

U-S.  COASTGUARD  (G“*MHM/83) 
WASHINGTON.  D C SOS90 

pmone;  (202) -426-2306 


* 10330/4-2/1 
4 February  1977 


To:  Distribution* 

Gentlemen: 

I am  enclosing  a copy  of  my  draft  report,  "Predictability 
of  LNG  Vapor  Dispersion  from  Catastrophic  Spills  onto 
Water;  An  Assessment,*  prepared  for  the  Cargo  and 
Hazardous  Materials  Division,  Office  of  Merchant  Marine 
Safety,  U.  S.  Coasc  Guard.  This  report  represents  my 
understanding  and  assessment  of  frequently  cited  predictions 
of  hazardous  vapor  cloud  travel  which  might  occur  in  the 
event  of  a catastrophic  accident  involving  a marine  LNG 
carrier.  My  description  of  the  technique  used  to  make  these 
predictions,  and  the  calculations  based  on  those  techniques, 
are  based  on  reports  prepared  by  the  investigating  groups 
which  developed  the  modeling  techniques.  In  some  instances, 
these  groups  have,  at  my  request,  provided  assistance  in 
this  effort.  Such  assistance  involved  discussions  to  clari- 
fy questions  which  I had  based  on  my  review  of  the  published 
reports  cited  in  the  report,  as  well  as  provision  of  computer 
codes  allowing  me  to  make  predictions  of  my  own  utilizing 
each  of  the  models.  However,  the  description  of  the  models 
and  the  associated  predictions  were  prepared  by  me.  I have 
purposely  not  included  the  Conclusions  and  Recommendations 
sections  and  the  Summary  (which  includes  same) . I consider 
these  sections  tentative  until  such  time  as  I have  received 
your  comments  on  the  accuracy  of  my  technical  review  of  this 
problem. 

It  is  my  intention  to  recommend  the  release,  by  the  U.  S. 

Coast  Guard,  of  the  completed  report  to  all  interested  parties. 

I hope  that  it  will  be  helpful  in  answering  some  of  the 
questions  which  prevail  in  the  area  of  safety  management  in 
LNG  transportation. 

1 respectfully  request  your  review,  as  a representative  of 
the  investigating  groups  whose  work  I have  discussed,  of  the 
technical  and  interpretive  accuracy  of  my  description  of  your 
model  and  the  associated  predictions.  It  is  my  intention  to 
make  your  comments,  and  any  revisions  or  rebuttals  which  may 
be  indicated,  a part  of  the  final  report. 

In  the  interest  of  releasing  the  final  report  as  soon  as 
• possible,  please  send  me  your  written  comments,  in  form  suitable 
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for  subsequent  inclusion,  by  25  February  1977.  Please  feel 
free  to  call  me  if  I can  clarify  any  point  in  the  report  or 
its  intended  development  to  final  form. 


Sincerely, 


A*<a4 


JERRY  H 


HAVENS 
Technical  Advisor 
Cargo  and  Hazardous 
Materials  Division 


Enel:  (1)  Draft  Report,  "Predictability  of  LNG  Dispersion 

from  Catastrophic  Spills  bnto  Water:  An  Assess 

ment" 


^Distribution: 

Dr.  David  S.  Burgess 

Pittsburgh  Mining  Safety  and  Research  Center 
U.  S.  Bureau  of  Mines 
4800  Forbes  Avenue 
Pittsburgh,  PA  15213 

Dr.  Walter  May 

Exxon  Nuclear  Co.,  Inc. 

Plains  Road 

Balls ton  Spa,  New  York  12020 

Mr . Don  Oakley 
Distrigas  Corporation 
125  High  Street 
Boston,  MA  02110 

Professor  James  Fay 
Room  3-246 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

Dr.  Theodore  Needles 
Federal  Power  Commission 
825  North  Capital  Street,  N.E. 

Washington,  D.C.  20426 

Dr.  Walter  England 
Science  Applications,  Inc. 

1200  Prospect  Street 
P.O.  Box  2351 
La  Jolla,  CA  92037 


United  States  Department  of  the  Interior 

BUREAU  OF  MINES 

4800  FORBES  AVENUE 
PITTSBURGH.  PENNSYLVANIA  1521S 

Pittsburgh  Mining  and  Safety  Research  Center 


February  25,  1977 


Dr.  Jerry  Havens 
Technical  Advisor 

Cargo  and  Hazardous  Materials  Division 
U.  S.  Coast  Guard  (G-MHM/83) 

Washington,  DC  20590 

Dear  Jerry: 

I have  read  your  manuscript  quite  carefully  and  find  nothing 
to  which  I can  object.  You  have  performed  a useful  service 
for  the  many  people  who  still  ask  questions  about  atmospheric 
dispersion. 


Sincerely  yours, 

David  Burgess 
Research  Supervisor 
Fires  and  Explosions 
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E?£ON  NUCLEAR  COMPANY,  Irvo. 

777  ■ 106th  Avenue  N E.,  Bellevue,  Washington  98004,  Telephone  (206)  455-5130 

Malta  Enrichment  Program,  Plains  Road,  Ba'ston  Spa,  New  York,  NY.  12020  .Telephone  (518  ) 899-2947 


March  7,  1977 


Mr.  Jerry  Havens  1 

Technical  Advisor  \ 

Cargo  and  Hazardous  Materials  Division  i 

Department  of  Transportation 
U.S.  Coast  Guard 
Washington,  D.C.  20590 

Dear  Mr.  Havens: 


I appreciate  The  opportunity  of  commenting  on  your  Draft  Report 
on  I.NG  Vapor  Dispersion.  My  comments  are  attached.  j 


1 feel  that  the  review  and  comparison  you  have  carried  out  is 
a very  useful  study. 

I indicated  during  our  telephone  conversation  that  I *iod 
written  a paper  on  the  "unmixed"  character  of  an  LNG  vapor/air 
plume.  A copy  is  attached.  I also  attach  a handwritten  deriva- 
tion of  the  Pasquill  equation,  as  promised  in  our  phone  conver- 
sation. J hope  you  can  read  it. 
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Very  truly  yours. 


c~r- 


W.G.  May 
Senior  Sc 
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Advisor 
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COMMENTS:  "Predictability  of  LNG  Vapor  Dispersion  from 

Catastrophic  Spills  Onto  Water:  An  Assessment". 


First,  it  is  interesting  to  see  a direct  comparison  of  the 
models  that  have  been  suggested  for  LNG-di spers ion , particularly 
the  comparisons  of  Table  I V - 1 . .The  number  that  is  quoted  for 
our  analysis,  5.2  miles,  is  a little  shorter  than  we  would  have 
anticipated  (7.6  miles).  We  never  calculated  results  for  such 
a large  spill  however,  25000  m3;  our  largest  (calculation)  was 
a 4000  m3  spill.  ’ believe  that  Mr.  Havens  has  followed  our 
procedure  correct, y,  so  that  his  number  of  5.2  miles  is  to  be 
preferred  over  our  simple  extrapolation. 

I agree  with  Dave  Burgess  that  the  lack  of  historical  perspec- 
tive in  the  Report  is  unfortunate.  I am  proud  of  the  work  that 
we  did  and  feel  strongly  that  some  of  the  later  analyses  borrowed 
heavily  from  it.  But  in  the  absence  of  any  historical  comment, 
the  old  analyses  have  to  suffer, 

I have  several  comments  concerning  the  presentation  of  our 
analysis,  where  either  I don't  agree,  or  I feel  that  useful 
comparisons  with  others  could  have  been  drawn, 

1.  The  Report  assigns  us  to  the  category  of  those  using 
a "steady-state''  model,  and  states  that  we  have 
used  our  measurements  of  "maximum  vapor  flow  rate  - as 
an  estimate  of  Q in  Eq.  VII-5".  I believe  that  this 
doesn't  do  justice  to  what  we  did,  and  certainly  carries 
an  implication  that  I don't  agree  with. 

First,  it  is  important  to  understand  Eq.  VII-5,  the 
“Pasquill"  equation.  It  is  easily  derived,  starting 
from  an  assumption  that  concentrations  follow  a Gaussian 
distribution.  The  equation  is  simply  a. material  balance, 
which  relates  the  amount  flowing  at  any  instant  to  the 
plume  velocity  and  the  concentration  level.  [The  amount 
flowing  is  calculated  as  the  amount  crossing. a plane  at 
right  angles  to  the  direction  of  the  wind.  Axial  diffusion 
is  neglected  - an  important  consideration  discussed  later]. 

The  important  point  is  that  the  "Q"  in  the  Fasquill  equation 
i s the  vapor  flow  rate.  The  vapor  flow  rate  is  not  an 
"estimate"  for  Q;  it  i_s  Q.  We  didn't  estimate  "Q",  we 
measured  it.  People  who  equate  the  evaporation  rate  to  Q 
are  making  an  "estimate". 

For  steady  state  spills,  of  course,  the  value  of  "Q"  (vapor 
flow  rate)  is  the  same  as  evaporation  rate.  But  for  instan- 
taneous spills  it  is  not;  measurements  showed  this,  and 
elementary  thinking  suggests  it  is  not.  [What  value  do  you 
assign  to  Vapor  Flow  Rate  as  wind  speed  approaches  zero?]. 


2.  I think  it  worthwhile  to  compare  our  use  of  the  Pasquill 
equation  with  the  "Puff"  model  (Eq.  V - 1 0 ) . I do  not 
believe  that  the  major  difference  is  "unsteady  state"  vs. 
"steady-state".  Our  analysis  was  certainly  not  a "steady- 
state  analysis.  Our  values  for  "Q"  were  transient  (i.e. 
never  reached  steady-state),  and  were  in  general  a long 
way  from  the  evaporation  rate  (i.e.,  the  value  that  would 
be  used  in  a steady-state  analysis). 

The  important  distinction,  it  seems  to  me,  is  the  way 
that  axial  dispersion  is  handled.  The  "Puff"  model  incor- 
porates an  axial  dispersion  factor;  that  is,  the  plume 
lengthens  out  as  it  goes  downwind  (while  spreading  as  well). 

Ne  did  not  use  an  axial  dispersion  factor;  our  plume  was 
pictured  as  keeping  the  same  length  that  it  had  at  the 
point  where  we  measured  it  (i.e-t  at  our  line  of  sensors). 

I consider  this  a weakness  in  our  model,  an  advantage  for 
the  "Puff". 

There  are'pros'  as  well  as  'cons',  however.  There  is 
difficulty  getting  the  Puff  model  started  in  a sensible  way. 
Undoubtedly,  the  plume  from  a large  instantaneous  spill  will 
be  stretched  out  alot  immediately  downwind  of  the  spill  point. 
The  wind  simply  cannot  drag  away  the  gas  as  fast  as  it  is 
evolved,  and  a large  cloud  accumulates;  particularly  at  low 
wind  speed.  An  instantaneous  spill  would  never  start  off 
downwind  as  a round  Puff.  Our  use  of  a measured  "Q"  has  an 
advantage  in  this  respect;  the  measured  "Q"  Is  the  result- 
ant of  alot  of  complex  interactions:  the  wind  attacks  the 
accumulation  at  the  spill  point;the  spill  is  not  quite  in- 
stantaneous; some  stretching  of  the  plume  has  occurred  in 
flowing  from  the  spill  point  to  the  point  where  the  measure- 
ment was  made. 

I believe  that  most  of  the  stretching  out  of  the  plume  is  a 
consequence  of  the  initial  conditions  (the  accumulation  of 
vapor  over  the  spill  point)  - particularly  at  low  wind 
speed.  I don't  think  the  Puff  Models  handle, this  very  well, 
while  some  evaluation  of  It  appears  automatically  in  our 
measurement  and  use  of  vapor  flow  rate  Q. 

I hadn't  seen  the  Germeles  and  Drake,  and  the  Fay,  Puff 
Models  before  the  Coast  Guard  write-up.  I have  tried  to 
compare  them  with  ours  (particularly  the  Germeles  and 
Crake  model),  and  some  comments  are  given  below.  But 
first,  it  seems  to  me  helpful  to  give  a little  historical 
persoecti ve. 

There  were  two  general  approaches  for  plume  analysis  at 
the  time  we  did  our  work.  One  attributed  the  plume  spread 
to  gravity  effects.  The  assumption  was  generally  made  that 
the  density  of  LNG  vapor/air  mixtures  was  about  linear  with 
composition  (for  adiabatic  conditions).  The  other  approach 
simply  used  the  standard  dispersion  due  to  the  weather. 


We  and  the  Bureau  of  Mines  (at  about  the  same  time),  calcu- 
lated the  effect  of  air  humidity  on  the  mixture  density; 
the  profound  effect  that  was  found  made  the  first  assumption 
(above)  untenable.  Simply  assigning  the  observed  results  to 
weather  was  also  untenable,  however;  the  plumes  were  much 
too  low  and  wide.  This  led  us  then  to  the  analysis  which 
includes  an  effect  of  both;  the  gravity  effect  controls 
initially,  but  as  the  plume  is  diluted  and  its  density 
approaches  that  of  air,  the  final  mixing  is  assigned  to  the 
weather . 

Apparently  the  Germeles  and  Drake,  (as  well  as  Fay?)  models 
follow  this  same  plan.  There  are  differences  in  details  - 
e.g.,  the  criterion  used  for  switching  from  "gravity"  spread 
to"weather"  spread,  and  others  - but  the  general  approach  is 
the  same. 

Some  of  the  similarities  and  differences  between  our  analysis 
and  the  later  Germeles  and  Drake  analysis  are  outlined  below: 

a.  The  gravity  spread  relationships  used  are  essentially 
the  same  (they  differ  by  a constant  coefficient). 

We  started  with  a relationship  (Fannelop  and  Waldman): 
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We  differentiated  with  respect  to  x,  to  get  the  change 
of  width  with  distance,  while  allowing  for  changing 
conditions  in  the  plume  (i  e. , the  analysis  keeps  track 
of  plume  temp,  density,  composition,  dimensions,  just 
as  the  Germeles  and  Drake  model  does). 
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If  instead,  we  had  differentiated  with  respect  to  time, 
while  making  tne  assumption  that  velocity  is  constant  so 
that  distance  and  time  are  related  (x  = ut),  we  would 
have  obtained  the  Germeles  and  Drake  gravity  spread 
(Eq.  VI-11) 
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The  only  difference  between  our  analysis  and  that  of 
Germeles  and  Drake  is  the  constant  coefficient;  [we 
would  get  18.5  for  the  G&D  spread  equation,  compared 
with  their  29.11].  I believe  that  there  are  questions 
about  the  Fannelop  and  Waldman  analysis  and  its  applica- 
tion to  this  work,  so  that  the  "right"  coefficient  is 
uncertain  [see  the  article  by  Hoult  in  Rev.  of  FI.  Mech.). 
[Incidentally,  the  Coast  Gurrd  Report  gives  different 
equations  for  the  G&D  spread  rate,  pp.  49  and  51..  I 
assume  the  one  on  p.  51  is  correct]. 

b.  In  our  analysis,  we  maintain  the  plume  height  constant 
during  the  gravity  spread.  This  has  the  effect  of  speci- 
fying, indirectly,  a mixing  coefficient.  The  procedure 
may  appear  arbitrary  but  was  based  on  experimental  obser- 
vation; our  plumes  appeared  to  rise  very  little  during 
their  downwind  travel.  The  Germeles-Drake  analysis 
introduces  a mixing  coefficient.  I note  however,  in  the 
example  given  in  the  Coast  Guard  report,  that  the  plume 
height  for  a 25000  m3  spill  increased  only  a very  small 
amount  during  the  gravity-spread  portion  of  its  travel; 
the  increase  was  from  13  m to  18m,  over  a downwind  distance 
that  was  presumably  several  miles.  I consider  that  to  be 
in  very  close  agreement  with  our  observation  of  constant 
plume  height.  The  mixing  rates  of  the  two  studies  must 
be  very  close. 

3.  We  assign  a downwind  speed  to  the  plume  which  varies  with 
concentration;  my  impression  - perhaps  incorrect  - is  that 
the  Germeles-Drake  analysis  assigns  constant  (wind)  speed. 

The  Coast  Guard  report  unfortunately  does  not  give  our 
rationale  for  varying  plume  speed:  we  calculated  conserva- 
tion of  momentum,  assuming  no  pressure  effects.  At  high 
concentration,  near  the  source,  the  velocity  is  therefore 
low,  but  approaches  wind  speed  at  large  dilution. 

The  experimental  data  confirm  this  type  of  effect.  The 
average  plume  speeds  measured  have  always  been  lower  than 
the  wind  speed,  by  substantial  factors,  e.g.,  3.  I consiaer 
our  approximation  more  acceptable  than  an  assignment  of  a 
constant  speed,  equal  to  the  wind. 

4.  We  have  apparently  used  somewhat  different  coefficients  than 
Germeles  and  Drake  for  our  analysis  of  the  effect  of  weather. 
The  important  coefficient  (at  least  in  cur  analysis),  is  the 
vertical  coefficient,  az.  The  plumes  are  already  so  wide  at 
the  end  of  the  gravity  spread  that  further  spreading  at  the 
edges  due  to  the  weather,  is  not  very  significant.  Our  ver- 
tical coefficient  (Singer  and  Smith  "D"),  is  a little  more 
stable  than  the  Gi f ford-Pasqui 1 1 "F"  category  used  by 
Germeles  and  Drake.  [Incidentally,  the  recent  Brookhaven 
data  show  stabilities  that  exceed  Pasquill  F by  a factor  of 
2,  for  example]. 
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I The  Coast  Guard  report  comments  on  our  use  of  mixed 

\ weather  coefficients.  When  we  did  our  tests  there  were 

essentially  no  data  for  mixing  over  water.  If  we  had 
applied  the  meteorological  data  in  the  usual  way  (i.e., 
lapse  rate,  wind  speed,  etc.),  the  weather  during  our 
tests  would  have  generally  been  classified  category  "C", 
slightly  unstable.  But  we  recogni2ed  that  this  was  not 
the  case;  our  plumes  behaved  as  though  the  weather  was 
much  more  stable  than  that.  We  concluded  that  weather 
over  water  was  simply  different  than  over  land,  generally 
much  more  stable.  The  Brookhaven  data*  have  become  avail- 
able since  that  time  and  confirm  the  high  degree  of 
stability.  The  BNL  report  also  points  out  that  application 
of  land-based  weather  correlations  to  predict  stability  over 
water  will  give  large  errors  - just  as  we  had  concluded 
earl i er . 

We  chose  the  spread  coefficients  that  we  did  simply 
because  we  thought  they  matched  our  data  best.  We  would 
not  want  to  claim  any  great  generality  for  them,  [See  p.  72 
of  the  Coast  Guard  Report].  If  we  were  doing  the  work  again, 
we  would  make  use  of  the  data  that  have  since  been  measured 
over  water. 

5.  A few  other  general  comments  on  the  Coast  Guard  Report: 

a.  Table  VI-1  shows  correlations  for  maximum  pool  sizes 
I suggest  that  my  ASME  paper  would  be  a useful  addition. 
Opinions  may  vary  as  to  the  quality  of  that  paper,  but 
it  does  have  one  major  advantage  - it  contains  data, 
all  of  the  'ata  that  were  available  at  the  time  it  was 

• written.  Further,  the  data  cover  the  impressive  range 
of  spill  size  from  5 to  10,000  lbs. 

I'm  not  familiar  with  all  the  references  cited  in 
Table  VI-1,  but  those  that  I recognize  represent  early 
theoretical  studies,  done  before  data  was  available. 

Some  of  the  early  theory  has  been  proven  wrong  by  later 
experimental  work. 

b.  The  SAI  analysis  appears  to  be  a significant  contribution 
but  really  needs  critical  evaluation. 

It  seems  obvious  to  ask  that  the  analysis  should  be 
checked  against  experimental  data.  To  my  knowledge,  it 
has  been  used  to  check  our  Run  11  (relatively  high  wind 
speed,  18  mph),  and  gives  a fair  check  - not  an  exact 
check  bj  any  means.  I would  like  to  see  it  checked 
against  other  data,  particularly  at  lower  wind  speed 
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(e.g.,  our  Run  10  at  5 mph).  Interestingly,  the  SAI 
calculations  show  an  effect  of  wind  speed  that  is 
opposite  to  the  small  scale  experimental  data.  It 
is  important  to  see  if  this  occurs  in  the  analysis  only  * 
for  larger  spills. 

I have  to  say  that  the  SAI  calculated  result  presented 
in  Table  IV-1,  a distance  of  1.2  miles  for  a 37,500  cu.m, 
spill,  seems  highly  improbable  to  me.  It  just  loaks  un- 
reasonable when  plotted  alongside  the  measured  data. 

Their  value  of  3.75  miles  for  a 15  m/s  wind  speed  appears 
much  more  acceptable.  Check  calculations  against  the 
existing  data  would  lend  some  confidence.  Final  answers 
will  probably  await  larger  tests. 
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March  25,  1977 


Dr.  Jerry  R.  Havens 
c/o  U.  S.  Coast  Guard 
Hazardous  Materials  Division 
Room  8308 

Washington,  D.  C.  20590 
Dear  Dr.  Havens: 

Please  find  enclosed  comments  on  your  draft 
report  entitled  "Predictability  of  LNG  Vapor  Dispersion 
from  Catastrophic  Spills  onto  Water:  An  Assessment”, 

January  1977. 

The  comments  were  prepared  in  the  main  by  Drs. 

A.  E.  Germeles  and  F.  Feakes.  We  thank  you  for  the 
opportunity  to  comment. 


Donald  W.  Oakley 
President,  Distrigas 


Corporation 
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COMMENTS  BY  A.  E,  GERMELES  AND  F.  FEAXES  ON  THE  DRAFT  REPORT 
"PREDICTABILITY  OF  LNG  VAPOR  DISPERSION 
FROM  CATASTROPHIC  SPILLS  ONTO  WATER:  AN  ASSESSMENT", 

PREPARED  BY  J.  A.  HAVENS  FOR  U.S.C.G.  (JANUARY  1977) 


We  wish  to  thank  the  Coast  Guard  for  the  opportunity 
to  comment  on  the  draft  report  "Predictability  of  UNG  Vapor 
Dispersion  from  Catastrophic  Spills  onto  Water:  . An  Assessment" 
by  Jerry  A.  Havens.  Dr.  Havens'  contribution  helps  reduce 
the  confusion  that  has  developed  concerning  this  subject. 

It  is  welcomed  because  it  is  clear  that  the  differences 

and  errors  in  models  have  lead  to  overemphasis  of  the  question 

of  "How  far  will  LNG  vapor  travel?" 

We  have  one  major  comment#  and  several  of  a more  technical 
nature.  The  main  comment  relates  to  Tables  IV-1  and  VI-3 
where  the  Cabot  model  based  on  the  work  of  Germeles-Drake 
is  used  with  F weather  to  compute  a downwind  distance  of  about 
10  miles  for  a 25,000  m3  instantaneous  spill.  Cabot's  ex- 
peri  nee  is  clear.  The  U.  S.  Coast  Guard  permits  LNG  ships 
to  enter  Boston  Harbor  only  in  daylight.  The  worst  meteoro- 
logical condition  that  is  reasonably  applicable  du.ving  the 
day  is  D weather.  For  D weather  the  Germeles-Drake  model 
gives  a maximum  downwind  distance  of  about  3 miles.  If  the 
Coast  Guard  maintains  its  present  rules,  we  believe  that  for 
spills  onto  water, D weather  is  the  worst  applicable  stability 
class  and  about  3 miles  is  the  maximum  downwind  travel  distance 

The  other  areas  of  comment  pertain  to: 

• The  sensitivity  of  the  Germeles-Drake  model 
to  chosen  values  of  the  entrainment  constant  a. 


• Comments  on  the  Fay-Lewis  model. 

• Comparison  of  Fay-Lewis  and  Germeles -Drake 
models . 

• Recommendations  for  further  work  on  the  SAI 
model. 

• Further  comments  on  the  choice  of  applicable 
weather  stability  classes. 

% 

These  subjects  are  considered  in  more  detailed  below. 

Entrainment  Constant  Sensitivity.  - On  page  54  and 
in  Figure  VI-5,  Dr.  Havens  presents  the  results  of  parametric 
studies  in  which  he  varied  the  entrainment  constant  a over 
the  range  from  0.01  to  0.5  in  the  Cermeles-Drake  model.  The 
resulting  large  sensitivity  led  him  to  conclude  on  page  117 
that  there  is  too  much  uncertainty  in  the  Germeles-Drake 
model.  It  is  not  reasonable  to  consider  values  of  a as 
large  as  0.5.  There  is  no  known  nonenergetic  entraining 
system  that  entrains  such  large  amounts.  As  pointed  out  in 
Reference  4,  a reasonable  value  for  a is  0.1;  a value  as 
large  as  0.15  might  be  possible,  but  surely  nothing  larger  than 
about  0.2.  The  conclusion  that  there  is  too  much  uncertainty 
in  the  Germeles-Drake  model  is  therefore  not  warranted.  As 
can  be  seen  from  Figure  VI-5,  for  reasonable  values  of  a 
(about  0.1),  the  uncertainty  in  the  values  predicted  by  the 
Germeles-Drake  model  is  relatively  small. 

Comments  for  Fay-Lewis  Model.  - Two  important  parameters 
in  the  atmospheric  dispersion  phase  of  the  model  are  the 
radius  (re)  and  the  height  (he)  of  the  cloud  at  neutral  buoyancy. 


For  a spill  of  25,000  m3 , Dr.  Havens  has  used  the  values 
given  by  the  Lewis  thesis  * re  - 816  m and  hc  - 2.9  in. 

Fay  and  Lewis  claim  that  these  dimensions  are  the  dimensions 
of  the  cloud  at  neutral  buoyancy.  Simple  arithmetic  will 
show  that  a pure  methane  cloud  of  mass  equivalent  to  that 
from  a 25,000  m3  LNG  spill  must  be  at  -259°F  in  order  to 
have  these  dimensions  — and  therefore  the  cloud  is  not 
neutrally  buoyant.  Apparently,  there  is  a basic  physical 
inconsistency  in  the  Fay-Lewis  model.  This  raises  very 
serious  doubts  about  the  values  used  for  re  and  he  and 
about  the  credibility  of  the  entire  analysis. 

Another  fundamental  question  on  the  Fay-Lewis  model 
involves  the  true  asymptotic  behavior  of  Equation  (VI- 5) , 
which  has  been  proposed  by  Fay  and  Lewis  for  calculating  LNG 
vapor  dispersion  in  the  atmosphere.  Fay  and  Lewis  have  claimed 
that,  under  certain  conditions,  this  equation  is  asyr;\otic 
to  Pasquill  dispersion  equations  (namely.  Equations  ;\Li-3)  and 
(VI-4)),  thus  leaving,  perhaps,  the  impression  that  Equation 
(VI-5)  is  not  that  different  from  classical  Pasquill  dispersion. 
However,  the  required  conditions  are  not  met  by  the  UJG  vapor 
dispersion  cases  considered  in  this  work.  This  point  is 
illustrated  by  the  following  table,  which  was  derived  for  the 
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Downwind 
Distance (Km) 

Stability 

oyI(m> 
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VI-3  VI-4 
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27 

Very  Stable 

176 

25.2 

0.99 

0.092 

0.050 

51 

Very  Stable 

310 

37.2 

0.22 

0.062 

0.025 

2.3 

Neutral 

74.3 

33.8 

4.13 

0.068 

0.050 

4.8 

Neutral 

146 

56.6 

0.64 

0.041 

0.025 
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Fay  and  Lewis  state  that  Equation  (VI-5)  is  asymptotic  to 
Equation  (VI-3)  for  oyI  >>  re  and  oZI  >>  he  (large  distances), 
and  to  Equation  (VI-4)  for  o x <<  re  and  oZI  >>  he  (inter- 
mediate distances) . The  above  table  shows  that  these  conditions 
are  not  met.  The  results  from  the  various  equations  differ 
considerably.  The  only  place  where  the  asymptotic  criteria 
seem  to  be  met  is  in  the  third  line,  but  even  here  the  results 
from  Equations  (VI-4)  and  (VI-5)  differ  by  36%  because,  evidently, 
the  criteria  for  intermediate  distances  are  not  satisfied  in 
the  required  sense.  Contrary  to  the  claims  of  Fay  and  Lewis, 
Equation  (VI-5)  is  not  close  to  Pasquill  dispersion  for  LNG- 
cloud  dispersion  calculations  of  practical  interest. 


Dr.  Havens  has  also  shown  the  large  difference  between 
Equation  ( ri-5)  and  classical  Pasquill  dispersion  equations 
(see  Figure  vi-i) , but  does  not  point  out  the  implication 
stated  above,  that  the  Fay-Lewis  dispersion  model  is  not  of 
the  Pasquill  type.  Instead,  Dr.  Havens  states  repeatedly 
that  dispersion  in  the  Fay-Lewis  model  is  based  on  Pasquill 
dispersion  (see  pages  26,  35,  40  and  46).  Still  another 
question  that  might  be  raised  is:  If,  indeed,  dispersion  in 
the  Fay-Lewis  model  is  not  of  Pasquill  type,  then  are  dispersion 
coefficients,  formulated  and  quantified  for  Pasquill  dis- 
persion, applicable  to  non-Pasquill  dispersion  techniques? 

(The  Slade  coefficients  used  by  Fay  and  Lewis  are  Pasquill 
dispersion  coefficients.) 


Comparison  of  Fay-Lewis  and  Germeles-Drake  Models.  - On 
page  56,  Dr.  Havens  states  that  predictions  of  downwind  vapor 
travel  from  the  Fay-Lewis  and  Germeles-Drake  models  "during 
neutral  and  stable  atmospheric  conditions  are  in  close  agreement 
According  to  Tables  VI-3  and  VI-4,  predictions  from  the  two 
models  differ  by  about  40%  to  over  60%  under  all  conditions 
analyzed.  This  is  not  close  agreement.  Further,  in  view 
of  the  several  inconsistencies  in  the  Fay-Lewis  approach,  any 
agreement 'between  the  two  models  is  fortuitous. 

Recommendations  on  SAI  Model.  - The  SAI  approach  to 
calculating  downwind  travel  distances  has  a number  of 
attractive  features  from  a theoretical  point  of  view,  we, 
however,  have  not  been  able  to  either  check  their  estimates 
or  to  ascertain  the  relative  importance  of  the  differences  be- 
tween the  SAI  model  and  the  Germeles-Drake  approach.  We 
recommend  that  efforts  be  made  to  make  a more  direct  com- 
parison than  Dr,  Havens  has  made  and  that  the  SAI  mode]  be 
tested  for  its  sensitivity  to  important  parameters.  It 
would  be  of  great  interest  to  us  if  the  SAI  model  confirmed 
that  our  estimate  of  about  3 miles  for  maximum  downwind  dis- 
tance is,  in  fact,  conservative. 

Applicable  Weather  Stabilities.  - For  his  comparisons. 

Dr.  Havens  uses  stability  classes  D (neutral)  and  F (most 
stable)  with  the  Germeles-Drake  model  and  computes  maximum 
downwind  travel  distances  of  about  3 and  10  miles,  respectively, 
for  5%  average  concentrations.  Dr.  Havens  states  on  page  13 
that  he  used  "the  "worst  applicable"  meteorological  conditions 
suggested  by  the  groups".  Lest  the  impression  is  left  that 


-6- 


flammable  concentrations  might  travel  downwind  to  10  miles 
and  more,  we  must  emphasize  that  we  consider  calculations 
based  on  class  F weather  as  academic  and  of  little  practical 
significance.  In  our  opinion,  the  "worst  applicable"  class 
is  D and,  therefore,  about  3 miles  is  an  upper  bound  for  the 
downwind  travel  of  5%  average  concentrations.  Our  reasons  for 
considering  class  F (and  even  class  E)  as  inapplicable  are  as 
follows : 

(i)  According  to  Turner  (Reference  27,  page  6)  classes 
E and  F are  possible  only  during  nighttime, 

(ii)  Current  Coast  Guard  regulations  require  that  UJG 
tankers  come  into  port  only  during  daytime. 

(iii)  According  to  Turner  (Ibid) , the  standard  Gifford 
Pasquill  classes  have  been  defined  for  "open 
country  or  rural  areas".  It  iB  important  to 
keep  in  mind  that  in  calculating  LNG  vapor  dis- 
persion as  a part  of  safety  analyses  for  metropolitan 
areas,  more  unstable  classes  should  be  used  because 
of  "the  larger  surface  roughness  and  heat  island 
effects"  of  such  areas. 


From  a practical  point  of  view,  Cabot  does  not  regard 
large  vapor  cloud  travel  distances  as  a reasonable  possibility. 
The  conditions  specified  by  Dr.  Havens  on  page  15  of  the 
draft  are  extremely  unlikely,  if  not  impossible,  if  one  takes 
into  account  the  strict  Coast  Guard  rules  that  have  been 
applied  in  Boston  Harbor.  Massive  spills  from  LUG  tankers 


would  require  a highly  energetic  collision  with  a large  ship 
and  we  believe  that  this  possibility  is  eliminated  by  Coast 
Guard  rules  creating  a traffic-free  harbor*  Under  these  cir- 
cumstances , the  estimation  of  vapor  travel  distances  is  an 
interesting  mathematical  exercise  done  in  response  to  the 
National  Environmental  Policy  Act. 

The  ^results  of  these  vapor  travel  calculations!  however, 
need  to  be  placed  in  the  proper  perspective  as  one  element  of 
a careful  risk  analysis.  Risk  is  a function  of  both  the 
probability  and  the  consequences  of  an  undesired  event.  The 
relevant  probability  includes  early  ignition  as  part  of  an 
event  leading  to  a large  spill  as  well  as  the  likelihood  of 
ignition  by  land-based  and  water-based  sources.  And  the 
resulting  risks  can  only  be  assessed  on  a realistic  basis  if 
the  fire  hazards  are  compared  with  other  flammable,' although 
less  volatile,  fu_l  substitutes  for  XJNG. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
DEPARTMENT  OF  MECHANICAL  ENGINEERING 
CAMBRIDGE  MASSACHUSETTS  02159 


JAMES  A FAY 
nortttot 


March  10,  1977 


Dr.  Jerry  Havens 

Office  of  Merchant  Marine  Safety 
U.S.  Coast  Guard 
400  7th  St. , S.W. 

Washington,  D.C. 

Dear  Dr.  Havens: 

Enclosed  are  my  comments  on  your  draft  report.  At  the  very  end 
are  two  detailed  comments  which  you  may  wish  to  delete  if  you  make 
corrections  to  the  necessary  parts  of  the  draft  paper. 

Thanks  for  the  courtesy  of  asking  for  my  comments. 


Sincerely  yours, 


JAF:daf 

cc:  H.  Walter 

enc. 
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COMMENTS  ON  "PREDICTABILITY  OF  LNG  VAPOR  DISPERSION 
FROM  CATASTROPHIC  SPILLS  ONTO  WATER:  AN  ASSESSMENT" 

James  A-  Fay 

Massachusetts  Institute  of  Technology 

This  excellent  review  serves  two  very  important  purposes.  It 
compares  the  predictions  of  the  several  LNG  vapor  cloud  dispersion  theories 
on  a common  basis,  i.e.,  a given  spill  size  and  distances  to  given  ground 
level  concentrations.  More  importantly,  it  explains  in  full  detail  the 
various  assumptions  and  calculation  procedures  used  in  each  theory 
which  are  not  always  adequately  described  in  the  original  publications. 

The  disparate  predictions  of  the  various  approaches  are  well  illustrated 
and  some  of  the  intermediate  steps  in  determining  the  downwind  concen- 
trations are  usefully  contrasted.  It  will  be  very  helpful  in  clarifying 
the  state  of  knowledge  regarding  vapor  cloud  dispersion  and  suggesting 
further  analytical  and  experimental  approaches  to  a more  reliable  method 
of  prediction. 

The  discussion  in  the  introduction  (pp.  15-17J  of  the  probability 
of  various  accident  scenarios,  which  is  clearly  not  an  aspect  of  the 
scientific  review  of  the  various  dispersion  theories  but  more  nearly 
a policy  statement  regarding  risk,  unfortunately  tends  to  denigrate  the 
value  of  this  analysis.  The  reader  may  wonder  whether  *-he  assessment 
is  to  be  taken  Beriously,  or  has  been  carefully  made,  gxvori  the  asserted 
unlikelyhood  of  the  process  being  discussed.  But  if  one  ignores  the 
casuistry  of  this  portion  of  the  introduction,  the  subsequent  analysis 
is  scientifically  useful  and  more  than  worth  the  effort  to  have  performed 


it. 
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On  p.  35  (and  inferentially  in  Table  VI-2)  it  is  stated  that  the 

instantaneous  spill  models  assume  that  the  vapor  cloud  radius  and  the 

✓ 

liquid  pool  radius  are  equal  at  the  end  of  the  evaporation  period. 
Neither  Fay  (24)  nor  Fay  and  Lewis  (6)  make  this  assumption,  nor  is  it 
necessary  for  the  determination  of  vapor  cloud  spread  according  to  their 
analysis,  but  this  assumption  is  used  by  Germeles  and  Drake  (4)  (see 
Fig.  VI-3)  and  the  FPC  staff  (7)  (see  p.  80)  as  an  intermediate  step 
in  determining  vapor  cloud  spread.  I expect  that  the  vapor  cloud  will 
extend  beyond'  the  edge  of  the  pool  at  the  end  of  the  evaporation  period, 
but  this  point  obviously  deserves  further  investigation. 

In  determining  the  molar  concentration  of  LNG  according  to  the 
method  of  Fay  and  Lewis  (6) , the  instantaneous  source  strength  Q in 
Eq.  (VI-2)  and  the  equation  following  it  should  be  the  volume  of  pure 
vapor  at  atmospheric  temperature,  i.e.,  590  times  the  liquid  volume 
or  2.45  times  the  saturated  vapor  volume  at  atmospheric  pressure  when 
the  air  temperature  is  0°C.  Since  the  Gaussian  puff  vapor  dispersion 
equations  conserve  the  partial  volume  of  the  dispersing  material  and 
hence  assume  a constant  temperature  dispersal  process  for  the  gas  con- 
taminant being  dispersed,  the  equivalent  source  strength  should  be  the 
constant  partial  volume  of  vapor  at  atmospheric  temperature.  Thus  the 


initial  vapor  cloud  height  hy  in  Eq.  (VI-3)  and  subsequent  equations 

should  be  determined  from  this  instantaneous  source  strength  according 
2 

to  h * Q/TTr 
v * vm 

For  the  purpose  of  determining  the  maximum  radius  of  the  vapor 
cloud  at  the  point  of  neutral  buoyancy.  Fay  (24)  assumed  that  the  vapor 
cloud  motion  would  be  the  same  as  that  of  an  adiabatic  cloud  of  saturated 


vapor  equivalent  to  the  spill  volume  and  spreading  for  a time  needed 
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to  transfer  heat  sufficient  to  render  it  neutrally  buoyant.  The  height 
of  this  adiabatic  cloud  at  the  end  of  the  cloud  spread  period  should 
not  be  used  as  the  starting  point  for  the  dispersion  calculation  for 
the  reasons  given  above. 

The  comparison  shown  in  Fig.  VI -1  is  very  useful  in  that  it  shows 
that  the  Fay  and  lewis  (6)  dispersion  model  cannot  be  adequately  rep- 
resented by  a point  source  model  even  at  the  largest  distances  of  interest. 

This  is  a consequence  of  the  effect  of  gravity  spread  of  the  vapor 
cloud  in  the  early  stages  which  produces  an  initial  shape  quite  differ- 
ent from  that  which  eventually  ensues  far  downstream.  It  is  also  apparent 
from  a close  examination  of  this  figure  that  a "virtual  source"  solution 
of  the  form  of  Eq.  (VI-9)  will  also  fail  to  match  the  Fay  and  Lewis 
solution  over  most  of  the  region  of  interest,  for  the  same  reason. 

The  inability  of  virtual  source  models  to  account  for  the  initial 
cloud  shape  is  well  illustrated  by  the  Germeles  and  Drake  (4)  solution 
for  neutral  stability  (p.  52,54).  At  the  actual  source,  where  the 
gravity  spread  model  concentration  is  matched  to  the  virtual  source  puff 
model  concentration,  the  gravity  cloud  radius  is  750  m and  the  height 
is  18.4  m (Fig.  vi-3) . In  contrast,  the  virtual  source  model  horizontal 
and  vertical  deviations  at  this  same  point  can  be  found  from  Appendix  I 
to  be  280  m and  82  m respectively.  Thus  the  puff  model  aspect  ratio 
(width/height)  is  nearly  one  twelfth  that  of  the  calculated  cloud  at 
the  point  where  the  former  is  supposed  to  depict  the  beginning  of  the  dis 
persion  process.  It  would  seem  that  the  virtual  source  models  are 
inappropriate  for  describing  dispersion  ol  clouds  of  such  unusual  shape. 
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illustrated  by  the  CHRIS  model  (p.  77) . The  concentration  calculated 
from  Eq.  (VII1-15)  evaluated  at  the  actual  source,  for  the  case  of  F 
stability  (for  which  5%  concentration  is  reached  at  86,000  m)  is  2226% l 
Thus  this  model  predicts  concentrations  exceeding  100%  for  very  long 
distances  from  the  spill. 

In  addition,  the  use  in  CHRIS  of  a steady  source  model  for  a source 
of  finite  lifetime  leads  to  quite  peculiar  results.  Since  streamwise 

t 

diffusivity  does  not  enter  the  calculation,  the  streamwise  cloud  dimen- 
sion would  approximately  equal  the  wind  speed  times  the  evaporation  • 
time  or  about  500  m.  But  for  F stability,  the  plume  transverse  and  ver- 
tical deviations  at  5%  concentration  are  19000  m and  90  m respectively. 
Such  an  odd-shaped  cloud,  with  the  transverse  dimension  4C  times  the 
streamwise  dimension,  does  not  seem  consonant  with  known  dispersion 
characteristics. 

dne  calculation  according  to  the  Germeles-Drake  model  of  the  effect 
of  different  entrainment  rates  on  vapor  concentration  during  the  gravity 
spreading  phase  of  the  vapor  cloud  motion,  as  depicted  in  Fig.  VI-5, 
clearly  indicates  the  significance  of  assumptions  regarding  the  magnitude 
of  this  process,  as  the  author  emphasizes  on  p.  117.  These  assumptions 
are  equally  important  to  the  SA1  model.  In  my  opinion,  the  very  rapid 
dilution  calculated  by  SA1  is  directly  related  to  their  assumed  (and  pre- 
sumably high)  values  of  vertical  momentum  diffusivity. 

Entrainment  coefficients  rarely  exceed  0.1,  and  then  only  for  mixing 
processes  across  gravitationally  unstable  interfaces.  Indeed,  the 
observation  that  intrusions  exist  for  layers  having  very  large  values 
for  the  ratio  of  horizontal  to  vertical  dimensions  indicates  that 


entrainment  coefficients  must  be  very  small  for  gravitationally  induced 
motion  of  this  type.  While  the  parametric  study  of  the  effect  of  various 
entrainment  coefficients  on  gravitational  spread  is  a useful  analytical 
tool,  it  is  doubtful  that  the  calculations  for  high  entrainment  coeffi" 
cients  are  describing  physically  realizable  processes. 

The  vertical  momentum  diffusivity  used  in  the  SA1  model  appears 
to  affect  the  results  significantly  since  j.t  determines  the  entrainment 
rate  during  gravitational  spread.  The  explanation  given  on  p.  108 
of  the  choice  of  stability  parameter  (which  affects  the  choice  of  diff- 
usivity) is  not  sufficient  to  enable  a reader  to  reproduce  the  SA1  pre- 
scription. A more  precise  explanation  is  very  desirable. 

Important  information  for  the  FPC  model  appears  to  be  lacking. 

The  heat  transfer  coefficient  h used  in  Eq.  (VII-25)  is  not  specified 
nor  is  it  explained  how  it  is  to  be  determined.  Similarly,  the  origin 
of  the  heat  transfer  coefficient  used  in  Eq.  (VII-28)  should  also  be 
explained.  It  would  also  be  important  to  obtain  from  the  FPC  staff 
a physical  explanation  of  the  vapor  release  process  calculated  on  p.  84 
if  an y serious  consideration  is  to  be  given  to  this  model. 

It  would  be  more  accurate  to  describe  the  dispersion  models  (p.  113) 
as  including  the  effects  of  gravitational  spread  as  a precursor  to  neu- 
tral buoyancy  dispersion  or  as  the  determinant  of  the  initial  conditions 
for  the  latter.  The  model  may  or  may  not  include  entrainment  during 
the  spreading  process,  but  if  it  does  the  mixing  .is  related  to  spreading 
speeds  and  not  to  atmospheric  turbulence.  These  models  also  conserve 
mass,  momentum  and  energy  (to  various  approximations)  as  does  the  diff- 
erential equation  (SA1)  approach. 


The  author's  suggestion  that  heat  and  momentum  transfer  effects  on 
the  vapor  dispersion  from  small  scale  spills  is  unimportant  (p.  114-116) 
may  be  true,  but  the  evidence  in  support  of  it  given  by  the  author  is 
far  from  convincing.  First  of  all,  it  is  a confined  land  spill  (in  con- 
trast to  the  water  spills  exclusively  treated  in  this  paper)  which  is 
considered.  Secondly,  given  the  kind  of  disagreement  between  the  models 
for  a large  water  spill,  all  of  which  (including  Burgess'  model)  in- 
clude gravity  effects  to  a greater  or  lesser  degree,  the  comparison 
in  Fig.  IX-1  is  probably  fortuitous.  But  certainly  this  is  a matter 
deserving  further  thought  and  amalysis. 

The  advantages  of  the  differential  equation  model,  such  as  SA1 
model,  are  not  so  one-sided  as  the  author  suggests  on  p.  120.  For 
example,  such  models  will  not  predict  the  observed  dispersion  in  homo- 
geneous turbulent  flow.  But  since  the  vapor  cloud  is  being  dispersed 
in  the  atmospheric  shear  layer,  the  approach  of  relating  the  local 
diffusivity  to  the  distance  from  the  surface  and  the  local  gravitational 
stability  parameter  may  be  a reasonable  approximation.  Nevertheless,  it 
would  be  very  desirable  to  compare  such  solutions  with  measurements  of 
dispersion  of  passive  trace  diluents.  Also,  there  are  other  practical 
disadvantages  to  such  models,  for  example,  expense  of  obtaining  solutions 
and  hence  testing  for  the  sensitivity  to  various  assumptions. 

In  summary,  these  comments  are  made  to  elaborate  and  develop  several 
of  the  points  raised  by  the  author  and  thus  to  improve  the  general 
level  of  understanding  of  this  difficult  problem. 
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DET AILED  COMMENTS 


On  p.  41,  the  argument  of  the  exponential  term  in  Eq.  (rv-2)  should 
2 2 

read  (-z  /2CZ  ).  This  follows  from  Eq.  (1)  of  Fay  and  Lewis  (6)  for 
the  case  of  r * 0. 

For  the  reasons  explained  above,  regarding  the  determination  of 
hv,  the  height  in  item  II,  Tables  VI-3  and  VI-4  should  read  7.1  m under 
column  one.  The  corresponding  distances  in  items  IV  and  V of  column 
1 should  be  28.0  miles  and  47.2  miles  in  Table  VI-3  and  3.0  miles  and 
5.3  miles  in  Table  VI-4.  Line  5 in  Table  IV-1  should  thus  read  28.0 
miles.  Also,  Fig.  VI-1  should  be  modified  accordingly. 
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KEDERAL  POWER  COMMISSION 
Washington,  D.C,  20426 

March  3,  1977 


Mr.  Jerry  A.  Havens 

Cargo  & Hazardous  Materials  Div. 

U.  S.  Coast  Guard  (G-MHM/83) 

Washington,  D.  C.  20590 

pear  Sir: 

I appreciate  the  opportunity  to  review  and  comment 
on  your  draft  report  "Predictability  of  LNG  Dispersion  from 
Catastrophic  Spills  into  Water  --  An  Assessment."  It  is 
useful  to  have  a review  of  these  dispersion  models  under 
one  cover,  and  you  have  done  a good  job  of  placing  them 
in  perspective.  I agree  completely  tnat  "The  SAI  model  . . . 
is  a significant  advance  ..."  (page  119).  I would  hope 
that  in  the  future  no  U.  S.  Government  agency  would  support 
any  more  work  or  spend  any  more  staff  effort  on  development 
of  any  LNG  dispersion  model  less  adequate  or  complete  than 
this . 

In  general,  the  report  appears  to  be  heavy  on  providing 
everyone's  equations,  but  light  on  why  the  models  are  in- 
adequate ("assessment"  is  in  the  title).  Thus,  given  the 
evident  quality  of  the  SAI  model  and  its  reasonable  limits 
for  downwind  vapor  travel,  I would  suggest  a concise  summary 
of  reasons  why  those  models  that  produce  much  longer  plumes  -- 
by  factors  of  10  tc  AO  times  too  much  (page  1A  and  page  3A)  — 
are  so  erroneous.  Such  a summary  of  how  each  model  is  deficient 
compared  to  SAI  would  be  helpful,  particularly  in  FPC  cases 
involving  LNG  applications.  These  models  have  caused  con- 
siderable confusion  and  delay  in  hearings.  Having  such  material 
available  before  hand  could  marked] y shorten  the  hearing  process. 

It  is  most  unfortunate  that  vou  have  placed  so  much 
emphasis  on  LNG  spills  of  25,000  m3,  which  is  a size  that 
is  probably  too  large  ever  to  be  observed.  The  following 
calculation  illustrates  this  point. 
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The  probability  of  a water  spill  from  an  LNG  ship 
accident  is  the  product  of  the  probability  of  an  accident 
times  the  probability  of  a spill  given  an  accident.  A way 
of  estimating  the  accident  probability  is  to  consider  LNG 
tanker  operating  experience.  I have  estimated  recently 
that  since  1964  there  have  been  about  1,600  LNG  tanker  voyages 
worldwide,  or  3,200  tanker  transits  of  ports,  harbors,  or 
piers  while  loaded  with  cargo,  without  a major  accident  or 
major  spill.  From  this  observed  excellent  accident  rate  an 
estimate  of  the  true  accident  rate  may  be  made  using  standard 
statistical  techniques.  The  result  is  about: 

1.5  x 10"3  accidents /transit 
for  all  types  of  accidents . 

In  order  to  get  a feeling  of  the  probability  of  a 
25,000  m3  spill  from  a tanker  accident  consider  that  this 
is  equivalent  to  about  6.6  x 10^  gallons  of  oil  in  volume. 

The  probability  of  such  an  oil  spill  in  U.  S.  ports,  harbors, 
or  piers,  based  on  data  from  the  Oceanographic  Institute  of 
Washington,  1974,  is  about 

-25.000  -16 

e 1,5/0  “ e of  1.1  x Hr  spills  per 

accident 

This  is  discussed  more  fully  in  the  FPC  Final  Environmental 
Impact  Statement  on  Alaska  Natural  Gas  Transportation  Systems, 
April  1976,  Vol.  Ill,  page  410-413.  This  important  reference 
does  not  seem  to  have  been  included  in  your  list.  Thus  the 
probability  of  such  a large  spill  is  about  1.7  x 10-1*  per 
transit,  which  is  indeed  negligible  for  any  foreseeable  annual 
rate  of  LNG  deliveries  worldwide. 

This  view  is  underscored  by  the  FPC  Administrative  Law 
Judge  in  his  "Initial  Decision  on  Proposed  Alaska  Natural 
Gas  Transportation  Systems."  FPC,  February  1977,  in  which  he 
states  in  part:  (page  94)  "The  fear  raised  by" those  opposing 
LNG  facilities  in  populated  areas  requires,  therefore,  certain 
assumptions.  First,  there  must  be  a large  spill."  (Page  95) 
"In  order  to  achieve  the  large  size  vapor  cloud  necessary  to 
create  even  measurable  risks  for  people  located  some  distances 
away,  an  assumption  has  to  be  made  that  a high  volume  of  LNG 
be  released  instantaneously  ..."  (page  95)  "LNG  is  hazardous 
and  must  be  treated  with  respect.  The  risks  associated  with 
its  use  must  be  analyzed.  But,  they  must  be  done  so  on  a 
credible  basis  with  assumptions  that  are  in  themselves 
credible,  and  much  of  the  risk  analysis  has  not  been  done  on 
that  basis."  Your  otherwise  fine  report  may  therefore  permit 
misleading  information  because  it  analyzes  essentially  im- 
possible events. 
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Obviously,  plumes  from  small  water  spills  of  LNG  (6100  m^), 
while  presumably  more  prevalent,  do  not  represent  a public 
hazard  either.  The  range  of  interest  encompasses  spills  that 
are  large  enough  to  be  a hazard  but  small  enough  to  possibly 
occur.  Based  on  probability  studies  at  FPC  (in  the  first 
reference  above),  this  range  is  believed  to  be  500-3,000  nH. 

What  are  the  comparative  results  from  the  models  you  have 
analyzed  for  this  spill-size  range? 

I note  on  page  116  that  “sufficient  accuracy  may  be 
obtainable  from  classical  dispersion  models  for  the  prediction 
of  LNG  vapor  dispersion  from  small  spills  on  land  or  water," 
where  "small"  is  not  otherwise  defined.  It  would  be  helpful 
to  know  the  accuracy  expected  from  these  models  compared  to 
the  SAI  model  in  the  above  spill-size  range.  The  loss  in 
accuracy  may  be  more  than  offset  by  the  substantially  lower 
cost  incurred  from  their  use.  Likewise,  it  would  a pear  to 
be  worthwhile  to  perform  sensitivity  analyses  on  the  SAI 
model  in  this  spill-size  range  in  order  to  reduce  it's  com- 
putation cost  without  sacrificing  significant  accuracy.  Such 
sensitivity  analyses  should  be  supported. 

The  discussion  and  rationale  on  page  15  for  performing 
this  assessment  of  downwind  vapor  plumes,  in  spite  of  the 
probability  of  prompt  ignition  of  the  LNG  vapor  being 
"extremely  likely,"  sterns  shaky.  The  truth  is  that  the 
available  accident  reports  from  the  Coast  Guard  show  that 
the  probability  of  prompt  ignition  after  oil  tanker  accidents 
is  not  known  accurately,  but  appears  to  be  reasonably  large  — 
perhaps  as  large  as  90  percent.  Many  witnesses  at  FPC  hearings, 
including  the  writer,  believe  that  this  is  probably  true  for 
LNG  tanker  collisions  also  (but  here  again  there  are  no  data). 
Extrapolating  to  LNG  tankers,  a 10  percent  chance  of  non- 
ignition with  a subsequent  hazardous  plume  is  sufficient 
cause  for  your  assessment,  I should  think,  without  confusing 
the  reader  further.  In  short,  I do  not  equate  "highly  un- 
likely" with  a 10  percent  probability. 

Let  me  compliment  you  on  a thorough,  clear,  and  timely 
report.  It  is  most  important  that  someone  from  outside  the 
LNG  community  perform  such  an  assessment  at  this  time. 


Environmental  Specialist 


February  16,  1977 
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Dr.  Jerry  Havens 
Technical  Advisor 
Hazardous  Materials  Division 
United  States  Coast  Guard 
Washington,  D.C.  20590 

Dear  Jerry, 

I want  to  thank  you  for  the  extra  effort  it  took  in  seeing  that 
I got  a copy  of  your  draft  report  on  the  "Predictability  of  LNG 
Vapor  Dispersion".  We  still  don't  know  what  happened  to  the 
first  copy  that  was  sent. 

I want  to  compliment  you  on  the  very  excellent  report  which  you 
have  prepared.  I certainly  feel  it  was  an  extremely  worthwhile 
effort  for  someone  who  was  not  directly  involved  in  the  LNG 
community  to  evaluate  the  various  models  which  have  been  used 
to  produce  the  wide  range  of  answers  often  quoted  or  misquoted 
by  non-technical  individuals. 

We  have  a very  limited  number  of  comments  which  you  might 
consider  incorporating  into  your  report.  They  are  as  follows: 

Page  13 

In  the  final  line  of  the  first  full  paragraph,  you  note 
properly  that  the  worst  applicable  meteorological  condition 
suggested  by  the  groups  is  not  necessarily  the  worst  that 
might  have  been  assumed.  I wonder  if  it  might  not  be  worth- 
while siting  the  fact  that  the  SAI  result  quoted,  for 
example,  is  far  from  its  own  predicted  worst  case  associ- 
ated with  a high  wind  condition. 

Page  15 

In  the  sixth  and  seventh  lines  of  the  first  paragraph,  the 
point  is  made  that  ignition  will  probably  occur  in  a high 
energy  collision  because  of  frictional  heating  anticipated. 
While  we  agree  that  this  could  play  a role,  it  is  our  per- 
sonal feeling  that  sparks  and/or  broken  electrical  lines 
or  connections  would  produce  an  even  more  reliable  ignition 
source  than  that  associated  with  frictional  heating. 
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Page  21 

The  first  accumulation  term  should  contain  a p,  i.e.,  • 


Page  95 

Under  the  section  entitled  "Accountability  of  Momentum", 
the  statement  is  made  that  equation  VIII-2  is  expanded  with 
vertical  accelerations  and  viscous  forces  neglected  in  the 
equation  for  accountability  of  vertical  momentum.  Another 
view  of  this  would  be  to  say  that  the  vertical  accelerations 
and  the  viscous  forces  are  assumed  equal  and  opposite.  This 
results  in  the  same  equation  VIII-5c,  but  does  not  give  the 
sometimes  mistaken  impression  that  vertical  velocities  are 
set  equal  to  0,  which  in  fact,  as  you  realize  from  the  other 
equations,  are  not. 

Pages  102-103 

There  is  a spurious  p in  the  horizontal  diffusion  terms 
of  the  conservation  equations  for  u,^,H,^,i.e. 

3 3 

"3x  KH  3x  and  similarly  for  the  y diffusion  terms. 


A definition  of  4>  should  also  be  provided  as  » gz 
(the  geopotential  height). 

Also,  the  substantive  derivative  ^ takes  on  a new  form 
in  the  a coordinate  system  as 


D 

Dt 


On  page  103,  the  p equation  should  contain  an  RT  factor 
in  the  denominator. 


Page  104 

Modifications  to  this  page  are  shown  in  the  attached  copy. 
Page  109 

The  final  sentence  states  that  the  ratios  Table  VIII-3 
are  based  on  proprietary  field  data  obtained  by  SAI.  This 
is  correct , as  stated  in  any  of  our  reports  completed  for 
Western  LNG  Terminal  Company.  However,  it  may  also  be 
stated  that  they  were  compared  with  similar  data  published 
by  Lantz,  and  the  SAI  results  are  conservative,  i.e.,  they 
would  produce  smaller  dif fusivities  in  the  horizontal  direc- 
tion than  those  using  the  results  of  Lantz. 
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Page  111 

In  the  third  line  of  the  first  full  paragraph,  it  is  stated 
that  SAI  has  not  published  their  calculated  results  for  a 
25,000  m3  spill.  I think  it  is  more  properly  stated  that 
SAI  has  not  calculated  the  dispersion  associated  with  a 
25,000  m3  instantaneous  spill. 

For  your  additional  information,  I have  enclosed  a response, 
which  we  prepared  for  Western  LNG  Terminal  Company,  to  a question 
from  the  FPC  regarding  more  detailed  information  on  the  numerical 
methods  used  in  the  SIGMET  code.  I trust  it  will  be  of  some 
value  to  you . 

Again,  let  me  congratulate  you  on  a very  excellent  report. 
Sincerely , 

SCIENCE  APPLICATIONS,  INC. 

Walter  G.  England 
Manager 

Environmental  Sciences  and 
Safety  Division 
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end 
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March  4,  1977 


• Dr.  Jerry  Havens 

Technical  Advisor 

Cargo  and  Hazardous  Materials  Division 
Office  of  Merchant  Marine  Safety 
United  States  Coast  Guard 
Washington,  DC  20590 

Dear  Dr.  Havens: 

Thank  you  for  the  opportunity  to  review  your  comp rehens ive  assessment 
of  the  various  models  in  use  for  prediction  of  vapor  dispersion  from 
potential  LNG  spills  on  water.  Unfortunately,  1 did  not  have  time  to 
make  a detailed  review  of  the  report,  so  my  comments  are  primarily 
based  on  impressions  from  a once-through  reading  of  the  report  and  on 
my  own  concerns  about  the  strengths  and  weaknesses  of  the  various 
models . 

My  colleague,  Dr.  Germeles  of  Cabot  Corporation,  has  discussed  with 
me  and  sent  me  a copy  of  the  comments  he  is  sending  you  based  on  his 
quite  thorough  review  of  your  analysis  of  our  model  and  the  Fay -Lewis 
model.  I concur  with  essentially  all  of  his  remarks. 

My  chief  concern  is  that  the  draft  report  may  give  the  impression  to  a 
reader  who  is  not  thoroughly  versed  in  the  technical  Issues  that  you 
are  recommending  the  SA1  model  as  the  best  available.  While  I concur 
that  their  approach  is  the  most  rigorous,  I strongly  share  the  concerns 
you  express  later  in  your  report  that  the  SAI  model  neede  much  additional 
checking,  sensitivity  testing  and  verification  before  it  should  become  a 
recommended  method.  The  model  developed  by  Dr.  Germeles  and  myself  is 
simplified  to  the  point  of  Including  physically  unrealistic  assumptions; 
however,  these  assumptions  can  be  defended  as  being  conservative.  While 
we  would  expect  our  model  to  overpredict  downwind  hazard  distances,  I 
doubt  that  one  could  sort  through  the  many  assumptions  Incorporated  in 
the  SAI  model  and  say  whether  the  net  effect  is  conservative  or  optimis- 
tic. (This  gets  back  to  your  concerns  about  the  sensitivity  of  the  model 
to  key  parametric  assumptions.) 
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Dr.  Jerry  Havens 
United  States  Coast  Guard 


My  major  concerns  with  the  SAI  model  are  in  the  following  areas: 

1.  Since  the  SAI  model  is  proprietary  and  very  expensive  to  execute, 
it  has  not  yet  been  extensively  studied  by  an  independent  expert  like 
yourself. 

2.  The  SAI  model  should  be  as  sensitive  to  choice  of  mixing  param- 
eters as  the  GD  model  since  the  same  basic  phenomena  are  Involved. 

3.  The  SAX  model  check  with  "data"  are  net  a real  verification  since 
the  small  land  spill  test  used  in  the  comparison  was  not  large  enough  to 
have  discernible  gravity  spreading  behavior.  In  fact,  a simple  line 
source  Gaussian  model  also  is  in  good  agreement  with  the  data. 

4.  Any  large  computer  programs  are  difficult  to  verify  since  they 
may  contain  insidious  errors  or  be  subject  to  subtle  numerical  insta- 
bilities. Only  by  extensive  sensitivity  testing,  comparison  with  analyti- 
cal solutions  for  simple  test  cases,  comparison  with  any  pertlne  t experi- 
mental data  available  (e.g.  gravity  spreading  from  small  spllle  under  calm 
wind  conditions) , and  careful  selection  of  values  (and  uncertainty  bands) 
for  important  parameters  can  one  gradually  build  confidence  in  a complex 
computer  model. 

5.  The  turbulent  mixing  parameters  used  by  SAI  are  based  on  their 
own  data.  When  these  are  compared  with  the  widely  used  Pasqulll-Glfford 
coefficients,  it  appears  that  the  SAI  parameters  themselves  may  partially 
be  responsible  for  minimizing  the  effects  of  atmospheric  stability  on 
downwind  travel. 

6.  The  increasing  downwind  travel  distances  with  increasing  wind 
speed  are  physically  possible  but  have  never  been  observed  in  practice. 
Whether  or  not  this  is  real  could  probably  be  shown  only  by  a series  of 
extremely  large  (and  costly)  experiments. 

Given  these  uncertainties,  I would  have  preferred  you  to  emphasize  that 
1)  the  G-D  model  gives  simplified,  but  conservative  estimates  of  downvjnd 
travel  and  2)  that  the  SAI  approach  1b  an  attempt  to  obtain  a more  physi- 
cally realistic  answer  but  that  the  model  Itself  still  requires  further 
testing  and  scrutiny  before  it  can  be  recomended  per  ee. 

The  ironic  part  is  that  ail  these  models  are  being  developed  for  use  in 
risk  assessment  studies  or  for  definition  of  some  maximum  accident  scen- 
ario. In  fact,  the  present  ranges  of  uncertainties  in  the  models  are  not 
large  compared  to  uncertainties  in  defining  spill  scenarios  (quantity  and 
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rate)  or  in  the  distributions  of  population  and  Ignition  sources  in  the 
, path  of  a hypothetical  vapor  cloud.  In  the  SAI  risk  assessment,  it  is 

their  conservative  estimate  of  ignition  source  distribution  — not  their 
vapor  cloud  analyses  — - that  determine  the  potential  hazards.  With  un- 
• certainties  such  as  these  in  other  assumptions  about  accident  scenarios, 

perfection  of  LNG  vapor  dispersion  models  seems  to  me  to  be  more  of  a 
technically  interesting  goal  than  an  urgently  needed  effort.  (Even  then, 
for  risk  studies,  cloud  width  is  a much  more  important  parameter  than 
maximum  downwind  travel.) 

I'm  enclosing  a few  quick  calculations  using  the  CHRIS  model  for  an  in- 
stantaneous 25,000  spill  which  are  slightly  different  from  those  in 
your  report. 

Please  phone  me  if  you'd  like  to  discuss  any  of  these  points  further. 

Your  report  is  an  excellent  contribution  and  will  be  widely  disseminated, 
so  I'm  sure  we  all  would  like  to  see  it  as  fair  and  easily  understandable 
as  possible. 

With  best  regards. 


Sincerely , 
Elisabeth  M.  Drake 


EMD:km 

cc:  Dr.  A.  Germeles /Cabot  Corporation 

D.  S.  Allan/Arthur  D.  Little,  Inc. 


